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ABSTRACT 


Solid  rocket  motor  nozzle  environment  was  studied  in  cold  flow  modeling 
and  in  motor  firings.  The  emphasis  in  the  work  was  on  the  effects  of  the 
condensed  phase  in  well  impingement  and  radiation.  In  motor  firings  using 
H2  -  O2  -  AI2O3  ~  H2O,  emissivity  values  for  the  alumina  wero  determined.  Then 
radiation  heat  flux  scaling  relationships  were  developed  for  application  to  any 
size  motor.  Particle  impingement  was  studied  in  firing  witi  a  molybdenum  nozzle 
fitted  with  thermocouples  from  which  local  transient  heat  flux  data  were  developed 
High  heat  flux  rates  due  to  impingement  were  determined.  In  cold  flow  studies 
particle  impingement  rates  and  static  pressure  values  for  input  to  corrective 
analysis  were  determined.  Analyses  of  heat  flux  modes  in  a  large  motor  nozzle  and 
in  a  gimballed  nozzle  were  determined.  A  large  part  of  the  work  involved  techniqu 
development  in  handling  micron- s  ized  particles  in  firings  and  in  cold  flow.  A 
particle  sampler  was  developed  for  sampling  motor  exhausts. 
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SECTION  1 


INTRODUCTION 


This  report  sumnarlzes  the  accomplishments  of  a  1-year  study  of  the 
nozzle  environment  of  solid  rocket  motors.  The  work  was  sponsored  by  the  Rocket 
Research  Laboratories,  Edwards  Air  Force  Base,  under  the  technical  cognizance  of 
Mr.  Leonard  Siegers,  Emphasis  in  the  work  was  aimed  toward  developing  new, 
practical  understanding  of  the  environmental  effects  of  the  condensed  phase. 

The  work  performed  on  this  contract  provides  much  useful  information 
which  the  rocket  industry  can  apply  in  (1)  formulating  cold  flow  studies  of  new 
designs,  (2)  determining  particle  sizes  in  engine  effluents  of  alumina  or  other 
condensed  compounds,  (3)  performing  particle  cloud  radiation  studies,  and 
(4)  analyzing  heat  fluxes  in  specific  designs  of  motor  nozzles. 

The  work  focuses  on  the  environment  rather  than  on  the  complicated 
environment -wall  interract ions .  However,  the  techniques  and  new  understandings 
of  this  work  can  now  be  applied  to  studies  where  walls  of  graphite  or  plastics 
are  involved  which  interract  to  modify  the  conditions  at  the  environment -wall 
interface . 


The  radiation  and  particle  impingement  effects  of  alumina  were  studied. 
The  emissivity  of  alumina  particles  was  determined  from  the  melting  point  up  to 
5200  F,  using  a  H2“02“H20-Al202  combustion  system  in  a  2-inch  chamber  without  a 
nozzle.  The  emissivity  data  was  then  arranged  in  a  form  to  allow  calculations 
of  radiative  heat  flux  for  specific  design  cases.  The  radiant  energy  emitted 
from  a  particle  cloud  was  derived  from  a  volume  element  model  developed  by 
R.  A.  Dobbins  (see  Appendix  A) .  General  expressions  representing  the  radiant 
Interchange  between  a  surface  and  absorbing  media  were  developed  using  the  method 
devised  by  J.  T  Bevans  and  R.  V.  Dunkle  (5.4)  using  A.  K.  Oppenhelm's  (5,5) 
network  method. 
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Particle  impingement  on  the  nozzle  wall  was  studied  both  in  cold  flow 
modeling  and  in  motor  firings  (1,25-inch  throat),  Conversion  of  the  hot  firing 
nozzle  surface  thermocouple  measurements  into  heat  flux  data  was  done  using  a 
digital  computer  program  presented  in  Jet  Propulsion  Laboratory  Report  No.  32-257 
by  Powelly  Howell,  and  Irving  (3.1).  It  is  basically  an  analysis  for  radial  modes 
of  transient  conduction.  Thus,  imposed  heat  fluxes  from  actual  impingement  were 
measured  during  firings  and  Indicated  impingements  in  large  motors  were  determined 
from  cold  flow  studies. 

The  cold  flow  studies  used  scale  models  of  nozzles  with  tubular  inserts 
representing  certain  grain  burn-time  instants.  Nitrogen  and  sulfur  hexafluoride 
gases  were  flowed  in  the  models.  Micron- sized  polystyrene  spheres  were  the  model 
particles.  The  analysis  and  assumptions  used  to  specify  modeling  conditions  are 
a  more  thorough  and  sophisticated  development  of  the  work  of  earlier  studies  at 
Philco  in  cold  flow  modeling  and  gas  particle  flow  analysis  (References  2.1,  2,2, 
8,1). 


In  the  case  of  both  radiation  and  particle  impingement,  a  large  part 
of  the  effort  was  devoted  to  technique  development.  The  developments  included 
radiancy  measurement  and  analysis;  special  rocket  engine  and  alumina-water 
slurry  combustion  operations;  micron- size  particles  sampling,  counting,  and 
measuring;  uncooled  molybdenum  calorimetric  nozzle  design;  nozzle  surface 
temperature  measurement ;  nozzle  wall  transient  thermal  analysis  to  determine 
instantaneous  local  heat  fluxes;  cold  flow  modeling  facility  establishment; 
rocket  nozzle  model  design;  model  particles  handling,  injection,  impingement 
pickup,  and  c  anting;  model  wall  static  pressure  contour  measurement;  time 
Integrated  particle  impingement  rate  calculations. 

A  significant  achievement  of  the  work  was  the  development  of  a  particle 
sampler.  This  device  collected  by  high  velocity  impact  all  the  particles  from 
a  very  small  sample  of  the  motor  effluent  directly  from  the  plume,  and  displayed 
the  particles  in  the  complete  size  range  distribution  on  one  microscope  slide. 

ConVvj"'  vVe  heat  transfer  analyses  were  performed  using  a  computer 
program  of  the  t  '  iott ,  Bartz,  Silver  method  (7.4).  These  were  applied  to  the 
hot  firings  for  comparison  with  the  total  heat  fluxes  determined  from  a 
calorimetric  nozzle.  Using  the  cold  flow  modeling  wall  static  pressure  data, 
they  were  applied  t  infer  the  convective  heat  flux  component  for  nozzles  of 
large  motors. 
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Experience  in  the  rocket  industry  with  the  expensive  and  time  consuming 
cut-and-dry  approaches  to  nozzle  design  has  demonstrated  the  need  for  more  ana¬ 
lytical  engineering  design  approaches.  The  lack  of  precise  design  inputs  creates 
both  a  high  testing  cost  and  a  high  burnout  weight  for  reliability  in  large  motors, 
especially  where  man-rated  requirements  exist. 

Test  data  taking  to  contribute  to  analytical  understanding  has  usually 
been  very  limited  for  two  reasons:  (1)  It  is  difficult  to  Justify  instrumenta¬ 
tion  for  analysis,  if  it  imposes  any  Jeopardy  to  the  immediate  hardware  develop¬ 
ment  test  objectives;  and  (2)  The  severe  environment  of  the  motor  often  obviates 
meaningful  measurements.  Nevertheless,  the  rocket  engineering  industry  is 
obligated  to  develop  the  same  sophisticated  analytical  approaches  used  to  great 
effectiveness  in  other  engineering  sciences. 

Subjective  materials  screening  tests  and  subscale  motor  tests  have  little 
practical  meaning,  if  there  is  a  lack  of  scaling  laws  or  aerothermochemical  detail. 
However,  the  level  of  effort  devoted  to  analytical  understanding  has  been  small  in 
the  R  &  D  programs  for  rocket  motors,  and  it  has  been  negligible  in  comparison 
with  the  magnitude  of  the  major  missile  and  booster  programs.  This  Is  evident  from 
the  fact  that  the  level  of  effort  has  never  been  large  enough  to  warrant  a  session 
to  present  such  work  at  a  national  technical  meeting  on  solid  rocket  motors. 

Cold  flow  studies  conducted  by  contractors  in  support  of  their  major 
hardware  programs  have  yielded  interesting  results,  but  never  quantity  tive  data  of 
more  value  than  post  firing  inspections*  The  work  reported  here  should  open  the 
way  to  more  complete  analyses  and  studies  of  current  and  future  designs.  The 
detailed  results  of  such  studies  could  be  the  basis  foe  new  understandings  of  de¬ 
sign  parameters,  which  could  be  usefully  applied  to  advanced  designs. 
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SECTION  2 


COLD  FLOW  MODELING 


I . 1  INTRODUCTION 

There  are  two  main  objectives  in  the  cold  flow  modeling  of  a  rocket 
engine  nozzle.  The  first  is  to  determine  the  alumina  particle  impingement  on 
the  nozzle  wall.  This  points  out  areas  of  potentially  severe  erosion  or  heat 
transfer  and  provides  a  quantitative  figure  for  a  detailed  heat  flux  analysis. 

The  second  objective  is  the  development  of  static  pressure  contours  over  the 
nozzle  for  input  to  connective  heat  transfer  analysis.  In  order  to  accomplish 
this,  a  modeling  analysis  was  derived  to  allow  the  setting  of  valid  modeling 
conditions  using  certain  assumptions ,  and  a  v/r.-j  .'esic,ned  and  built  to 

supply  the  required  modeling  flow  conditions.  Finally,  a  technique  was  developed 
for  handling,  injecting,  picking  up,  and  measuring  the  model  particles. 

Modeling  equations  were  derived  from  an  analysis  of  the  equations  of 
change  for  a  compressible  fluid  and  condensed  particles  in  two-phase  flow.  The 
equations,  written  in  parametric  form,  represent  conditions  within  the  motor 
and  the  model.  For  valid  modeling,  these  parameters  in  the  rocket  and  model 
must  be  equal. 

From  the  modeling  equations,  it  was  observed  that  for  a  reasonable 
model  size,  model  chamber  pressure,  and  particles  larger  than  submicron  sizes, 
the  particles  specific  gravity  must  be  on  the  order  of  one.  This  restriction 
limited  the  choice  of  particles  to  a  low  vapor  pressure  liquid,  atomized  into 
small  droplets,  or  to  solid  polystyrene  spheres  (available  from  Dow  Chemical 
Company)  in  discreet  particle  sizes.  The  polystyrene  particles  were  available 
in  0.09,  0.26,  0.36,  0.56,  0.8,  1.3,  and  2.0  micron  diameters. 
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The  model  chamber  diameter  was  arbitrarily  set  at  eleven  inches.  The 
scale  factor  was  then  set  by  the  chamber  diameter  of  the  hot  rocket  being  modeled. 
It  varied  from  one  to  ten  for  the  casi s  studied.  The  ratio  of  particle  densities, 
model -to-rocket  motor,  was  set  using  the  polystyrene  spheres.  The  working  fluid's 
specific  heat  ratio,  viscosity,  molecular  weight  and  temperature  ratios,  model-to- 
rocket,  were  set  using  unheated  nitrogen  or  sulfur  hexafluoride  gas  in  modeling. 
Geometric  similitude  was  used  in  £^11  configurations.  Thus,  the  remaining  variables 
in  experiment  design  were  chamber  pressure  and  particle  size. 

In  order  to  develop  model  chamber  pressures  below  atmospheric  and  to 
study  exit  cones,  the  facility  included  an  air  ejector  which  could  pump  the  exit 
plane  region  to  1  psia  with  5  lb/ sec  of  nitrogen  flowing  through  the  chamber.  All 
runs  were  fed  by  blowdown  from  a  tank  farm  of  compressed  gas. 

Since  particle  Impingement  on  the  aft-closure  is  dependent  upon  the 
approach  velocity  developed  in  the  grain  port,  it  was  necessary  to  model  a  rocket 
at  several  different  burn  times  and  then  integrate  the  impingement  over  the  total 
firing  time.  At  early  burn  times,  the  approach  velocity  to  the  aft-closure  Is 
higher  than  at  later  burn  times  when  the  grain  port  is  larger  but  mass  flow  is  not 
greatly  different. 

The  main  objective  was  to  determine  the  total  impingement  from  all 
particle  sizes  found  In  a  rocket  engine.  However,  it  was  also  of  interest  to  know 
where  each  size  of  particles  would  Impinge,  and  how  the  Impingement  rate  varied 
from  one  particle  size  to  another.  Throe  polystyrene  particle  sizes,  with  differ¬ 
ent  discrete  sizes  in  each  batch,  were  used  in  different  runs.  The  impingement 
over  the  entire  particle  size  range  could  be  obtained  by  Integrating  over  the 
particle  size  distrlbutl<>n  in  a  solid  fuel  rocket. 

Two  fairly  important  areas  where  the  modeling  was  limited  were  the 
modeling  of  the  boundary  layer  and  the  grain  surface  burning.  Since  it  was  im¬ 
practical  to  have  equal  Reynold  numbers  in  the  model  and  rocket,  the  boundary 
layer  could  not  b?  modeled.  However,  the  Reynold  numbers  were  high  enough  so 
that  the  flow  was  highly  turbulent  in  both  cases,  and  the  stream  line  shapes 
would  be  essentially  the  same. 


The  flow  of  gases  from  a  burning  grain  surface  was  not  exactly  simulated 
in  the  modeling  of  the  grain  port.  However,  uniform  particle  concentration  distri¬ 
bution  across  the  model  grain  port  cross  section  approaching  the  model  nozzle  was 
an  Important  criterion.  The  system  geometry  upstream  of  the  model  grain  port  was 
modified  until  an  even  distribution  of  particles  was  achieved. 

The  polystyrene  particles  were  furnished  In  a  latex  form  with  10  percent 
particles.  By  weighing  a  sample  of  the  latex  for  a  test,  the  exact  number  of 
purtlcles  used  was  known.  The  latex  was  diluted  with  methanol  so  that  the  atomized 
mixture  added  to  the  modeling  gas  would  contain  approximately  one  particle  per  drop. 
The  diluting  liquid  evaporated,  leaving  the  particles  in  the  gas.  This  injection 
of  the  particles  Into  the  gas  was  sequenced  so  that  the  gas  flow  was  already  sta¬ 
bilized  in  the  nozzle  as  the  particles  began  flowing  through. 
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Tapes  that  were  sticky  on  both  sides  were  put  on  the  nosxle  wall  to  act 
as  collection  receptors  for  the  particles  that  Impinged.  After  each  run,  the 
tapes  were  removed  from  the  nozzle  wall,  placed  on  glass  slides,  and  counted  under 
a  microscope.  Using  the  same  two  sided  tape  technique,  bars  were  placed  across 
the  grain  port  to  check  the  particle  concentration  distribution.  Duplicate  runs 
were  made  to  verify  the  reproducibility  of  the  data. 

The  pressure  contours  were  indicated  on  a  manometer  board,  which  was 
photographed  when  the  manometers  had  stabilized.  Since  the  run  time  for  the 
blow  down  system  was  only  about  twenty  seconds,  a  preprassurized  manometer  system 
was  used  to  allow  rapid  response.  Prepressurizing  of  the  manometer  tubes  required 
the  use  of  an  auxiliary  pressure  supply,  three-way  solenoid  valves,  and  bleed 
valves.  Each  manometer  was  prepressurized  by  the  auxiliary  pressure  supply  to 
approximately  the  height  of  liquid  (water  or  mercury)  which  was  anticipated  for 
a  given  pressure  measuring  position  in  the  nozzle.  During  the  run,  the  three-way 
Solenoid  valves  were  switched  to  read  the  nozzle  static  pressures. 

Three  nozzles,  which  will  be  described  in  more  detail  later,  were  modeled. 
A  glmballed  nozzle  was  modeled  in  both  a  straight  position,  and  in  a  5-degree  canted 
position.  A  particle  Impingement  study  on  the  aft-closure  of  Aerojet-General's 
100  FW2  was  also  made. 

A  total  of  157  tests  were  made  on  the  facility.  Table  2.1  is  a  breakdown 
of  the  runs.  Seventy-six  of  the  runs  were  technique  development  runs,  and  produced 
modeling  data. 


2.2  COLD  FLOW  MODELING  ANALYSIS 

The  modeling  equations  were  derived  by  (1)  considering  the  equations  of 
change  through  the  nozzle  for  the  modeling  gas  and  particle  cloud,  (2)  putting  the 
equations  of  change  into  dimensionless  form,  (3)  setting  the  parameters  of  these 
equations  equal  in  the  model  and  rocket,  and  (4)  applying  assumptions  which  allow 
a  set  of  modeling  equations  to  be  derived  from  these  parameters. 

Please  turn  topi£e2-18  fot*  foldout  nomenclature  sheets. 


2-3 


TAiJLE  2.1 


MODELING  TEST 
TEST  OBJECTIVES 


Develop  counting  technique 
Determine  tape  cullecii(n\  cificLency 
Preliminary  check  out  runs 
Develop  pressure  reading  technique 

Check  for  uniform  particle  distribution  across  grain 

Straight  Glmballcd  No^ale 

Pressure  contour  using 

Pressure  contour  using  SF, 

o 

Particle  impingement  using 
Particle  impingement  using  SF^ 

5^  Canted  Gimballed  Nozzle 
Pressure  contour  using 
Particle  impingement  using 
100  FW2 

Particjc  impingement  using 

TOTAL 


Number  of  Tests 
20 
17 
8 
31 
29 

b 

4C. 

23 

2 

10 

3 

6 

157 
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Eqgatlons  of  Chanuc*  Tlirouuh  No/.iclc 

The  equations  of  change  used  were  formulated  first  by 
Marble  (1)  for  general  gas-particle  flow  in  polar  coordinates* 
Marble,  assuming  Stokes  law  to  be  valid,  defined  two  character¬ 
istic  lengths;  the  '’thermal  range”  and  the  ’’momentum  range.” 

These  nkiy  be  viewed  physically  as  tlu?  distances  covered  by  a 
particle  in  reducing  its  relative  temperature  and  velocity  to  e  ^ 
of  its  initial  value.  The  equations  have  subsequently  been  re¬ 
formulated  for  the  axisymmetric  case  and  tlie  cl\avacterislic 
lengths  redefined  to  account  for  deviations  from  Stoke 's  Law. 

These  deviations  were  taken  into  account  by  defining  f  (force 
exerted/part  icle)  =  K,  f..  ,  where  K,.  =5  /C„  ,  and  has 

been  determined  at  Philco  Research  Laboratories  <2)  to  be  given  by 


1  +  0.15  Re 


0.687 


/  L 


1  +  exp 


-0.427 
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4.63 


Re 


0.88 


M 

1  + 

Re 


3.82  +  1,28  exp  (-1.25 


Uith  tills  correction. 


momentum  range  r.  X 

m 


temperature  range  = 


M  a 
-JE _ 
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In  «ddiLlon,  a  correction  for  heat  transfer  effects  was 

introduced . 

Nu  -  2  K„ 

n 

where 


2  +  0.439 

'Si  “  M  1  ^ 

2  +  6.84  2  +  0.459  Re  * 

Re  p 

P  \  ^ 

The  axisynuneCric  equations  of  change  were  then  applied. 
The  continuity  equations,  the  momentum  equations  in  the  radial 
and  axial  directions,  and  tiie  energy  equations  of  the  gas  and 
particle  cloud  are  presented  for  steady*-state  conditions  and 
neglect  wall-gas  and  wall-particle  interactions. 

For  the  gas,  the  continuity  equation  i& 


(pvr)  +  r 


iL 

I'^X 


(/nu) 


0 


the  radial  momenCum  equation  is 


(1) 


f  V  -  V  \ 

;v  ^ 

p 

V 

nr 

r  u 

J 

I  *  —  "N  *  a 

I  tjr  ^  pc 

i  / 

(2) 


Che  axial  momentum  equation  is 
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ar 


+  u 


()u 

dx 


a 

pc 


(3) 
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and  the  energy  equation  is 


For  the  condensed  phase,  which  will  be  created  as  a 
continuum  free  of  particle  interactions,  the  continuity  equation 
is 


1. 

hr 


o  V  r 
/^pc  p  I 


-  0 


radial  momentum  equation  is 


(5) 


(6) 


the  axial  momentum  equation  is 


r3u  riu 


(7) 
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«nd  Che  energy  equation  is 


dr 

— E 
dr 


dX 

+  u 

P  dx. 


•Si 


JL 


*4) 


(8) 


In  addition,  perfect  gas  behavior  is  assumed, 

P  .  -R-  T 
p  M.W. 


(9) 


Dlmena ion less  Form  Developmeat 

The  equations  of  change  may  be  put  in  dimensionless  or 
parametric  form  by  making  the  following  substitutions: 
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ie 


u/U 

★ 

X 

-  x/D 

* 

X 

m 

-  \  /D 
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Up/U 

★ 

T 

P 

-  T  /T 

P  o 

★ 
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P/pU^ 
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V 

-  v/U 
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T 
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Xj./D 

* 

V 

P 

■ 
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The  dimensionless  forms  of  Equations  (1)  through  (8)  are  then 


fi  ^  *  A  ic 

(f,v  r  )  +  r  (pu  )  -  0 


dr 
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(11) 
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^P 
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IC  la  apparent  that  t!ie  solutlona  of  Equations  (10) 
through  (17)  will  ba  Identical  for  two  different  aeta  of  opera¬ 
ting  conditions  only  If  certain  conditions  are  satisfied.  The 
solution  of  Equations  (10)  through  (17)  aiay  be  expressed  by 


Solution  ■  0, 


♦  p£c  a  Ppc  a  U 


^pc  a 


U 


P  U  Cp  •  p  U  Kp 


.  c  T  ’ 
p  o 


a  c 


e  p  ‘Si 

U  *  U  c  Kp  »  «  »  t 


(18) 


But,  since  certain  of  the  slnllarlty  groups  on  the  right  hand 
side  of  Equation  (16)  are  related,  it  Is  possible  to  write 


teluCtOQ  -  I  t".  x*.  ^  ,  §.  ^ 


'S. 


U  ji  j 


(19) 


Equated  Modeling  Paraweters 

Thus,  If  the  parameters  arc  Individually  equal  In  two 
different  systems,  the  solutions  In  terms  of  the  dimensionless 
variables  are  equal.  In  order  to  simulate  rocket  performance  in 
a  model.  It  is  necessary  Chat 
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Four  of  these  similarity  groups  may  be  further  simplified.  Making 
use  of  the  ideal  gas  relationships  and  the  appropriate  defining 
relationships,  it  is  possible  to  demonstrate 
2 

“  '“p  ■'o  •  ^  “ 

a 

and,  since  ^  is  already  constrained,  we  may  state  that  7^  *  7d 
a  n  K 

a  necessary  condition;  a  is  also  a  function  of  7. 

K„/Kj  -  f  Mp) 

and  thus, 

«'p>M  ■  <%>R  •”'*  <"„)m  -  («p>R  • 

R  -y 

P  “  M.W.  <7  -  1)  ’ 
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therefore 


“I'm'I'r  ■ 


C 

Ji 

\  C  /  M 


!  c  1 

I  p  t 

^1  ^  may  be  replaced  by 


(M.W.)j^  -  (M.W.)j^  Cj^ 


*<  Vf  ’ 

Since  \  Is  constrained,  X  ■ 

m  ,  f  i  M 


X^l  may  be  replaced  by 


(Pr)^  -  (Pr)^ 


The  similarity  criteria  may  be  rearranged  and  presented 
in  the  following  form 
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from  (Rep)j^  -  (Re^)^  : 


from  -  (Mp)j^  ; 


'M  ‘M 


X*  M  with  (ROp)^  -  (Rep)j^  and  (Mp)„  -  (Mp)^: 


As*;umpt  ions 


It  is  apparent  from  these  conditions,  which  should  be 
satisfied  at  all  corresponding  points  in  model  and  rocket,  that  to 
achieve  an  exact  modeling  condi:  Ion,  the  modeling  fluid,  modeling 
particulate  phase,  and  model  geometry  must  satisfy  rigid 
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conditions.  To  satisfy  cliese  conditions  exactly  ntay  be  inipos- 
aibla,  but  an  attempt  was  made  to  satisfy  those  which  appear  to 
be  laost  critical  for  effective  modeling. 

Similarity  condition,  Equation  (24),  must  be  aatiafied 
in  order  to  model  fluid  flow  conditions.  If  ■  7j^,  model  and 
rocket  area  ratios,  pressure  ratios  and  temperature  ratios  are 
easily  related  when  simple  geometrical  scaling  is  used.  In  addi> 
tion,  taking 


Pr 


47 

97-5 


the  Prandtl  number  will  be  equal  in  systems  of  equal  7.  Thus, 

satisfying  Equation  (24),  satisfies  Equation  (29),  and  makes 

possible  direct  geometric  scaling.  Thus,  conditions,  Equations 

(20)  and  (21),  are  easily  met. 

The  specific  heat  ratios  encountered  in  hot  rocket 

firings  are  low  (7  -  1.18  for  gas  alone,  7  ^  *1.1  for  gas 

'  ©  «  'effective 

and  condensed  phases)  as  is  the  molecular  weight  of  the  liot  prod* 
ucts.  Unfortunately,  at  the  relatively  low  temperatures  (ambient 
or  slightly  higher)  desired  for  cold-flow  modeling,  only  high 
molecular  weight,  polyatomic  molecules  (e.g.,  CC* ,F2  and  SF^)  have 
specific  heat  ratios  of  the  necessary  value.  These  lugh  molecular 
weights  will  make  It  difficult  to  satisfy  conditions.  Equations 
(22)  and  (27),  since  high-particle  cloud  density  and  extreme  low- 
particle  specific  heat  will  be  required.  High-particle  cloud 
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density  requires  the  use  of  a  Large  number  of  particles  which  is 
an  undesirable  modeling  condition,  since  it  complicates  analysis. 

A  modeling  particle  with  low  specific  heat  (as  compared  to  that  of 
AI2O2  in  rocket  firings)  will  be  difficult  to  find,  since  it 
appears  that  the  necessary  molecular  weight  ratio 
will  be  on  the  order  of  1/6  in  order  to  satisfy  Equation  (24). 

The  effect  of  failing  to  satisfy  the  conditions  imposed  by  Equa¬ 
tions  (22)  and  (27)  appears  to  be  less  severe  than  if  the  equal 
7  condition.  Equation  (24),  were  neglected.  While  7  determines 
the  fluid  flow  field,  particle  cloud  density  and  specific  heat 
serve  to  modify  the  field.  In  addition,  if  it  may  be  assumed  that 
individual  particles  behave  independently  of  surrounding  parti¬ 
cles,  neglect  of  the  particle  cloud  density  as  a  parameter  is 
reasonable . 

The  final  governing  modeling  relationships  can  now  be 
stated  from  Equations  (20)  to  (29)  considering  the  limitations  of 
not  being  able  to  satisfy  Equations  (22)  and  (27). 

From  Equations  (20)  and  (21), 


Scale  Factor 


h 


(30) 


That  is,  there  is  geometric  similarity  in  all  directions. 
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By  equating  Equations  (25)  and  (26),  the  model  pressure 


can  be  expressed  as 


M 


\ 


1/2 


M.W.^ 


'R 


and  Equation  (28)  can  be  arranged  to  give 


M 


j! 


Pp 

R  R 


M 


M  M 


which  is  the  particle  scale  factor. 


(31) 


(32) 
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2,3  FACILITY 


To  accomplish  the  modeling,  a  facility  had  to  be  designed  and  built  which 
would  satisfy  the  previous  mentioned  requirements.  Figure  2-1  showp  a  schematic  of 
the  modeling  facility.  The  particles  were  mixed  with  35  cubic  centimeters  of 
methanol  and  fed  from  two  tanks  during  a  run.  The  indicated  valves  were  opened  for 
the  modeling  gas  with  the  pressure  regulator  controlling  the  flow  to  maintain  the 
desired  model  chamber  pressure.  After  the  flow  was  stabilized,  the  valves  on  the 
bottom  of  the  alcohol-particle  tanks  were  opened  and  the  mixture  wa£  atomized  into 
the  mixing  chamber.  The  particle  atomizing  system  consisted  of  two  sets  of  pneumatic 
atomizing  spray  nozzles,  one  set  of  six  on  each  side  of  the  "T"  shaped  spray  mixing 
chamber. 


As  can  be  seen  in  Figure  2-1,  there  were  baffle  plates  over  the  two 
primary  gas  inlets  located  with  the  atomizing  nozzles.  These  baffles,  shown  in 
greater  detail  in  Figure  2-2,  directed  the  primary  modeling  gas  to  envelop  each 
spray  and  keep  the  spray  from  wetting  the  walls. 

Figure  2-3  is  a  view  of  the  part  of  the  facility  just  described  and  also 
shown  is  the  ejector  gas  valve  which  regulates  the  flow  of  pumping  air  to  the  air 
ejector. 


In  the  spray  mixing  chamber,  the  methanol  was  evaporated  leaving  the  poly¬ 
styrene  spheres  entrained  in  the  modeling  gas.  The  particle  laden  gas  then  flowed 
on  to  the  propellant  grain  port  and  nozzle.  The  test  chamber  and  a  nozzle  model  are 
shown  in  Figure  2-4.  Figure  2-4  also  shows  the  manometer  prepressurizing  system, 
while  Figure  2-5  shows  the  manometer  board  and  camera  used  to  record  the  data.  If 
care  were  taken  in  prepressurizing  the  manometers,  the  time  for  the  manometers  to 
adjust  was  approximately  two  seconds.  A  typical  run  length  was  twenty  seconds. 

Prior  to  a  run,  the  sticky  tapes  for  collecting  impinging  particles  were 
placed  on  the  nozzle  wall  and  the  propellant  grain  port  put  in  place.  A  rubber  seal 
was  used  around  the  tubular  grain  port  to  ensure  a  leakproof  intersection  between 
the  grain  port  and  model .  The  grain  port  was  attached  to  the  model  nozzle  by  a 
system  of  springs  and  the  assembly  was  then  inserted  into  the  test  chanber. 

Figure  2-6  shows  a  grain  port  and  nozzle  being  assembled;  four  particle  collection 
tapes  can  also  be  seen.  The  air  ejector,  used  to  pull  a  vacuum  on  the  exit  cone, 
is  shown  in  Figure  2-7  along  with  the  other  components  of  the  cold  flow  modeling 
facility. 


All  modeling  was  accomplished  by  using  the  blowdown  method* 
hexafluoride,  and  nitrogen  were  used  for  the  modeling  gas.  Table  2.2 
farm  gas  storage  capacities. 


Air,  sulfur 
shows  the  tank 


:.-i7 


NOMENCLATUPJE 


A 

a 

c 


c 

P 


n 

DStokea 


D 

D 


f 

^Stokes 
sub  I 


M 

M 

P 

M.W 


Nu 


? 

Pr 


R 


Re 

Re 

1 

r 


T 

T 


Area 

Local  sonic  velocity 
Particle  specific  heat 
Gas  specific  heat 
Particle  drag  coefficient 

Particle  drag  coefficient 
predicted  by  Stokes  lav 

Characteristic  dimension 

Particle  diameter 

Force  exerted  on  a  particle 

Force  exerted  on  a  particle 
as  predicted  by  Stokes  law 

Heat  of  sublimation 
As  given 
As  given 

Fluid  Mach  number 

Particle  Mach  number  £  u  /a 

8 

Molecular  weight 

Nusselt  number 

Pressure 

Prandtl  number 

Gas  constant 

Fluid  Reynolds  number 

Particle  Reynolds  number  »  D  u  p /ju 

P  s 

Local  radial  dimension 

Local  static  temperature 

Stagnation  temperature  z  T  , 
chamber  temperature 


Tp  Pa rt i c  le  tempera ture 

U  Characterlatic  velocity 

u  Axial  component  of  fluid  velocity 

Axial  component  of  particle  velocity 
Particle  slip  velocity 

V  Radial  component  of  fluid  velocity 

Vp  Radial  component  of  particle  velocity 

W  Mass  flow  rate 

X  Local  axial  dimension 

Greek  Letters 

Specific  heat  ratio 

Momentum  range 

Thermal  range 

Fluid  viscosity 

Fluid  density 

Particle  density 

Density  of  particle  cloud 
(mass  of  solid/gas  volume) 

Subscripts 

c  Condition  in  chamber 

e  Condition  &t  exit  plane 

M  Model  condition 

R  Rocket  condition 

*  Condition  at  nozzle  threat 


7 

X 


P 

^pc 
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FIGURE  2-1.  COLD  FLOW  FACILITY 
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TABLE  2.2 


TANK  FARM  CAPABILITIES 


Gas 

Max.  Pressure 
Ipjt) 

No.  of  25  ft.^ 
Cylinders 

Total  Volume 

at.  3) 

Air 

1500 

21 

525 

Nitrogen 

2200 

18 

450 

300 

3 

75 

2.4  MODEL  PARTICLES 

At  th«  outset  of  the  program,  both  dyed  atomized  liquid  droplets  and 
polystyrene  spheres  were  considered  for  model  particles.  The  liquid  droplets  could 
be  picked  up  on  blotters  placed  along  the  nozzle  model  surface.  Particle  size 
sampling  and  analysis  would  have  been  required.  A  demister  would  have  been  used  to 
remove  over- sized  droplets. 

The  polystyrene  spheres  were  selected  because  their  size  was  predetermined 
no  sampling,  measuring,  or  counting  were  required  to  establish  this.  Sticky  tape 
was  found  to  be  very  effective  in  picking  up  the  impinging  polystyrene  spheres  for 
count ing . 


All  counting  of  the  polystyrene  particles  was  accomplished  with  the  use  of 
a  microscope  at  430X.  An  attempt  was  made  to  dye  the  particles  for  easier  counting. 
However,  it  wes  not  possible  to  get  the  particles  to  absorb  enough  dye. 

Figure  2-8  is  a  photograph  of  the  1.3  micron  particles  at  SOOX;  These 
particles  had  been  impinged  on  tape. 


2.5  TECHNIQUE  DEVELOPMENT  FOR  PARTICLE  IMPRINGEMENT 

Before  meaningful  data  could  be  obtained  from  the  modeling  facility, 
considerable  technique  development  was  required.  The  five  main  areas  of  technique 
development  were  (1)  dispersing  of  polystyrene  particles  In  the  modeling  gas, 

(2) . establ ishlng  the  correct  number  of  particles  In  each  run,  (3)  developing  an 
even  particle  distribution  across  the  grain  port,  (4)  obtaining  a  leak-free 
precision, grain  port-nozzle  Intersection,  and  (S)  counting  the  particles. 

As  mentioned  earlier,  the  polystyrene  particles  were  mixed  in  a  liquid  and 
sprayed  into  the  spray  chamber  where  the  liquid  evaporated  leaving  the  particles  in 
the  gas.  The  two  problems  in  accomplishing  this  were  (1)  spraying  small  enough 
liquid  droplets  to  preclude  settling  out  before  evaporation  and  (2)  obtaining 
unagglomerated  single  part  icier  in  the  gas.  To  get  a  fine  spray  into  the  spray 
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FIGURE  2-8.  1.3  MICRON  POLYSTYRENE  PARTICLES  IMPACTED  ON  TRANSFER 

TAPE,  500  X  MAGNIFICATION.  VIEW  REPRESENTS  0.012  SQUARE 
INCHES  ON  100-INCH  ROCKET  MOTOR  WALL 
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chamber,  twelve  Spraying  Systeras  Company  No.  lA  pneunnatic  atomizing  nozzles  were 
used.  In  the  beginning,  the  particles  were  mixed  with  water,  but  inspection  of  the 
collection  tapes  Indicated  that  not  all  of  the  water  was  evaporating.  This  problem 
was  solved  by  using  methyl  alcohol  instead  of  water  for  the  mixing  fluid  in  order 
to  increase  the  evaporation  rate. 

To  inciicace  whether  or  not  liie  spray  nozzles  were  spraying  excessively 
large  drops*  or  were  wetting  the  chamber  walls*  a  run  was  made  using  dye  in  the 

alcohol.  Blotters  were  attached  to  the  inside  of  the  spiay  chamber.  No  large  dye 

urops  were  found,  indicating  the  spray  was  very  fine,  and  no  appreciable  wall 

splashing  was  indicated.  It  was  interesting  to  \\nte  that  the  dye  gave  a  good  visual 

picture  of  particle  impingement  on  the  nozzle  aft  closure  and  indicated  clearly  the 
point  it  separation  in  the  exit  cone. 

In  order  to  have  no  particle  cjgglomerates  in  tlie  gas,  it  was  first 
necessary  to  assure  that  no  agglomerates  existed  in  the  aicohoi-part icle  mixture. 

It  was  found  that  by  placing  the  alcohol -part icle  mixture  in  an  ultra-sonic  vibrator 
directly  befert  each  run,  the  particles  in  the  mixture  would  be  deagglomerated ,  The 
mixture  of  alcohol  and  latex  was  proportioned  so  that  there  would  be  only  one 
particle  per  drop  of  atomized  liquid  assuming  an  average  drop  size  of  twenty  microns. 
Examination  of  the  tapes  then  indicated  no  appreciable  agglomeration. 

At  first,  the  number  of  particles  in  each  run  was  determined  by  weighing 
the  particle  latex  for  that  one  run  and  then  mixing  it  with  the  alcohol.  It  became 
apparent  that  small  errors  n  the  weight  of  the  sample  could  cause  rather  large 
errors  in  the  *mmber  of  particles.  Therefore,  it  was  necessary  to  mix  a  large 
quantity  of  particle  latex  with  the  alcohol  and  then  divide  the  alcohol- 
particle  mixture  into  smaller  samples.  In  this  manner,  small  errors  in  measurement 
T^ide  negligible  errors  in  the  number  of  particles. 


Obtaining  an  even  part ic le  distribution  in  the  grain  port  required  that 
the  particles  be  thoroughly  mixed  in  the  modeling  gas  and  then  having  a  flow 
path  which  would  not  cause  particle-gas  separation  ahead  of  the  grain  port  model. 


At  first,  the  particle  concentration  was  very  high  near  the  grain  port 
wall.  Three  changes  were  then  made  in  the  test  facility  to  correct  this.  The  first 
change  was  to  put  baf'^e  plates  (Figure  2-2)  over  the  two  gas  inlets  near  the  pneu¬ 
matic  spray  nozzles  lie  baffle  plate  divided  the  main  gas  inlet  stream  into  six 

each  spray  nuzzle.  The  second  Improvement  was  to  turn  the 
Lhe  rear  as  shown  in  Figure  2-1.  When  the  inlet  was  toward 
ntly  produced  a  core  of  part ic le-f ree  gas.  The  third  im- 
.  f fled  5  tubular  insert  into  che  entrance  to  the  modeling 
chamber,  also  shown  in  Figure  2-1.  This  device  was  designed  to  break  up  any 
boundary  layer  particle  concentrations  developed  in  the  6-inch  pipe  approaching  the 
test  section.  Tn^^se  three  improvements  resulted  in  a  very  even  particle  concentra¬ 
tion  across  the  grain  port  flow. 


streams,  one  sui'^ 
center  gas  inlet 
the  test  nozzle,  f 
provement  was  to  put 
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Two  grain  porL  wall  configurations  were  tried  before  the  above  improve¬ 
ments  were  made.  One  was  a  solid  wall  tube  and  the  other  was  a  perforated  wall 
vjhich  was  intended  to  simulate  more  closely  the  flow  from  a  burning  grain.  After 
an  even  concentration  distribution  was  found  in  the  solid  wall  grain  port,  it  was 
chosen  for  testing.  Furthermore,  it  was  much  easier  to  check  the  distribution 
across  this  type  of  grain. 

For  the  first  tests  little  attention  was  paid  to  the  grain  port  nozzle 
intersection.  It  soon  became  apparent  that  this  intersection  had  to  be  a  leak- 
free  precision  fit  to  avoid  secondary  flow  disturbances  which  would  upset  particle 
impingement,  A  leak-free  seal  was  obtained  by  placing  the  grain  in  the  nozzle  and 
then  pouring  in  a  polymerizing  liquid  rubber  behind  the  grain  port. 

Various  attempts  were  made  to  ease  the  job  of  counting  the  particles.  Ini¬ 
tially,  the  particles  were  photographed  by  use  of  a  microscope,  and  the  particles 
counted  on  the  photograph.  However,  this  left  too  much  room  for  error  since  a  par¬ 
ticle  slightly  out  of  focus  would  look  like  a  blotch  on  the  tape.  Also,  a  speck  of 
dirt  out  of  focus  would  look  like  a  particle.  For  accurate  counting,  the  plane  of 
focus  had  to  be  changed  while  counting  which  meant  the  particles  had  to  be  counted  by 
eye  directly  with  a  microscope.  The  microscope  used  had  a  43X  objective  with  lOX  eye 
pieces,  A  grid  in  one  eye  piece  was  used  to  give  a  constant  known  counting  area, 

2.6  PARTICLE  IMPINGEMENT  CALIBRATION  TESTS 

To  monitor  particle  flow  through  the  grain  port  for  uniformity,  impinge¬ 
ment  samples  were  obtained  along  l/^“inch  diameter  bars  placed  perpendicular  to  the 
flow,  Figure  2-9.  This  data  also  lent,  itself  to  a  comparison  of  particle  size  and 
velocity  effects  on  impingement. 

Figure  2-10  shows  data  obtained  from  the  vertical  and  short  horizontal  bars 
in  the  largest  diameter  grain  port  used  (7,8  inches),  which  gives  the  most  severe 
test  for  a  uniform  particle  concentration  across  the  port. 


Table  2.3  shows  conditions  for  five  of  the  runs  which  used  the  calibra¬ 
tion  bars.  The  actual  data  taken  from  these  bars  are  shown  in  Figures  2-11  through 
2-15.  For  all  five  runs,  the  total  number  of  particles  in  each  run  was  10  billion. 
Therefore,  when  Runs  041575  and  041576  are  coirpared  to  the  other  three  runs,  the 
difference  in  particle  flux  due  to  different  grain  ports  must  be  considered. 

By  comparing  Run  040971  to  041678  and  Run  041575  to  041576,  the  effect  of 
pal  tide  size  can  be  seen  on  impingement  at  the  stagnation  point  on  the  bar  at  the 
same  velocity.  This  shows  that  impingement  is  proportional  to  the  third  power  of 
the  particle  diameter.  The  velocity  effects  on  impingement  can  be  seen  by  comparing 
runs  at  different  velocity  with  the  particle  size  constai  By  comparing  Run  041374 
to  041678  and  Run  041574  to  041575,  it  will  be  seen  that  /ingement  is  approximately 
proportional  to  the  first  power  of  the  velocity. 
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FIGURE  2-9.  GRAIN  FORT  MODEL  (LOOKING  UPSTREAM  FROM  NOZZLE) 


FIGURE  2-10.  TYPICAL  PARTICLE  CONCENTRATION  DISTRIBUTION  IN  GRAIN 
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FIGURE  2-11.  IMPINGEMENT  AROUNb  CALIBRATION  BAR  FOR  RUN  040971 
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2-12.  mPINGEMEMT  AROUND  CALIBRATION  BAR  FOR  RON  041574 


FIGURE  2-14.  IMPIRGEMENT  AROUND  CALIBRATICW  FOR  RUN  041576 


TABLE  2.3 


CONDITION  OF  RUNS  WITH  CALIBRATION  BARS 


Run  No. 

Polystyrene 
Particle  Size 

Grain  Port  Size 
Dia .  ,  in . 

Grain  Port 
Velocity,  ft/sec. 

Impingement  at 
Stagnation  Point 

040971 

2.05  M 

5.2 

165 

160 

041574 

1.3  M 

5.2 

258 

75 

041575 

1 .3  |U 

7.8 

113 

18 

041576 

2.05  /i 

7.8 

113 

75 

041678 

1.3 

5  2 

165 

40 

2.7  PARTICLE  IMPINGEMENT  TESTS 

Particle  impingement  studies  vere  made  on  iour  nozzles:  (1)  a  model  oi 
the  nozzle  used  in  hot  tiring  studies  of  particle  impingement,  (2)  a  gimballed 
nozzle  (straight),  (3)  a  gimballed  nozzle  (canted  5  degrees)  and  (4)  a  model  of  the 
aft  closure  of  the  Aerojet-General  100  FW2  engine. 

As  mentioned  earlier,  both  nitrogen  and  sulfur  hexafluoride  (SFg)  were 

used  as  the  modeling  gas.  A  summary  of  N_  and  SF,  properties  is  shown  in  Table  2.4. 

z  o 


TABLE  2.4 

SUMMARY  OF  N^  AND  SF,  PROKERTIES 


Gs  s 


Temperature  (T  )^R 
7 

Molecular  Weight 
Viscosity  Ih/ft  sec 


480 

1.4 

28 

1.2  X  10 


500 

1.1 

146 

1.01  X  lO’^ 
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2,7,1  S IMULATOR  NOZZLE 


The  simulator  nozzle  (Figure  2-16)  was  a  full  scale  model  of  the  nozzle 
used  in  the  hot  firing  part  of  this  program.  All  impingement  found  on  this  nozzle 
was  negligible  when  3  fj,  particles  were  modeled.  This  indicated  that  large 
agglomerated  particles  must  have  been  produced  in  the  hot  firing  to  produce  the 
Impingement  indicated  in  those  tests. 

2.7.2  GIMBALLED  NOZZLE 

In  the  gimballed  nozzle  (Figure  2-17),  there  was  no  impingement  in  the 
gimballing  section  of  the  nozzle  in  either  the  straight  or  canted  5  degree  position. 
Therefore,  a  full  impingement  study  was  made  only  on  the  entrance  section  of  the 
straight  nozzle. 

It  was  questioned  whether  or  not  the  absence  of  impingement  in  the  gim¬ 
balling  section  was  due  to  no  particles  impinging  on  the  tapes  or  the  velocity  was 
too  high  and  the  particles  bounced  off  the  tape.  A  90  degree  wedge  and  a  bar  were 
placed  in  the  exit  cone  (Figure  2-18)  of  the  gimballed  nozzle  to  check  whether  or 
not  any  particles  could  be  picked  up  in  the  supersonic  region.  Two  runs  were  made, 
one  with  sulfur  hexafluoride  and  one  with  nitrogen  as  the  modeling  gas.  A  high 
concentration  of  particles  were  found  on  the  wedge  and  bar  for  the  sulfur  hexafluo¬ 
ride  run,  but  none  was  found  for  the  nitrogen  run.  Since  the  molecular  weight  of 
sulfur  hexafluoride  is  larger  than  nitrogen,  the  sonic  velocity  would  have  been 
lower  and  thus  the  velocity  at  the  wedge  smaller.  This  indicates  that  the  higher 
velocity  particles  in  the  nitrogen  flow  may  not  have  been  picked  up.  On  the  other 
hand,  impingement  was  also  negligible  on  the  gimballing  part  when  sulfur  hexafluo¬ 
ride  was  used.  (One  explanation  of  the  difference  in  impingement  on  the  shapes  in 
the  exit  cone  with  the  two  gases  is  that  the  differences  are  a  function  of  the 
differences  In  shock  lines  produced  and  their  effect  on  the  particles.) 

The  gimballed  nozzle  which  was  modeled  had  a  63-inch  diameter  chamber 
with  a  15-inch  gimballed  threat.  It  was  assumed  that  the  rocket  had  a  total  burn 
time  of  90  seconds,  a  burn  rate  of  500  lb/ sec.,  and  an  alumina  loading  of  32  percent. 
Since  for  subsonic  flow  the  impingement  is  nearly  the  same  when  nitrogen  and  sulfur 
hexafluoride  is  the  modeling  gas,  nitrogen  was  used  for  the  study. 

From  the  modeling  equations  (Equations  2.31  and  2.32),  the  particle  dia¬ 
meter  ratio  and  model  pressure  can  be  deter.iined  for  the  case  of  nitrogen  as  the 
modeling  gas  and  Dow  polystyrf'ne  spheres  as  the  particles.  Since  the  model  chamber 
is  11  inches  in  diameter,  a  scale  factor  of  =  5.73  is  dictated.  The  other 

valies  needed  to  s^lve  Equation8(2.3l)and(2.3^are: 

-  2.28  gm/^m^  T  -  6360°R 

p„  -  1.05  gm/cm^  T„  »  480°R 

Pm  “ 
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4100  fps 


MWj^  -  20.4 


a  -  992  fps 

■  6.7  X  10  ^  Ib/ft  sec. 


“  1.06  X  10  Ib/ft  sec. 
Pj^  “  613  psla 


HW^  -  28 
n 


7«  -  1.J8 


Therefore,  from  the  modeling  equations,  the  particle  diameter  ratio  and 
the  model  chamber  pressure  are 


_1 _  ^  2.28  4100 

5.73  1.05  992 


1.06  X  10 

6.7  X  10"' 


1/2 


0.496 


(2.33) 


and 


1 

0.496  * 


1.06  X  10 

6.7  X  10 


X 


480 

6360 


I 


-,1/2 


20.4  1.18 

28  1.4 


X  613  psla 


(2.34) 
42  psla 


In  Equation  (2.32),  p  and  D  correspond  to  the  conditions  in  the 

rocket  at  6360  R>and  when  a  certain  particle  diameter  Is  being  modeled,  it  refers  to 
the  diameter  of  the  particle  at  room  temperature  (530°R) .  For  example,  if  a  Zp 
diameter  particle  is  being  modeled,  the  diameter  of  the  particle  in  the  model  must 
be  determined  with  the  knowledge  that  the  hot  particle  in  the  rocket  Is  actually 
larger  than  2p.  The  hot  diameter  can  be  determined  by  conserving  mass. 


or 


Mass  of  hot  particles  ■  mass  of  cold  particle 


cold 


D_^|  (2.28).  -  (2)^  (3.7) 

}  hot 


jD  I  -  2.35M 

'  hot 


(2.35) 


2-3^ 


FIGURE  2-16.  SIMULATOR  MODEL 


FIGURE  2-18.  PARTICLE  IMPINGEMENT  TAPES  ON  WEDGE  AND  BAR  FOR  SUPERSONIC 
REGION  OF  NOZZLE 
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Therefore,  to  model  a  cold  2/^  dlamacar  particle  with  a  denaity  of 
3.7  gm/cm  ,  p  and  D  must  be  2.26  gm/cm^  and  2.3SP,  respectively. 

R  R 

The  density  for  aluminum  oxide  at  6360°R  was  taken  by  extrapolating  the 
data  of  Klrshenbaum  and  Cahill. 

For  the  three  sizes  of  polystyrene  particles  used,  Table  2.5  shows  the 
sizes  of  aluminum  oxide  particles  being  modeled. 


TABLE  2.5 


MODELED  PARTICLE  SIZES  FOR  GIMBALLED  NOZZLE 


Polystyrene 
Particle  Dia. 

Hot  Aluminum  Oxide  Particle 
Dla .  in  Rocket 

Cold  Aluminum  Oxide 
Particle  Dla.  from  Rocket 

0.8 

1  .DM 

1  .4m 

1.3  M 

2.6m 

2.2m 

2.05M 

4.1M 

3.5m 

As  shown  In  Figure  2'19,  three  grain  port  sizes  representing  three 
different  instants  In  burn-time  (0,  23,  and  58  percent  of  burn-out)  were  needed  to 
show  the  variation  of  Impingement  with  time.  Since  three  grains  and  three  particle 
sizes  were  used,  nine  runs  were  made  to  develop  the  overall  impingement  picture. 


Figures  2-20,  2-21,  and  2-22,  present  the  impingesMnt  data  of  the  1.4M  » 
2.2p,  and  3.5p  particles,  respectively.  The  ordinate  is  expressed  as  the  weight 
percent  of  that  particular  particle  size,  at  that  particular  Instant  of  tloM, 
Impinging  on  the  aft  closure.  Each  data  point  represents  the  average  of  from  6  to 
12  counts  of  impinged  particle  density  in  the  region  of  the  particular  position 
along  the  aft  closure. 


An  exsmple  is  given  to  illustrate  how  the  particles  which  are  impinged  on 
the  tape  of  the  model  can  be  converted  to  weight  percent  impinged  on  the  aft  closure 
of  the  rocket.  Assusw  counting  by  use  of  a  microscope  gives  N  particles/ sq.  in.  on 
the  tape  at  some  point  along  the  aft  closure.  Also,  it  is  known  that  10^^  particles 
were  originally  injected  Into  the  gas  and  the  scale  factor  from  the  model  to  the 
rocket  is  5.75.  Then  the  weight  percent  Impinged  is 


1 

(5.75)^ 


x 


100. 


It  is  also  interesting  to  note  how  reproducible  the  data  are  by  comparing 
one  person's  particle  count  with  another's  and  by  comparing  two  tapes  from  Che  same 
run  or  from  duplicate  runs.  These  comiarisons,  shown  In  Figure  2-23  indicate 
reproducibility  of  the  data. 
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OX  BURNOIT 


DISTANCE  FROM  (IN.)  R07249 


FIGURE  2-20.  WEIGHT  PERCENT  OF  PARTICLES  IMPINGING 

ALONG  THE  AFT  CLOSURE 
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WEIGHT  PERCElfT  IMPINGING  (7,/IN^)  x 


DISTANCE  FKOH  (INCHES)  R03398 

FIGURE  2-21.  WEIGHT  PERCENT  OF  2.2/i  PARTICLES  IMPINGING 
ALONG  THE  API'  CLOSURE 
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C  B  A  Tj  T2 

DISTANCE  FROM  T^,  (INCHES) 

R05597 

FIGURE  2-22.  WEIGHT  PERCENT  OF  3.5^  PARTICLES 
IMPINGING  ALONG  THE  AFT  CLOSURE 
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WSIGHT  PERCENT  IMPINGING 
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From  Figures  2-20,  2-21,  and  2-22,  ic  is  now  possible  to  develop  a  curve 
of  impingement  versus  grain  burn-time  for  each  part'cle  size  at  any  point  along  the 
aft  closure.  Figures  2-24,  2-25,  and  2-26  present  such  curves  for  points  A,  B,  and 
C,  of  Figure  2-19,  respectively.  The  ordinates  of  these  curves  have  been  changed  to 
weight  fraction  impinging  per  square  inch  Instead  of  weight  percent  impinging  per 
square  inch  as  in  Figures  2-20,  2-21,  and  2-22. 


From  the  curves  in  Figures  2-24  through  2-26,  an  average  weight  fraction 
impinged  can  be  found  from  each  particle  size  at  poinus  A,  B,  and  C  by  noting: 


100-tb 


6w.  d") 

1  ^ 

5  w 


100 
6  w. 


(2-36) 


where  D  is  a  particular  particle  size 


(°} 
6w  (d) 


average  weight  fraction  impinged  per  sq.  in.  for 
this  particle  size 


6w.  (ji) 

6«  {d> 


weight  fraction  impinged  per  sq.  in.  for  the  same 
particle  size  as  a  function  of  time 


t 


b 


time,  percent  of  burnout  when  the  point  on  the  aft 
closure  is  first  exposed  to  particle  impingement 


The  integral 


100 


Figures  2-24,  2-25,  and  2-26.  Table 


^t^  dt  is  the 


area  under  the  curves  of 


6wi  (d) 

2.6  is  a  tabulation  of  — ■ — ■r-—  . 

6w  {d} 


6Wj  (d) 

From  Table  2.5,  it  is  possible  to  obtain  - ®  function  of  par- 

6W 


tide  size  for  points  A,  B,  and  C.  These  curves  are  shown  in  Figure  2-27.  As  will 
be  explained  later,  the  particle  size  range  of  Interest  is  from  0.5/1  to  5.5/1  . 


2-50 


[ 

c 

c 

[ 

t 

[ 

I 

t 

I 

I 

I 


9OI 


{ai“g 

1^1  |at^9 


9NI0NI«(H1  N0U3VU  IHOIVl 


fk. 


2-51 


2-24.  WEIGHT  mCTION  OF  1.4^  .  2.2^  AHD  3.5^,  PARTICLES  mPIHGSD  AT  POIBT  A  VS.  PnCEMT  OF 


mtA 


I 


PARTICLE  SIZE  (D)  M  R05628 


FIGURE  2-27.  AT  POINTS  A,  B,  C:  ALUMINA  WEIGHT  FRACTION 
IMPINGED  AS  A  FUNCTION  OF  PARTICLE  SIZE 


TABLE  2.6 


'W  (D^ 
TABL«LATION  OF  — ^  '  ;; 

6W  {d> 


Point 

*^b  (Z  of 
Burnout ) 

D 

(u  ) 

100  ^ 

/  -  (t)  dt 

so} 

(Z  of  Burnout) 
in^ 

(d) 

6w  (d} 

A 

0 

1.4 

254 

2.54 

A 

0 

2.2 

285 

2.85 

A 

0 

3.5 

168 

1.68 

B 

23 

1.4 

49.6 

0.65 

B 

23 

2.2 

59 

0.77 

B 

23 

3.5 

115.4 

1.5 

C 

58 

1.4 

2.1 

0.05 

C 

58 

2.2 

4.6 

0.115 

C 

58 

3.5 

16 

0.4 

The  next  step  is  to  Integrate  impingement  over  the  size  range  of  particles 
found  in  a  rocket.  However,  to  do  this  the  weight  fraction  of  particles  as  a  function 
of  particle  size  be  known.  This  function  can  be  obtained  by  differentiating  a 

curve  of  Sehgal's^*^  data  for  a  500  psi  chamber  pressure.  The  differential  of  his 
percent  total  volume  below  curve  is  shown  in  Figure  2-28  with  units  of  the  ordinate 
as  I/M  since 

6W  <D>  \  T  ! 

"  dD 

where  is  the  total  weight  of  all  the  alumina  particles. 

By  multiplying  the  ordinate  of  the  curves  in  Figure  2-27  with  the  ordinate 
of  the  curve  in  Figure  2-28  at  each  particle  size  and  plotting  as  a  function  of 
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WEIGHT  mcnon  of  particles 


PARTICLE  SIZE.  D,  Qi) 


R0S624 


FIGURE  2-28.  PARTICU  SIZE  DISTRIBUTION  OF  ALIMINA  SAMPLE  OBTAINED 
FROM  MOTOR  FIRING  AT  SOO-FSI  CHAMBER  PRESSURE 
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particle  size  while  still  keeping  points  A,  B,  and  C  a^  paraoieters,  Figure  2-29  la 
obtained.  The  ordinate  of  Figure  2*29  will  be  ^  1;  ^  •-  ■  A  ■  with 

,  w  D  "x  "t 

units  of  l/jX  in**. 

For  any  one  of  the  three  points  (A,  B,  C)  It  is  now  possible  to  integrate 
iupingement  over  the  entire  particle  size  range.  The  area  under  these  curves  can  be 
expreseed  as 


•  5: 

T 


which  is  the  weight  fraction  of  alumina  that  impinged  at  any  point.  For  a  run  tine 
of  90  seconds^  a  burn  rate  of  500  Ib/sec..  and  32  percent  of  the  propellant  aa 
alumina,  the  total  impingement  in  Ib/aq.  in.  can  be  calculated  as  a  function  of 
position.  For  example: 


Impingement  "  ~ 


500  lb. \  (0..'^2)  (90  sec.) 


-  ~  (14400)  Ib/aq.  in. 

^T 


This  curve  for  impingement  along  the  nozzle  wall  integrated  for  the  total  firing  is 
shown  in  Figure  2-30.  This  is  the  final  result  desired  from  the  cold-flow  modeling 
test  data. 


2.7.3  100  EW2 

The  model  of  the  100  FW2  consisted  only  of  the  aft  closure.  Therefore, 
nitrogen  was  used  as  the  modeling  gas.  As  shown  in  Figure  2-31,  the  aft  closure  was 
modeled  at  two  burn  times,  16  and  54  percent  of  burnout.  Three  perticle  sizee  were 
modeled,  making  a  total  of  six  runs. 

Again,  from  the  modeling  equations  (2.31)  and  (2.32)  the  particle  dlanieter 
ratio  and  model  pressure  can  be  determined.  Since  the  model  chamber  la  11  inches  in 

Dr 

diameter,  a  scale  factor  of  ~  •  9.1  is  dictated.  The  other  values  which  are  input 

''m 

to  the  equations  are  as  follows: 

-  2.28  gm/cm^  T,,  -  6360°R 

Fr  R 
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2,7.1  SIMULATOR  NOZZLE 


The  simulator  nozzle  (Figure  2-16)  was  a  full  scale  model  of  the  nozzle 
used  in  the  hot  firing  part  of  this  program.  All  impingement  found  on  this  nozzle 
was  negligible  when  3  jj,  particles  were  modeled.  This  indicated  that  large 
agglomerated  part  ides  must  have  been  produced  in  the  hoi  firing  to  produce  the 
impingement  indicated  in  those  tests. 


2.7.2  GIMBALLED  NOZZLE 

In  the  gimballed  nozzle  (Figure  2-17),  there  was  no  impingement  in  the 
gimballing  section  of  the  nozzle  in  either  the  straight  or  canted  5  degree  position. 
Therefore,  a  full  impingement  study  was  made  only  on  the  entrance  section  of  the 
straight  nozzle. 

It  was  questioned  whether  or  not  the  absence  of  impingemeni  in  the  gim¬ 
balling  section  was  due  to  no  particles  impinging  on  the  tapes  or  the  velocity  was 
too  high  and  the  particles  bounced  off  the  tape,  A  90  degree  wedge  and  a  bar  were 
placed  in  the  exit  cone  (Figure  2-18)  of  the  gimballed  nozzle  to  check  whether  or 
not  any  particles  could  be  picked  up  in  the  supersonic  region.  Two  runs  were  made, 
one  with  sulfur  hexafluoride  and  one  with  nitrogen  as  the  modeling  gas,  A  high 
concentration  of  particles  were  found  on  the  wedge  and  bar  for  the  sulfur  hexafluo¬ 
ride  run,  but  none  was  found  for  the  nitrogen  run.  Since  the  molecular  weight  of 
sulfur  hexafluoride  is  larger  than  nitrogen,  the  sonic  velocity  would  have  been 
lower  and  thus  the  velocity  at  the  wedge  smaller.  This  indicates  that  the  higher 
velocity  particles  in  the  nitrogen  flow  may  not  have  been  picked  up.  On  the  other 
hand,  impingement  was  also  negligible  on  the  gimballing  part  when  sulfur  hexafluo¬ 
ride  was  used.  (One  explanation  of  the  difference  in  impingement  on  the  shapes  in 
the  exit  cone  with  the  two  gases  is  that  the  differences  are  a  function  of  the 
differences  in  shock  lines  produced  and  their  effect  on  the  particles.) 

The  gimballed  nozzle  which  was  modeled  had  a  63-inch  diameter  chamber 
with  a  13-inch  gimballed  threat.  It  was  assumed  that  the  rocket  had  a  total  burn 
time  of  90  seconds,  a  burn  rate  of  500  lb/ s^c . ,  and  an  alumina  loading  of  32  percent. 
Since  for  subsonic  flow  the  impingement  is  nearly  the  same  when  nitrogen  and  sulfur 
hexafluoride  is  the  modeling  gas,  nitrogen  was  used  for  the  study 

From  the  modeling  equations  (Equations  2.31  and  2  32),  the  particle  dia¬ 
meter  ratio  and  model  pressure  can  be  deter.iined  for  the  case  of  nitrogen  as  the 
modeling  gas  and  Dow  polystyrene  spheres  as  the  particles.  Since  the  model  chamber 
is  11  inches  in  diemeter  ,  a  scale  factor  of  =  5.73  is  dictated.  The  other 

val  les  needed  to  s^^lve  EquationsC  2 . 31)  and(  2 . 32)  are  : 

p„  -  2.28  gm/.m^  T  -  6360°R 

p_  ■  1.05  gm/cm^  T„  "  480°R 
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«  6.7  X  lo'^  Ib/ft  sec. 
Uj^  =  1.06  X  10  ^  lb/ ft  sec. 
Pj^  *  613  psla 


MWj^  -  20.4 


2S 

1.18 
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Therefore,  from  the  modeling  equations,  the  particle  diameter  ratio  and 
the  model  chamber  pressure  are 


D. 


M 


_1 2.28  4100  1.06  X  10 

5.73  ^  1.05  *  992  *  ,  ,  .^-5 


6.7  X  10 


1/2 


-  0.496 


(2.33) 


and 


M  0.496 


1.06  X  10 
X  - - - T—  X 


6.7  X  10 


"5 


480 

6360 


20.4  1.18 


28 


1.4 


1/2 


(2.34) 


X  613  psia  42  psia 


In  Equation  (2.32),  p  and  D  correspond  to  the  conditions  in  the 

rocket  at  6360  R,and  when  a  certain  particle  diameter  is  being  modeled,  it  refers  to 
the  diameter  of  the  particle  at  room  temperature  (530*^R)  .  For  example,  if  a  2/j. 
diameter  particle  is  being  modeled,  the  diameter  of  the  particle  in  the  model  must 
be  determined  with  the  knowledge  that  the  hot  particle  in  the  rocket  is  actually 
larger  than  2fi.  The  hot  diameter  can  be  determined  by  conserving  mass. 


or 


Mass  of  hot  particles  ■  mass  of  cold  particle 


(2.28).  -  (2)^  (3.7) 

JD  )  -  2.35P 

1  hot 

'  hot 

(2.35) 


O 
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FIGURE  2-16.  SIMULATOR  MODEL 
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FIGURE  2-18.  PARTICLE  IMPINGEMENT  TAPES  ON  WEDGE  AND  BAR  FOR  SUPERSONIC 
REGION  OF  NOZZLE 
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Therefore,  to  model  a  cold  2M  diameter  particle  with  a  denalty  of 
3.7  gm/cm^,  p_  and  D.  muat  be  2.28  gm/cwr  and  2.33P,  respectively. 


ie  density  £ui  aluniouir 


data  of  Klrshenbaum  and  Cahill 


o.xldc  at 


6360^R  was 


taken  by  excrapoiatlng  the 


For  the  three  slees  of  polystyrene  particles  used, 
sizes  of  aluminum  oxide  particles  being  modeled. 


Table  2.5  shows  the 


TABLE  2.5 


MODELED  PARTICLE  SIZES  FOR  GIMBALLED  NOZZLE 


Polystyrene 
Particle  Dla. 

Hot  Aluminum  Oxide  Particle 
Dia .  in  Rocket 

Cold  Aluminum  Oxide 
Particle  Dia.  from  Rocket 

0.8  n 

1  .6m 

1.4m 

1.3  ju 

2.6m 

2.2m 

2.05m 

4.1M 

3.5M 

As  shown  in  Figure  2' 19,  three  grain  port  sizes  representing  three 
different  instants  in  burn-time  (0,  23,  and  58  percent  of  burn-out)  were  needed  to 
show  the  variation  of  impingement  with  time.  Since  three  grains  and  three  particle 
sizes  were  used,  nine  runs  were  made  to  develop  the  overall  impingement  picture. 

Figures  2-20,  2-21,  and  2-22,  present  the  impingement  data  of  the  1.4/t(  , 
2.2ix,  and  3.Sfx  particles,  respectively.  The  ordinate  is  expressed  as  the  weight 
percent  of  that  particular  particle  size,  at  that  particular  instant  of  time, 
impinging  on  the  aft  closure.  Each  data  point  represents  the  average  of  from  6  to 
12  ccunts  of  impinged  particle  density  in  the  region  of  the  particular  position 
along  the  aft  closure. 


An  example  is  given  to  illustrate  how  the  particles  which  are  impinged  on 
the  tape  of  the  model  can  be  converted  to  weight  percent  impinged  on  the  aft  closure 
of  the  rocket.  Assume  counting  by  use  of  a  microscope  gives  N  particles/sq.  in.  on 
the  tape  at  some  point  along  the  aft  closure.  Also,  it  is  knotm  that  10*^  particles 
were  originally  injected  into  the  gas  and  the  scale  factor  from  the  aodel  to  the 
rocket  is  5.75.  Then  the  weight  percent  impinged  is 


_N 

10 


10 


(5./5)‘ 


100. 


It  is  also  interesting  to  note  how  reproducible  the  date  are  by  comparing 
one  person's  particle  count  with  another's  and  by  comparing  two  tapes  from  the  same 
run  or  from  duplicate  runs.  These  comiarisons,  shown  in  Figure  2-23  indicate 
reproducibility  of  the  data. 
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FIGURE  2-20.  WEIGHT  PERCENT  OF  1,4..  PARTICLES  IMPINGING 
ALONG  THE  AFT  CLOSURE 
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FIGURE  2-21.  WEIGHT  PERCENT  OF  2.2/1  PARTICLES  IMPINGING 
ALONG  THE  AFl*  CLOSURE 
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FIGURE  2-22.  WEIGHT  PERCENT  OF  3.5/i  PARTICLES 
IMPINGING  ALONG  THE  AFT  CLOSURE 
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FIGURE  2-23.  REPRODUCIBILITY  OF  IMPIM6BMBMT  DATA 


From  Figures  2-20,  2-21,  and  2-22,  Ic  Is  now  possible  to  develop  a  curve 
of  impingement  versus  grain  burn-time  for  each  part'cle  size  at  any  point  along  the 
aft  closure.  Figures  2-24,  2-25,  and  2-26  present  such  curves  for  points  A,  B,  and 
C,  of  Figure  2-19,  respectively.  The  ordinates  of  these  curves  have  been  changed  to 
weight  fraction  impinging  per  square  inch  instead  of  weight  percent  Impinging  per 
square  inch  as  in  Figures  2-20,  2-21,  and  2-22. 

From  the  curves  in  Figures  2-24  through  2-26,  an  average  weight  fraction 
impinged  can  be  found  from  each  particle  size  at  poin;:s  A,  B,  and  C  by  noting: 
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1  ^  ^ 

6  w 


(2.36) 


where  D  is  a  particular  particle  size 

6Wi  (d} 

.  ••  average  weight  fraction  Impinged  per  sq.  in.  for 

6w  this  particle  size 

6wi 

'  -■  ^t^  •  weight  fraction  impinged  per  sq .  in.  for  the  same 

6w  particle  size  as  a  function  of  time 


t^  ■  time,  percent  of  burnout  when  the  point  on  the  aft 

closure  is  first  exposed  to  particle  impingement 


The  integral 


f  6.^(d} 

(  (d) 


^t^  dt  is  the  area  under  the  curves  of 


Figures  2-24,  2-25,  and  2-26.  Table  2.6  is  a  tabulation  of 


6w 


«»i  (“) 

From  Table  2.5,  it  is  possible  to  obtain  - 7^”  ®  function  of  par- 

6W 


tide  size  for  points  A,  B,  and  C.  These  curves  are  shown  in  Figure  2-27.  As  will 
be  explained  later,  the  particle  size  range  of  interest  is  from  O.Sjj.  to  5.5/.1  . 
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FIG»E  2-26.  WEIGHT  FRACTION  OF  1.4/X,  2.2/1,  AND  3 .4/1  PARTICLES  IMPINGING  AT  POINT  C  VS.  PERCEKF  OF 
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FIGURE  2-27.  AT  POINTS  A.  B,  C:  ALUMINA  WEIGHT  FRACTION 
IMPINGED  AS  A  FUNCTION  OF  PARTICLE  SIZE 
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TABLE  2.6 


TABULATION  OF  — *  — 

6W  \p) 

100  .. 

Point 

*^b  (X  of 
Burnout) 

D 

jUii 

(X  of  Burnout) 
in^ 

(d} 

A 

0 

1.4 

254 

2. 54 

A 

0 

2.2 

285 

2.85 

A 

0 

3.5 

168 

1.68 

B 

23 

1.4 

49.6 

0.65 

B 

23 

2.2 

59 

0.77 

B 

23 

3.5 

115.4 

1.5 

C 

58 

1.4 

2.1 

0.05 

C 

58 

2.2 

4.6 

0.115 

C 

58 

3.5 

16 

0.4 

The  next  step  is  to  Integrate  Impingement  over  the  sice  range  of  particles 
found  in  a  rocket.  Kowevert  to  do  this  the  weight  fraction  of  particles  as  a  function 
of  particle  size  be  known.  This  function  can  be  obtained  by  differentiating  a 

curve  of  Sehgal's^*^  data  for  a  500  psl  chamber  pressure.  The  differential  of  his 
percent  total  volume  below  curve  is  shown  in  Figure  2-28  with  units  of  the  ordinate 
as  1/m  since 

6W  <D>  1  “t 

"  dD 

where  is  the  total  weight  of  all  the  alumina  particles. 

By  multiplying  the  ordinate  of  the  curves  in  Figure  2-27  with  the  ordinate 
of  the  curve  in  Figure  2-28  at  each  particle  size  and  plotting  as  a  function  of 
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FIGURE  2-28.  PARTICLE  SIZE  DISTRIBUTION  OF  ALUMINA  SAMPLE  OBTAINED 
FROM  MOTOR  FIRING  AT  500-PSI  CHAMBER  PRESSURE 
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particle  size  while  still  keeping  points  A,  B,  and  C  ms  parameters, 


obtained.  The  ordinate  of  Figure  2*29  will  be 
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Figure  2-29  la 
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For  any  one  of  the  three  points  (A,  B, 
impingement  over  the  entire  particle  size  range, 
expressed  as 


C)  It  is  now  possible  to  integrate 
The  area  under  these  curves  can  be 


00 


O 


which  is  the  weight  fraction  of  alumina  that  impinged  at  any  point.  For  a  run  time 
of  90  seconds »  a  burn  rate  of  500  Ib/sec.g  and  32  percent  of  the  propellant  aa 
aluminag  the  total  impingement  in  Ib/aq,  in.  can  be  calculated  as  a  function  of 
position.  For  example: 


W 

**i 

Impingement  •  ~ 


in 


1 _  1 500  lb. 

2  I  sec . 


(0.32)  (90  sec.) 


W 

~  (14400)  Ib/sq.  in. 

'^T 


This  curve  for  impingement  along  the  nozzle  wall  integrated  for  the  total  firing  is 
shown  in  Figure  2-30.  This  is  the  final  result  desired  from  the  cold-flow  modeling 
test  data. 


2.7.3  100  PM2 

The  model  of  Che  100  FW2  consisted  only  of  the  aft  closure.  Therefore, 
nitrogen  was  used  as  the  osodeling  gas.  As  shown  in  Figure  2-31,  the  aft  closure  was 
modeled  at  two  burn  times,  16  and  54  percent  of  burnout.  Three  particle  aizea  were 
modeled,  making  a  total  of  six  runs. 

Again,  from  Che  modeling  equations  (2.31)  and  (2.32)  the  particle  diameter 
ratio  and  model  pressure  can  be  determined.  Since  the  model  chamber  is  11  Inches  in 

diameter,  a  scale  factor  of  ~  *  9.1  is  dictated.  The  other  values  which  are  input 

Co  the  equations  are  as  follows: 
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FIGURE  2-31,  100  FW2  COLD  FLOW  MODELS 
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For  the  three  particle  sizes  used.  Table  2.7  shows  Che  sices  of  aluaina 
particles  being  tnodeled. 


TABLE  2.7 


PARTICLE  SIZE  RELATIONSHIPS  IN  MODELING  THE  100  FW2  NOZZLE 


Polystyrene 
Particle  Diameter 


Hot  Alumina  Particle 
Diameter  in  Rocket 


Cold  Aluaina  Particle 
Diameter  from  Rocket 


0.8  H 

2.0M 

1.7  M 

1.3  M 

3.3pt 

2.3  M 

2.05jLi 

CM 

lA 

4.4  fi 

For  the  100  FW2,  the  chamber  pressures  are  620  psia  and  643  psia  at  16  and 
54  percent  of  burnout,  respectively.  Therefore,  the  modeling  pressure  at  16  percent 
of  burnout  will  be  S3  psia  and  at  54  percent  of  burnout  will  be  55  psia.  Figure  2-32 
and  2-33  present  the  impingement  data  at  16  and  54  percent  burnout,  respectively. 

The  impingement  data  is  integrated  over  burn  time  and  particle  sizes  found  in  a 
rocket  in  exactly  the  same  manner  as  for  the  gimballed  noccle.  The  final 

result  is  again  the  total  impingement  of  alumina  along  the  aft  closure  as  shown  in 
Figure  2-34. 
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FIGURE  2-32.  PARTICLE  IMPINGEMENT  IN  100  FW2  AT  16%  OF  BURNOUT 
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FIGURE  2-34.  TOTAL  IMPINGEMENT  ALONG  100  FW2  AFT  CLOSURE 
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2.8  PRESSURE  CONTOURS 


Pressure  contours  were  sraasured  on  manoaeCers.  Figure  2*35  is  a 
photograj^  taitan  of  a  pressure  run  with  the  sancssetsrs  stabilised.  Pressure 
contours  were  measured  on  two  of  the  nossles,  (1)  the  gisiballed  nossle  (straight) 
and  (2)  the  gimballed  nozzle  (canted  5  degrees).  Three  pressure  contours  were 
measured  on  the  5  degree  canted  nozzle:  one  on  the  outside  of  the  turn,  one  on 
the  inside  of  the  turn,  and  one  half  way  between, where  no  turn  effects  would  be 
expected. 


Figure  2*36  shows  the  Mach  number  from  these  pressure  contours  by  using 
the  one -dimensional  flow  tables.  These  runs  were  all  made  using  nitrogen  as  the 
modeling  gas.  It  is  of  interest  to  note  that  in  the  straight  position,  the  flow 
did  not  go  sonic  at  the  first  throat,  but  it  did  in  the  case  of  the  5  degree  canted 
nozzle.  This  indicates  nmre  flow  is  realized  when  the  nossle  is  canted. 

These  data  were  used  to  determine  the  convective  heat  transfer  for  the 
gimballed  nozzle  in  section  7. 

The  gimballed  nozzle  in  straight  position  was  axiaymetric.  For  this 
case,  the  mach  numbers  corresponding  to  area  ratios  were  taken  from  the  cne- 
dlmensional  flow  tables  and  compared  to  the  experimental  data.  This  comparison 
is  shown  in  Figure  2-37.  When  sulfur  hexafluoride  was  used  as  the  modeling  gas 
instead  of  nitrogen,  there  was  negligible  change  in  the  Mach  number  contour. 


2.9  FUTURE  WORK  AND  APPLICATIONS 

In  the  future,  sooie  work  should  be  done  to  study  the  basic  parameters 
of  cold  flow  particle  impingement  modeling.  That  is,  the  effects  of  varying 
different  parameters  such  as  turbulence  levels,  boundary  layer  thicknesses,  and 
angle  of  impinging  gas-particle  stream.  Basic  studies  should  also  be  carried 
out  in  the  supersonic  regions  where  particle  impingesmnt  could  be  studied  on  a 
wedge,which  can  also  be  studied  analytically. 

In  the  course  of  the  particle  impingement  modeling,  the  possibility 
of  applying  it  to  other  problem  areas  was  demonstrated.  For  exasiple,  flow 
separation  in  the  nozzle  was  indicated  by  a  step  change  In  particle  impinge* 
ment .  Also,  impingeisent  patterns  on  obstructions  in  the  exit  cone  were  readily 
shown.  Of  course,  nozzles  other  than  the  shapes  studied  could  be  modeled. 
Extream  nozzle  shapes  and  ducting  for  thrust  vector  control,  where  impingement 
might  be  very  high,  could  be  studied. 
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FIGURE  2-35  .  PHOTOGRAPH  OF  MANOiETER  BANK 
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FIGURE  2-36.  MACH  NUMBER  IN  GIMBALLED  NOZZLE,  CAITIED  5  AND  STRAIGHT 
COLD  FLOW  DATA 
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FIGURE  2-37. 


COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  MACH 
NUMBERS  IN  STRAIGHT  GIMBALLED  NOZZLE 
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SECTION  3 


PARTICLE  IMPINGEMENT  IN  SIMULATOR  FIRINGS 


3.1  INTRODUCTION 

Rocket  motor  firings  with  a  calorimetric  molybdenum  nozzle  were 
conducted  with  the  following  objectives: 

(1)  To  study  the  heat  transfer  to  the  nozzle  wall  due 
to  particle  impingement. 

(2)  To  obtain  quantitative  data  on  the  influence  of 
various  concentrations  of  particles  in  the  flow 
on  heat  transfer  to  inlet,  throat,  and  exit  cone 
regions  of  nozzles. 

(3)  To  study  the  variation  in  particle  flow  heat  transfer 
with  nozzle  geometry  (particularly  the  canted 
configuration) . 

The  test  firings  were  conducted  using  a  rocket  motor  in  which  slurries 
containing  various  concentrations  of  solid  particles  were  injected  into  the 
combustor  along  with  hydrogen  and  oxygen  gases.  In  most  tests,  the  solid 
particles  consisted  of  alumina  between  1.0  and  5.0  microns  in  diameter.  Combus¬ 
tion  products  and  particles  were  exhausted  through  a  molybdenum  nozzle 
instrumented  inside  and  outside  with  surface  thermocouples  at  15  locations  to 
obtain  the  desired  heat  transfer  data. 

A  total  of  64  runs  was  conducted.  The  alumina  concentration  in  the 
exhaust  products  was  varied  from  0  co  35  percent.  Throughout  the  schedule  of 
tests  the  formulations  were  designed  for  a  constant  theoretical  chamber 
temperature  of  about  5100  F.  The  principal  results  of  the  tests  are  presented 
in  Figures  3-7  through  3-18, which  illustrate  wall  temperatures,  heat  fluxes, 
and  film  coefficients  for  various  propellants,  solids  concentrations,  and 
nozzle  locations  and  geometries. 
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3.2  DESCRIPTION  OF  TEST  APPARATUS  AND  RUN  PROCEDURE 


The  rocket  motor  and  the  accompanying  molybdenum  nozzle  are  Illustrated 
in  Figure  3-1.  The  slurry-gas  Injection  combination  allowed  a  wide  range  of 
changes  In  combustion  products.  For  example,  with  rearrangement  of  feeds  the 
alumina  content  of  the  combustion  gases  could  be  eliminated  or  varied  while 
allowing  very  little  change  in  gaseous  composition  or  temperature. 

Combustion  conditions  were  regulated  by  metering  the  gases  with  sonic 
orifices,  and  the  slurry  was  metered  by  a  positive  displacement  pump  pushing 
across  a  diaphragm  with  a  hydraulic  oil  output.  The  Injector  was  designed  for 
good  mixing  of  gaseous  and  slurry  feeds  with  a  mlnlmunr.  of  momenta  In  the  direction 
of  the  nozzle.  The  slurry  was  atomized  In  twelve  injector  Inlets  where  a  single 
stream  of  the  slurry  admixed  with  an  annular  stream  of  hydrogen.  This  stream 
Impinged  upon  the  oxygen  for  further  mixing  In  the  combustion  area  at  the  in¬ 
jector  end  of  the  chamber.  The  solid  copper  Injector  readily  accomnodated  the 
planned  firings  without  over-heating.  A  photograph  of  the  simulator  engine  from 
the  Injector  erd  Is  shown  In  Figure  3-2. 

Chamber  wall  design  (graphite  backed  with  asbestos  phenolic)  was 
patterned  after  other  Fhllco  motors  which  have  withstood  total  firing  times 
of  around  1000  seconds,  and  single  firing  durations  of  over  two  minutes. 

Heat  loss  from  the  chamber  was  large,  and  was  accounted  for  In  data  reduction. 

The  step  at  the  end  of  the  chamber  caused  boundary  layer  break-up.  Therefore, 
boundary  layer  conditions  resulting  from  simulator  chamber  operation  would  have 
reduced  Influence  on  the  downstream  flow  along  the  nozzle  section. 

The  greatest  number  of  runs  consisted  of  Injection  of  an  alumina-water 
slurry  into  the  combustion  gases.  Two  types  of  alumina  were  used: 

(1)  Linde  Alumina  Abrasive  -  type  1.0-C  Highly  Uniform, 

1.0  macron  diameter.  (Poor  atomization  forced 
eventual  elimination  of  this  type) 

(2)  Norton  Alumina  Powder,  38-1200  mesh,  2.5  micron 
average  diameter. 

Originally  a  gelling  agent  was  used  to  maintain  the  alumina-water  suspension. 

This  was  eliminated,  however,  because  the  gel  had  an  adverse  effect  on  the  dis¬ 
persion  of  particles  by  the  injector.  A  modified  test  procedure  was  adopted  in 
which  the  ball-milled  suspension  was  added  to  the  run  tank  immediately  before 
firing.  Any  pre-run  settling  effects  were  determined  to  be  negligible. 

Two  types  of  molybdenum  nozzles  were  tested:  an  axisymmetrical  and 
a  S- degree  canted  shape.  Figure  3-3  is  a  photograph  of  the  axisymmetrical 
nozzle.  The  throat  diameter  in  each  configuration  was  1.25  inches.  Figure  3-4 
illustrates  the  distribution  of  thermocouple  locations  on  the  axisymnecrical 
nozzle.  Accompanying  this  diagram  is  Table  3.1  which  is  a  tabulation  of  axial 
distance,  radial  distance,  and  nozzle  area  ratio  of  each  thermocouple  Installation. 
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FIGURE  3-1.  ROCKET  MOTOR  AND  MOLYBDENUM  CALORIMETRIC  NOZZLE 


FIGURE  3-2.  HEAT  TRANSFER  ROCKET  ENGINE,  INJECTOR  END 


FIGURE  3-3 
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SURFACE  THERMOCOUPLE  LOCATIONS 
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TABLE  3.1 


MOLYBDENUM  CALORIMETRIC  NOZZLE  GEOMETRY 
AND  THERMOCOUPLE  LOCATION 


Axial  Distance 

Z.  Inches 

Radius.  Inch 

Area  Ratio 

Thermocoup  le 
Number 

0.418 

2.405 

14.8 

1 

0.618 

2. 119 

11.5 

2 

1.290 

1.159 

3.44 

— 

1.430 

0.991 

2.5 

6 

1.990 

0.667 

1.14 

8 

2.314 

0.625 

1.0 

9 

2.638 

0.667 

1.14 

10 

2.680 

0.680 

1. 18 

— 

3.318 

0.867 

1.92 

12 

3.988 

1.063 

2.9 

14 

4.328 

1.162 

3.47 

15 

Thermocouples  were  installed  on  the  inside  (hot)  surface  of  the  nozzle 
at  the  locations  shown  on  Figure  3-4.  During  testiM,  many  of  these  malfunctioned. 
On  a  typical  test,  data  were  gained  from  four  to  seven  thermocouples. 

Backwall  thermocouples  were  spot-welded  to  the  outside  (cool)  surface 
of  the  nozzle  at  four  locations.  Temperatures  for  other  positions  were  determined 
by  interpolation. 

The  majority  of  the  thermocouples  were  Nanmac  Model  P,  but  a  few  of  the 
same  basic  design  ntade  by  the  Thermotron  Company  were  also  used.  No  difference  in 
the  output  of  the  two  brands  could  be  detected.  Figure  3-5  shows  the  Nanmac 
thermocouple  design.  The  junctions  on  the  Inside  thermocouplea  were  made  flush 
with  the  nozzle  to  preclude  surface  flow  disturbance.  The  junctlo.  was  amde  by 
erosion  or  abrasion  across  the  0.0001-inch  mica  insulation  between  the  two 
0.001-inch  thermocouple  wires  as  shown  in  Figure  3-5.  At  hottest  locations, 
tungsten,  tungsten-rhen  um  junctions  were  used.  Remaining  thermocouple  Junctions 
were  chrome 1-alumel. 

In  a  typical  test,  the  slurry  was  introduced  2.0  to  2.5  seconds  after 
the  start  of  hydrogen- oxygen  combustion.  The  initial  combustion  period  served  to 
preheat  the  wall  surfaces.  The  combustion  restabilized  rapidly  after  slurry 
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FIGURE  3-S.  NANMAC  MODEL  P  THERMOCOUPLE 


Injection  and  the  run  was  normally  continued  for  an  additional  two  seconds. 
Continuous  nozzle  surface  temperatures,  chamber  pressure,  and  flow  rate  data 
were  recorded  throughout  each  run . 

To  summarize  the  test  apparatus  description  and  run  procedure,  a  brief 
tabulation  is  presented  below: 


Engine  size 

6"  I.D.  and  55"  long 

Throat  diameter 

1.7.5" 

L* 

1280" 

Chamber  pressure 

275  psia 

r  *1 

Injector 

fl2  H2  +  12  slurry  1  on  2A  O2 

Propf Hants 

Gaseous 

H2.  O2 

Slurry 

(1)  AI2O3  -  H2O 

(2)  Zr02  -  H2O 

(3)  A1  -  CH3OH 

Exhaust  Solids  Concentration 

0  -  35X 

Corresponding  Solids/ 

0  -  0.55 

Gas  Ratio 


3.3  DATA  REDUCTION  AND  TEST  RESULTS 

Tne  major  step  In  data  processing  was  the  reduction  of  the  transient 
thermocouple  data  to  total  heat  flux  in  the  nozzle.  This  was  done  on  a  digital 
computer  using  a  modified  analytical  technique  presented  in  Reference  3.1  from 
JPL.  The  equation  used  in  the  JPL  analysis  is  for  a  cylindrical  shape.  A 
technique  was  developed  to  convert  the  analysis  to  conical  geooieCry.  Also, 
because  the  JPL  program  assumes  radial  conduction  only,  a  complete  3-dlmeaslonal 
transient  conduction  analysis  of  the  molybdenum  nozzle  was  conducted  for  typical 
run  conditions  as  an  error  check.  A  very  versatile  digital  computer  program 
termed  "thermal  analyzer"  was  used,  which  operates  by  differencing  the  general 
partial  differential  equation  for  conduction  in  solids.  Results  showed  that 
axial  conduction  upstream  and  downstream  from  the  throat  region  of  the  nozzle 
amounted  to  only  about  5  percent  of  the  radial  heat  conduction.  Accordingly, 
it  was  concluded  that  the  JPL  one-dimensional  analytical  approach,  stodlfied 
for  c  mical  geometry,  was  sufficiently  accurate  for  data  reduction. 

The  radiation  heat  flux  was  estimated  according  to  procedures  explained 
in  Section  5,  and  subtracted  from  the  total  to  obtain  heat  flux  due  to  gaseous 
convection  and  particle  impingement.  A  "film  coefficient"  was  then  detcneined 
from  the  latter  heat  flux  and  the  difference  in  fluid  stagnation  and  %«11  surface 
temperatures.  The  stagnation  temperature  used  is  the  theoretical  flasM  temper¬ 
ature  of  the  particular  propellant  combination,  corrected  for  heat  loss  to  the 
combustor  walls. 
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The  data  presented  is  based  mainly  on  the  runs  conducted  during  the 
fourth  quarter  of  the  contract  period.  In  order  to  obtain  proper  comparisons, 
it  vas  necessary  to  select  runs  on  the  basis  of  conuTiuu  thetmocoupie  data.  Some 
of  the  runs,  therefore,  were  not  considered  for  data  reduction  and  presentation 
because  of  Incomplete  thermocouple  recordings. 

Figures  3*6,  3-7  and  3-8  show  direct  temperature  readings  at  thermocouple 
locations  6  (inlet),  9  (throat),  and  12  (exit  cone),  respectively,  versus  time 
elapsed  from  start  of  slurry  injection.  The  slurry  composition  in  these  tests  was 
alumina  and  water.  Note  that  the  wall  temperatures  remain  well  below  the  alumina 
solidification  temperature  of  3700^F.  In  general,  the  temperature  response  is  more 
erratic  as  the  solids  loading  ratio  , 0  ,  is  increased.  The  sharp  peak  in  Figure  3-8 

for  0  ■  0.2  indicates  the  presence  a  local  heat  transfer  phenomenon  of  an 

interesting  transient  nature  due  to  the  particles. 

Figures  3-9,  3-10  and  3-11  show  the  film  coefficients  corresponding  to 
the  previous  temperature  data.  In  the  inlet  and  throat  regions,  there  is  a 
definite  relationship  between  film  coefficient  and  alumina  solids  loading  radio. 

At  the  throat,  the  film  coefficient  increases  by  20  to  30  percent  with  an  Increase 
in  0  of  0.10  to  0.20,  At  the  inlet,  the  increase  in  film  coefficient  with  0  is 
clearly  observed  through  most  of  the  run  times.  The  film  coefficient  increases 
over  the  0  «  0  run  by  40  percent  to  100  percent  when  injecting  alumina  at  0  *  0.30. 

At  the  exit  cone  location  (thermocouple  12),  there  appears  to  be  no 

general  correlation  between  solids  loading  ratio  and  film  coefficient,  as  shown 
in  Figure  3-11.  There  is,  however,  a  very  sharp  momentary  peak  in  film  coefficient 
which  calls  for  further  investigation.  Accordingly,  Figure  i-12  is  presented  which 
shows  the  corresponding  total  peak  heat  flux  at  thermocouple  12.  In  0.10  second, 
the  heat  flux  Increases  from  4.0  to  13.5  Btu/sec  in^ .  It  subsequently  decays  to 
3.5  Btu/  sec  in^  during  the  remaining  0.20  second  of  the  test. 

Similar  rapid  increases  in  total  heat  flux  were  observed  at  inlet  and 
throat  regions  during  other  test  runs,  as  shown  In  Figure  3-13.  It  is  hypothesized 
that  this  response  is  characteristic  of  the  start-up  period  of  solid  propellant 
rocket  operation  when  the  nozzle  walls  are  hot  but  at  a  temperature  less  than  the 
alumina  melting  point.  The  impinging  particles  will  converge  and  solidify  on  the 
wall  surface  until  a  sufficient  depth  of  coating  is  formed  to  cause  a  local  sur¬ 
face  **wave**  of  collected  molten  alumina  to  flow  downstream.  On  flowing  down¬ 
stream,  the  molten  alumina  will  often  pass  onto  an  uncoated  section  of  wall  sur¬ 
face,  causing  an  extreme  temperature  gradient  and  momentarily  higher  heat  flux. 
After  the  alumina  solidifies  on  the  wall  surface,  the  heat  flux  will  decrease 
because  of  an  Insulating  effect.  It  is  probable  that  this  process  has  been  well 
represented  by  the  heat  flux  data  shown  in  Figures  3-12  and  3-13, 

Figure  3-14  is  a  comparison  of  film  coefficients  at  two  circumferential 
locations  on  the  throat  of  a  canted  nozzle.  For  this  run,  aluminum-methanol  slurry 
was  used  to  make  an  overall  propellant  aluminum  content  of  18.8  percent.  In 
describing  the  two  thermocouple  locations,  one  should  picture  the  canted  nozzle 
as  having  one  wall  surface  with  a  generally  larger  radius  of  curvature  than  the 
opposite  wall.  The  thermocouple  on  this  wall  surface  is  therefore  located  at  an 
**outslde  point.**  Conversely,  the  remaining  thermocouple  is  located  at  an  ''Inside 
point."  As  shown  In  Figure  3-14,  the  film  coefficient  at  the  outside  point  is 
about  50  percent  greater  than  at  the  inside  point. 
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FIGURE  3-6.  MEASURED  IMLET  WALL  TEMPERATURES  DtJRlllG  SLURRf  INJECTION 
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FIGURE  3-9.  EXPERIMENTAL  INLET  Flli.  COEFFICIENT  DURING  PERIOD  OF  SLURRY  INJECTION 
FOR  SOLIDS  LOADING  RATIOS  OF  0.0,  O.IO,  0.20,  AND  0.30 


IG1JRE  3-10.  EXPERIMENTAL  THROAT  FILM  COEFFICIENT  VERSUS  TIME  FROM  SLURRY  INJECTION 
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FIGURE  3-14.  THROAT  FILM  COEFFICIENTS  FOR  CANTED  NOZZLE 
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FIGURE  3-15.  Fim  COEFFICIENTS  AT  PEAK  MAGNITUDES  OF  HEAT  FLUX 
VERSUS  thermocouple!  NUMBER 
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FIGURE  3-17.  TOTAL  HEAT  FLUX  AT  ENTRANCE  AND  EXIT  CONE  LOCATIONS  FOR 
ALUMINUM  BURNING  PROPELLANTS  (18.8X  ALUMINUM) 
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SECTION  4 


ALUMINA  PARTICLE  EMISSIVITY  -  EXPERIMENTAL 


4.1  INTRODUCTION 

In  order  to  evaluate  the  radiative  heat  transfer  characteristics  of 
gas-particle  clouds,  knowledge  of  the  particle  emissivitles  at  elevated  temper¬ 
atures  is  required  One  of  the  objectives  of  this  study  has  been  to  determine 
experimentally  the  spectral  emlssivity  of  alumina  at  tenperatures  ranging  from 
the  meltlnp.  point  (2320^K)  to  about  2900^K.  The  experiments  were  conducted 
using  a  carbon-free  system  in  which  a  slurry  containing  only  AI2O3  and  water 
was  adt'ed  to  a  gaseous  hydrogen-oxygen  flame  (details  of  the  experimental 
procedure  are  presented  in  Subsection  4.3). 

Measurement  of  alumina  emittance  was  approached  In  two  different  «rays. 
First,  an  attempt  was  made  to  oieasure  the  radiation  emitted  from  the  gao- 
partlcle  cloud  at  a  location  within  the  chamber  of  the  test  motor.  It  was 
felt  that  the  thermodynamic  state  of  the  system  could  be  more  accurately  de¬ 
fined  at  this  position  due  to  the  lack  of  thermal  lag  between  the  particles  and 
the  gas.  Also,  higher  temneratures  could  be  obtained  in  the  chamber  than  at 
the  exit  of  an  expansion  nozzle.  The  second  method  consisted  of  measuring  the 
emitted  radiation  at  the  exit  of  a  throat  less  motor  operated  so  as  to  produce 
a  hot  subsonic  flow.  This  motor  was  designed  to  produce  approximately  100  pounds 
of  thrust  at  a  chamber  pressure  of  600  psia  with  the  nozzle  In  place;  firings 
made  without  the  nozzle  were  at  atmospheric  pressure . 

In  addition  to  the  alumina  emlssivity  tests,  a  scrl  .  of  runs  was  made 
for  which  the  emission  spectra  of  other  radiators  were  observed.  That  Is,  sub¬ 
stances  such  as  carbon.  Iron  oxide,  and  aluminum  were  added  to  the  flow,  and 
their  effects  on  the  emission  spectra  noted.  The  carbon  was  introduced  both 
in  particulate  form  (carbon  black)  and  as  liquid  hydrocarbons  (UgH^  and  CH3OH), 
while  the  iron  oxide  (Fe203)  and  aluminum  were  Introduced  as  solid  particles  only. 
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The  spectra  resulting  when  carbon  was  added  as  a  solid  differed  markedly  from 
those  produced  when  a  hydrocarbon  was  added,  but  since  little  is  known  of  the 
combustion  processes  involved,  it  is  difficult  to  attribute  the  results  to  a 
specific  phase  of  the  element.  The  effects  of  varying  the  concentration  of 
solid  Fe203  in  the  flow  were  also  shown. 


4.2  THEORY* 

The  measurement  of  the  a.lumina  emittance  is  accomplished  by  measuring 
the  spectral  steradiance  of  the  alumina  cloud  and  comparing  this  with  the  cali¬ 
bration  sceradiancy  of  a  tungsten  lamp.  The  cloud  steradiance  may  be  converted 
to  the  V verage,  or  weighted,  emittance  of  the  individual  particles  through  knowl¬ 
edge  of  the  particle  size  distribution,  nass  flow  rate,  and  gas  properties  at 
the  point  of  measurement . 

For  the  present  study,  measurements  were  made  at  wavelengths  tor  which 
no  gas  emission  bands  exist  (\  "  1.4,  1.78,  and  2.33/i).  At  other  wavelengths, 
account  must  be  taken  of  the  emission  and  absorption  characteristics  of  the  gas, 
and  the  analysis  necessarily  becomes  mere  complex.  An  account  of  the  theoretical 
expressions  for  this  case  is  presented  in  the  Appendix. 

The  steradiancy  of  an  optically  thin,  homogeneous  (in  temperature 
and  concentration)  cloud  over  a  one  dimensional  path  of  unit  area  and  length  L 
is,  at  wavelength 


(r)  N 
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3  ^ 
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where  e  (r )  is  the  size  (and  temperature)  dependent  spectral  emittance  of  an 
indit'idual  particle  of  radius  t.  The  size  distribution  may  be  characterized  by 
a  distribution  functiou  (j)  (r),  where 


<!'  (r) 


so  that  Equation  (4.1)  becomes 


(4.2) 


L  N  (r)  (J)  (r)7Tr  dr 


watts 


cm  str 


^Please  turn  to  page  4-29  for  foJdout  page  of  symbols. 
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The  total  number  of  particles  of  all  sizes  per  unit  volitaw,  N,  Is  t^ven  by  the 
formula 
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where  the  average  weight  w  Is 
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Th«  raittaacc  of  the  individual  particles  is  a  function  '«f  particle 
flisa  as  shovn  by  Eqwttion  (4.3).  In  order  to  obtain  a  clearer  presentation  of 
eiqperisMntal  results,  an  average,  or  weighted,  enittance,  C\  •  Is  defined: 
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E<!uation  (4.3)  becoaws,  with  Equations  (4.4)  and  (4.6) 
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Equation  (4.7)  is  the  bat>is  for  determination  of  €  with  the  measureisent  of  the 
cloud  sterediance,  R.  The  q'tantities  0  and  L  are  obtained  from  flow  rate  and 
geometry  measurements,  p  and  R^  from  calculations  of  the  state  of  the  gas  at 
the  point  of  amasurement,  and  the  size  distribution  integral  from  prerun  size 
distribution  Masurements .  The  actual  data  reduction  techniques  and  equations 
are  presented  in  another  report 

The  methods  of  obtaining  p,  R^,  and  size  distribution  information 
were  as  follows.  The  state  of  the  gas  is  calculated  from  theoretics!  data  modi¬ 
fied  for  heat  transfer  losses.  Figure  4-1  shows  the  theoretical  chamber  temper- 
aturfi  of  Che  H2  -  02  -  H2O/AI2O3  system  plotted  against  oxidizer-to-fuel  (0/F) 
ratio  for  various  values  of  the  parasmters  0  and  0 .  The  flow  is  kept  from 
accelerating  and  thus  causing  an  expansion  temperature  drop  by  lowering  the 
flow  rate  until  the  pressure  drop  along  the  chamber  is  small.  Experimental 
aMMsurUviwiat  of  tha  haat  transfer  losses  was  made  by  measuring  chandaar  coolant 
inlet  and  outlet  temperatures*;  tha  chamber  temperature  drop  due  to  heet  transfer 
then  could  bo  calculated.  Soom  independent  tests  of  this  calculation  were  made 
by  using  ^he  oodium  line  reversal  technique  to  measure  the  gas  temperature 
directly.  The  vesulta  of  both  the  calculation  and  the  measuramant  showed  a 
relative i.jr  small  loss  in  temperature,  about  lOO^K  in  3000^.  Tharafora,  as  a 
smttar  of  routine,  the  chamber  tempereture  was  taken  to  be  the  theoretical 


* 

This  swasuresmnt  also  sufficed  to  indicate  when  steady  state  was 
reached . 


FIGURE  4-1.  CHAMBER  TEMPERATURE  VERSUS  0/F  FOR  K„  -  0  -  H  C/Al-0-  SYSTEM 


t6Dper«tura  less  lOO^K.  This  embodies  the  a88umpt;.ou  that  the  heat  transfer 
dnes  not  change  slgnlflcanMy  when  partlCi.e8  are  addled,  the  chamber  pressure 
meaffureirent  was  mude  for  thi3  gas-^oniy  system*  This  should  be  a  valid  a8Sia4>tion 
sl:ice  the  coolant  temperature  did  not  rise  signlf icontly  between  the  gas-only 
and  the  gas-particle  cas';s,  and  also  because  the  radiative  transfer  could  not 
be  too  much  different  for  the  two  cases  as  can  be  neen  by  examining  the  spectra. 

In  order  to  determine  the  gas  density  from  the  equation  of  state, 
molecular  weight  values  computed  from  an  equilibrium  flow  thevmochemical 
computer  pr'igram  were  used  in  conjunction  with  ^he  measured  static  pressuie  and 
temperature  values  determined  as  described  above.  This  computer  program  also 
produced  the  temperature  data  of  Figure  4-1, 


The  size  distribution  was  examined  In  three  different  way?.  First, 
samples  of  the  unflred  powder  were  run  through  a  Micromerograph  size  analyser 
which  indicated  cumulative  weight  vexeus  size,  and  these  results  were  analyzed 
by  successive  graphical  differentiation  to  get  a  number  histogram.  This  was  in 
turn  normalised  by  dividing  the  total  area  under  the  curve  to  get  ^  (r) .  The 
resultant  distribution  function  pertinent  to  these  testa  is  shown  in  Figure  4-2. 

Secondly,  the  particles  were  collected  by  a  shutter  probe  nlaced  on 
the  exhaust  axis  about  15  fe^t  from  the  exit  plane.  The  collected  samples  were 
observed  under  an  optical  raicrcscope  and  qualitative  Infcrmation  concerning  the 
size  and  shape  of  the  particles  was  recorded.  In  most  cases,  the  collected 
particles  were  spherical  in  shape.*  The  largest  collected  particles  appeared 
to  be  about  10  microns  in  diameter,  which  agreed  with  the  distribution  of 
Figure  4-2. 


The  third  check  of  parctcle  size  distribution  was  obtained  by  measuring 
the  attenuation  of  light  incident  upon  the  particle  cloud  (transmission  tests). 
The  transmitted  lamp  spectral  steradiancy  is 
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which  becomes,  with  Equations  (4.2)  and  (4.4) 
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This  also  provided  an  independent  check  of  the  theoretical  temperature  since 
particles  fired  at  O/F's  giviug  chamber  tempezatures  below  the  melting  point 
of  alumina  were  not  spherical,  but  angular  like  those  put  into  the  slurry. 
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FIGURE  4-2.  PARTICLE  DISTRIBUTION  FUNCTION  VJ2R3US  PARTICLE  RADIUS 


Cc) 

Infonaation  on  Q  (r)  near  the  melting  point  of  elumina  is  available  from  re¬ 
sults  cf  Burcli  and  Gryvnak4.2  used  in  conjunction  with  the  Philco  Hie 
scattering  program.  In  analogy  to  Equation  (4.6),  an  expression  for  a  weighted 
axr faction  efficiency  factor  may  be  written; 
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-and  Equation  (4.9)  may  be  written 
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<4.11) 


The  quantity  Q<^)  is  not  too  sensitive  to  size  distribution  for  the 
large  particles  (geometrical  optics  region)  of  this  vestlgation,  and  thus  the 
noatlnal  values  of  Q  (r>  Sased  on  Burcli's  data  anu  Equation  (4.9)  provide  a 
measure  of 

oc 


71  r  ur 


0 


That  is  to  say,  the  original  distributlcn.  Figure  4-?.,  is  used  with  q(^)  (r)  to 

get  Q  and  Q  is  used  in  conjunction  with  Equation  (4.11)  and  the  ^^xperi- 

mencal  transmission  tu^ixd\xx:^^nt  to  check 
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dr 


In  the  M:ke  region,  this  would  not  be  possiole  as  Q  would  be  highly  depende.it 
upon  the  distt  ibution.  The  actual  data  reduction  ?cheipe  used  in  «.n.' lysis  of  the 
transmission  tests  is  also  described  in  Reference  4.1. 
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Note  that  the  expressions  derived  in  this  subsection  are  actually  for 
throatless  motor  experiments  and  for  the  chamber  radiancy  measurements  in  which 
a  no-lag  condition  exists.  However,  by  proper  application  of  the  gas-particle 
velocity  ratios  to  Equation  (4.4)  the  correct  value  for  the  particle  number 
density,  N ,  can  be  determined  for  the  expansion  norsle  case  in  which  considerable 
velocity  lag  may  exist  between  the  gas  and  the  particles.  Accounting  for  the 
thermal  lags,  which  also  exist  for  this  case,  is  not  trivial,  however.  For  this 
reason  the  experiments  were  conducted  under  the  no-lag  conditions  where  the 
state  of  the  gas  and  particles  could  he  more  easily  determined. 

A  systematic  analysis  of  the  experimental  errors  involved  has  been 
performed  and  details  of  the  results  may  be  obtained  in  Reference  4.1  • 


4.3  EXPERIMENTAL  EQUIPMENT 

During  the  course  of  this  contract,  another  study  was  under  %9ay  at 
Philco  Research  Laboratories,  under  the  sponsorship  of  Advanced  Research  Pre  • 
ects  Agenry  (ARPA)  (Contract  No.  NOnr  36  ’6( 00) directed  toward  further 
understanding  of  radiative  properties  of  elevated  temperatures  (specif¬ 

ically,  above  the  melting  temperature).  The  Interests  of  both  studies  and  the 
equipment  necessary  for  their  completion  were  compatible.  This  allowed  a  larger 
effort  in  the  development  of  the  spectroscopic  equipment  and  the  optical  system, 
and  also  permitted  a  more  extensive  experimental  program.  Some  of  the  results 
reported  here  were  also  reported  in  the  NOnr  3606  final  report .^*1 

The  optical  arrangement  i  led  in  the  infrared  radiation  experiments 
Is  shown  schematically  in  Figure  4-3,  and  the  spectrometer  and  light  pipe  con¬ 
figuration  are  shown  in  Figure  4-4.  The  spectrometer  is  a  Perkin-Elmer  Kodel  108, 
rapid-scan  instrument  utilizing  a  lithium  fluoride  prism  and  liquid  nitrogen 
cooled  indium  antimonide  photoconduct ive  detector.  In  order  to  provide  lower 
scan  speeds  than  the  1-1/2  cycles  per  second  originally  available,  a  variable 
speed  electric  motor  was  connected  directly  to  the  drive  shaft  of  the  spectrom¬ 
eter  internal  urive  motor.  The  original  drive  system  provides  for  14  discrete 
scan  speeds  ranging  from  3  to  150  scans  per  second,  and  the  interchangeable 
variable  speed  motor  provides  rates  from  nearly  zero  to  10  scans  per  second. 

Low  scan  ratei?  improve  tne  data  readout  and  available  resolution  (allowing, 
smaller  slit  openings)  for  long  duration  firings. 

Two  diff«^rent  test  motor  configurations  were  employed  during  the 
course  of  the  study,  one  oeing  used  for  the  chamber  radiation  measureok^nts ,  and 
the  other  throatless  ujotor  previously  mentioned.  Although  the  physical 
dimensions  of  the  first  motor  (chamber  studies)  are  described  In  detail  in 
Section  3  of  this  reporc,  the  characteristics  of  importance  to  the  radiation 
studies  will  be  mentioned  here.  Thf  motor  consisted  of  an  exterior  steel  Jacket 
surrounding  a  graphite  heij  sink  liner  (l.D.  •  6,0  inches).  The*,  chav^ber  view 
pnrtfi  were  arranged  at  the  nozzle  end.  The  position  is  indicated  in  Figure  4-5. 
Tue  chamber  port  and  window  view  are  shown  in  Figure 3-1.  A  single  wlnduw  was 
used  for  a  tots"  i.4ialatlon  thermocouple -type  *^adioineter ,  Preliftinary  spectral 
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views  wore  also  taken  through  a  single  window  with  focal  point  Just  Inside  the 
chamber  wall.  The  setup  for  calibrated  spectral  radiancy  required  focusing  at 
the  center  of  the  chamber  while  the  tungsten  filament  comparison  source  was 
focused  through  the  opposite  window.  Window  diameter  was  0.2  inch. 


The  smaller  engine  and  the  500  mm  Jarre 11-Ash  scanning  grating  mono¬ 
chromator  are  shown  in  Figure  A-6.  Thi.s  spectrometer  was  used  to  make  the 
sodium  line  reversal  measurements  of  chamber  temperature  , and  has  been  described 
elsewhere The  engine  uses  gaseous  hydrogen  and  gaseous  oxygen  (as  does  the 
larger  engine  described  in  the  preceding  paragraph) .  An  H^O/Al202  slurry  can 
be  injected  by  hydraulically  forcing  a  piston  to  squeeze  trie  slurry  Into  the 
chamber  at  the  desired  flow  rate.  The  slurry  flow  rate  Is  measured  in  two  ways: 
by  monitor  ng  the  motion  of  :he  piston,  and  by  recording  the  hydraulic  fluid 
flow  with  a  turbine  meter.  To  insure  good  particle  dispersion,  the  slurry  is 
injected  at  four  locations,  each  one  aligned  with  hydrogen  and  oxygen  injectors. 
Hydrogen  and  oxygen  flow  races  are  obtained  by  measuring  the  pressure  upstream 
of  a  calibrated  sonic  orifice , 


Chamber  pressure  in  the  smaller  engine  is  measured  at  two  locations: 
near  the  Injector  head,  and  near  the  chamber  exit.  Mass  flow  rates  are  adjusted 
so  that  the  chamber  pressure  drop  is  kept  small  and  the  gac  state  is  near  the 
stagnation  value;  l.e.,  a  minimum  amount  of  combustion  energy  has  appeared  in  the 
form  of  kinetic  energy  of  the  exhaust  products.  This  is  desirable  since  the 
data  reduction  requires  knowledge  of  the  particle  concentration  and  temperature, 
and  these  quantities  depend  upon  the  gas  flow  properties.  A  typical  total  mass 
flow  rate  Is  0.4  Ib/sec,  and  chamber  pressures  are  about  18  and  14  psia  for 
upstream  and  downstream  values,  respectively. 


4.4  RESULTS 

Certain  difficulties  were  encountered  while  attempting  to  measure  the 
spectral  emlttance  of  an  alumina-containing  flame  at  a  location  within  the 
chamber  of  the  larger  rocket  motor  Calculations  had  been  made  showing  that 
the  surface  temperature  of  the  heat  sink  liner  would  not  be  sufficient  to  develop 
an  appreciable  amount  of  wall  radiation.  Erosion  of  the  chamber  wall  apparently 
resulted  In  the  formation  of  small  particles  which  were  then  able  to  attain 
temperatures  comparable  to  the  gas  chamber  temperature.  This  was  made  evident 
when  attempts  were  made  to  record  the  emission  spectra  (in  the  chamber)  of  a 
pure  hydrogen-oxygen  flame.  Each  trial  resulted  in  the  observation  of  an 
overwhelming  continuum  emission  which  is  not  characteristic  of  the  H2-O2  flame. 
Since  no  solid  species  were  added  to  the  flow,  one  possible  source  or  trie 
observed  radiation  was  the  presence  of  wall  carbon  particles.  A  visual 
inspection  of  the  liner  revealed  roughness  and  erosion  of  the  material  at 
some  points. 
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FIGURE  4-5  ENGINE  AND  SPECTROMETER  ACCESSORIES 


FIGURE  4-6.  SMALL  SLURRY  ROCKET 


Satisfactory  oporatlon  the  system  for  obtaining  chamber  combustion 
spectra  was  critically  dependent  upon  proper  window  performance.  This  required 
that  the  windows  stay  clean,  intact,  and  that  the  window  purge  gas  not  be  ex¬ 
cessive,  prior  to  the  spectral  sladles.  the  same  window  design  was  used 
to  house  an  asymptotic  calorimeter  for  total  radiation  data.  The  data 
were  too  scattered  to  show  any  reliable  relationship  between  alumina  loadings 
and  total  radiation.  The  causes  of  this  could  only  be  guessed  to  be  such 


;:tuLS  as  exct^bslve  purge  ga s^  lecliculat Ing 


to  the  window,  or  variable  effects  of  the  graphite  wall. 


The  final  opt Ics  system,  designed  to  properly  focus  the  comparison 
source  and  chamber  emission,  produced  a  negligible  amount  of  signal  from  the 
scanning  spectrometer.  The  need  for  much  larger  windows  was  indicated.  This 
factor,  together  with  the  occurrence  of  continuum  from  gas-only  combustion 
and  a  number  of  window  failures,  led  to  the  decision  to  focus  efforts  on  the 
throat  less  chamber  combustion. 


A  series  of  experiments  designed  to  measure  the  cloud  emissivlty  of 
an  alumina-containing  flame  are  described  in  detail  in  the  Appendix.  The  re¬ 
sults  of  these  tests  were  the  basis  for  the  de\elopment  cf  a  new  set  of  experi¬ 
ments  designed  to  determine  the  spectral  emissivlty  of  Al20'j  at  elevated 
temperatures.  The  hoc  subsonic  flow  of  the  throatless  rocket  proved  to  be  a 
valuable  tool  as  nubiy  be  seen  by  inspection  of  the  results. 

Figure  4-7  shows  the  experimentally  determined  spectral  emissivlty 
data  due  to  alumina  as  a  function  of  particle  temperature,  while  Figure  4-8 
shows  the  experimental  results  of  Burch  and  Gryvnak.^'^  The  curve  In  Figure  4-7 
at  the  lower  temperatures  (below  2320^K)  Is  the  result  of  Mie  theory  calculations 
based  on  alumina  optical  properties  measured  by  Burch  and  Gryvnak  and  are  aver¬ 
aged  over  wavelengths  between  one  and  three  microns.  The  brackets  shoim  on 
three  of  the  data  points  in  Figure  4-7  represent  the  root  mean  square  error 
calculated  from  best  estimates  on  measurement  of  data  reduction  and  experimental 
uncertainties.  Examination  of  Figure  4-8  indicates  that  extension  to  1/2  and 
3-1/2  microns  may  not  be  too  bad  for  some  practical  applications.  At  microns 

V  ic  nr^  nvf  (■  1 ^  r%€  lO  a  a  .s  4  o»n  1 1  tj  4  t'  Vi 

attempted  extrapolation.  For  that  data,  however,  one  can  assume  that  c  Is  pro¬ 
portional  to  n**,  or  k,  as  long  as  2  ix  r /X  >>  1  and  n’*<<  n*  to  a  good 
approximation.  Therefore  the  Mie  theory  calculations  of  Figure  4-7  could  be 
extended  to  wavelengths  of  4,  5,  or  6  microns  by  multiplying  the  value  shown 
in  Figute  4-7  at  T,  by  the  ratio  of  k  at  4,  5,  or  6  microns  to  k  at  two  microns 
for  the  same  value  of  T, 


Similar  extrapolation  of  the  data  above  the  melting  point  Is  not  possible 
because  long  wavelength  measurements  were  not  obtained  there.  The  high  temper*^- 
ture  data  of  Figure  4-7  indicate,  however,  that  alumina  approximates  a  graybody 
at  \  -  1-4,  1.73,  and  2.33  ju.  That  is,  the  departures  from  graybody  behavior  are 
so  small  that  they  cannot  be  resolved  within  the  limits  set  by  the  experimental 
error  . 
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Rc^;nrdii\[';  Ihc  bclwivior  oi  the  results,  we  note  that  a  diatonl iiuious 
Jui.ip  in  appears  at  the  x  ltJnp>  point.  Burch  aiul  Gryvnal;  reported  ohr.crving  a 
hrigiitnoss  increase  of  one  to  two  orders  of  magnitude  as  a  portioi\  of  the  alumina 
I'lied  and  formed  a  molten  bead  in  the  optical  path.  (They  could  not  obtain 
1^  for  the  1  id ,  however  |  as  the  bead  thicscness  was  unknown.^  This  observation 
supports  tite  experimental  results  in  Indicating  a  discontinuous  increase  in  emit- 
tance  of  alumina  as  the  phase  changes  from  solid  to  liquid.  Siticc  absorption  ami 
emission  of  elect roiuapnetic  waves  depend  upon  conduction  electrons  present  in  the 
material,  it  can  be  specul  ‘  ed  that  the  number  of  conduction  elect ronr.  in  the 


liquid  phase  is  somewlmt  larger  than  the  number  present  in  th 
qualitative  observation  may  be  compared  with  results  of  Keyes 


solid  phase.  This 
^  who  recorded  a 


step  increase,  by  a  factor  of  13,  in  the  RF  conductivity  in  germanium  as  it 


melted.  It  remains  for  others,  however,  to  explain  the  details  of  the  physics 


of  these  experimental  results. 


The  results  of  the  small  motor  tests  fulfill  the  olijectivcs  of  this 
study.  However,  an  attempt  was  made  to  determine  the  effects  due  to  smaller 
particles  then  those  used  (Norton,  type  38-900).  Three  runs  were  made  using 
Norton  alumina  type  38-1200  (essentially  n  1200  mesh  particle)  and  Linde  type 
1,0-C  nominally  one  micron  particles  as  stated  by  the  distributor).  All  three 
runs  resulti'd  in  no  measurable  continuum  emis.vlon.  The  conditions  of  these  runs 
may  be  found  in  Table  4.1.  Two  possible  explanations  have  been  suggested  for 
the  lack  of  emission  from  these  runs.  The  first  is  that  the  size  dependency  of 
the  (i.ilssivity  nv  ■  be  important  for  these  smaller  particles;  i.e.,  the  emisslvity 
becomes  very  small  for  very  small  particles.  Electron  microscope  photographs 
were  made  of  particles  collected  during  these  three  runs  The  results  showed 
that  the  particles  ranged  from  O.I  to  3  microns  in  diameter,  with  the  majority 
being  IpL  or  less.  The  second  possible  explanation  is  that  the  exit  temperatures 
were  lover  than  expected.  The  photographs  also  showed  that  about  30  percent  of 
the  particles  collected  were  unmelted,  indicating  temperatures  very  near  the 
melting  temperature  of  alumina.  In  view  of  the  excess  water  present  in  these 
runs  (V'^  0..45),  it  may  be  that  the  100®K  temperature  drop  observed  in  the 
earlier  runs  was  not  valid  for  these  particular  tests.  No  attempt  was  made  to 
exp'^rimentally  determine  the  heat  losses  for  these  experiments. 


Finally,  a  series  of  nine  runs  was  made  to  determine  qualitatively 
the  effects  of  adding  solids  other  than  AI2O3  to  the  flow  (namely  C,  Fe203, 
and  Al) .  The  following  paragraphs  and  figures  describe  briefly  the  more  notable 
results. 


Benzene  (C,HJ  (4.12) 

o  o 

The  primary  feature  of  this  run  is  the  appearance  of  a  very  strong 
emission  in  the  4.3MC02  band  region.  There  is  no  measurable  increase  in  con¬ 
tinuous  radiation  in  the  1-3/11  region,  and  the  only  change  in  the  appearance  of 
the  spectrum  after  the  addition  of  CgHg  (other  than  the  CO2  band)  Is  the  reduc¬ 
tion  of  the  emission  in  the  l-3fX  region  (presumably  due  to  the  lower  temperature 
after  the  addition  of  the  benzene). 


Figure  4-9;  Run  5-B,  j 
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(4.13) 


Figure  4-10;  Run  3-8,  K^-O^  / Al-CH^OH 


This  run  also  exemplifies  the  effect  of  adding  carbon  to  the  flow  by 
showing  the  strong  4.3^  CO2  emission  characteristic  of  the  vaporised  species. 

In  strong  contrast  to  the  previous  run,  the  gas^only  emission  spectrum  has  been 
altered  considerably  by  the  appearance  of  a  relatively  strong  continuum  in  the 
1,3/J  region.  This  additional  radiation  may  possibly  be  ascribed  to  either  the 
AI2O3  resulting  from  the  combustion  of  aluminum,  or  to  some  interaction  between 
the  carbon  and  the  aluminum,  or  the  unburned  caibon  ir  the  system. 


Figure  4-11; 


Run  1-B,  H^-0^ 


/ 


Al-CH^OH- 


(4.14) 


The  primary  difference  between  the  conditions  of  this  run  and  the  pre¬ 
ceding  one  Is  in  the  presence  of  solid  carbon  in  the  flow.  The  carbon  was  in¬ 
troduced  in  the  form  of  a  powder  (carbon  black)  added  to  the  slurry  of  CH3OH 
and  A1 ,  Of  primary  interest  is  the  increase  in  radiation  in  t!ir  l-3|Li  renge  with 
the  simultaneous  decrease  in  CO2  emission.  Again,  it  is  not  possible  to  a.'^ertbe 
the  emission  to  any  one  constituent,  but  it  seems  likely  that  the  solid  carbon 
in  the  flow  contributed  considerably  to  the  observed  changes. 

/ 

Figure  4-12;  Run  10-B,  H2-O2  /  H^O-Al^O^-C  (4.15) 


For  this  run  carbon  black  was  added  to  an  AI2O3-H2O  slurry,  and  the  re¬ 
sults  are  shown  in  the  figure.  The  continuum  emission  level  (l‘3/i)  produced  by 
this  run  was  so  intense  that  the  recording  system  was  overdriven  and  complete 
details  of  the  spectrum  could  not  be  obtained.  It  is  interesting  to  note,  how¬ 
ever,  that  at  the  saturation  level,  the  emission  was  at  least  five  times  that 
from  the  gas.  Unburned  carbon  is  again  believed  to  be  the  strong  continuum 
emitter.  Another  feature  worth  noting  is  the  complete  absence  of  radiation  in 
the  4.3p  CO2  band  region.  This  spectrum  is  typical  of  those  obtained  when  measure¬ 
ments  were  atten^ted  in  the  chamber  of  the  larger  motor. 

Figure  4-13;  Run  12-B,  H«-0«  / (4.16) 

A  slurry  mixture  consisting  of  AI2O3,  H2O,  and  0.5  percent  by  weight 
of  ^^203  was  added  to  the  H2-O2  flame  during  this  run.  The  amount  of  Fe203 
corresponded  to  1/2  percent  by  weight  of  the  total  mass  of  the  exhaust.  Of 
primary  interest  is  the  appreciable  increase  in  emission  in  the  1-2. 3pt  region 
and  the  absence  of  any  measurable  change  beyond  the  2.3^  point.  This  is  some¬ 
what  similar  to  the  pure  AI2O3  emission,  although  the  increased  emission  in 
th^s  vase  cannot  be  ascribed  with  certainty  to  either  constituent  alone. 
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FIGURE  4-10.  SAMPLE  SPECTRA  FROM  RON  3B 
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FIGURE  4-11.  SAMPLE  SPECTRA  FR(M  RUN  IB 
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FIGURE  4-12,  SAMPLE  SPECTRA  FROM  RUN  lO-B 
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FIGURE  4-13.  SAMPLE  SPECTRA  FROM  RUN  12-B 
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(4.17 


Figure  4-14;  Run  11-B,  /H20-Al202-Fe20, 


C  AWp  U  J.vm  V/  X 


Conditions  for  this  run  were  similar  to  the  preceding  one  with  Che 

7a. n  TKi-fa  r'nno'icftftH  nf  rkn«k  nprrAnf 

.X I IV  X  t€»  A  xvrw  V»  X»*^<  »•:>  A.  •»«.««»«-  —  —  —  —  —  —  ^  ... - 

FC2O3  by  weight  of  the  exhaust  products.  The  resulting  Increase  in  emission 
may  he  easily  seen.  In  this  instance,  the  change  in  radiancy  level  at  the 
2  ,33/Li  wavelength  is  also  more  apparent. 


4  5  CONCLUSIONS 

Based  upon  the  experiments  performed  during  the  coui se  oi  these  inves- 
tigations,  the  following  conclusions  are  drawn 


4.5-1  SPECTRAL  EMITTANCE  TESTS 

Emittance  of  liquid  phase  alumina  particles  ranging  in  diameter  from 
1  to  10  microns  is  of  the  order  0,01  to  0,1  in  the  1  to  3  micron  wavelength 
region  for  temperatures  between  2300  and  2900^K- 

Based  upon  results  of  Reference  2,  the  tmittance  of  solid  AI2O3 
particles  at  the  same  wavelengths,  and  of  the  same  size,  varies  from  about 
5  X  10’^  at  1500°K  to  1  2  X  at  2300°K, 

A  discontinuous  jump  in  spectral  emittance  in  the  1  to  3  micron  wave¬ 
length  region  occurs  as  alumina  changes  from  the  solid  to  the  liquid  phase. 

Significant  increases  in  alumina  cloud  emission  are  produced  by  the 
addition  of  irinor  amounts  of  carbon  or  iron  oxide  impurities  such  as  could  occur 
from  the  burning  of  solid  propellants 


4,5 >2  CLOUD  EMITTANCE  TESTS 

In  small  size  rockets  the  particle  impingement  on  the  nozzle  is  so 
serious  that  the  simulation  of  the  rocket  plume  of  a  large  motor  is  not  achieved 

The  throatless  H2Cg)-02tg)  rocket  with  aqueous  slurry  provides  poten¬ 
tially  useful  devices  for  determining  the  radiation  emitted  by  AI2O3  particles 
present  in  rocket  plumes.  In  particular,  it  allows  tests  to  be  conducted  at 
high  particle  temperatures  and  without  dynamic  and  thermal  lag  between  particle 
and  gas  - 


A  technique  was  established  to  determine  plume  radiation  properties 
by  first  measuring  with  only  the  fuel  and  oxidizer  burning  and  then  with  the 
particulate  substance  present  without  interruption  of  the  firing.  This  pro¬ 
cedure  provides  an  excellent  demonstration  of  the  effect  of  the  addition  of  the 
part  ides  , 
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FIGURE  4-14.  SAMPLE  SPECTRA  FROM  RUN  11-B 
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The  radiation  einltte:^  by  the  plume  of  a  small,  100-pound  thrust  rocket, 
with  one  micron  AI2O3  particles  added  was  small  compare  to  the  radiation  emitted 
by  products  (mainly  H2O  at  atmospheric  pressure)  of  the  burning  of  the  fuel  and 
oxldlser  alone. 


4.5.3  CHAMBER  RADIATION  TESTS 

The  use  of  a  graphite-lined,  heat-sink  chamber  is  not  a  suitable  instru¬ 
ment  for  studying  the  isolated  radiative  characteristics  of  materials  such  as 
alumina . 


The  view  window  design  employed  is  not,  In  general,  satisfactory  for 
hlc;h  pressure  experiments  of  the  type  attempted.  A  sturdier  construction  should 
be  devised  to  prevent  damage  by  overpressures,  and  larger  limiting  aperture^  would 
be  desirable. 
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NOMENCLATWr; 


CransmliLed  apecLral  Irradlance 
incldani  spectral  liradlance 

ipectral  abaorptlon  coefflclei.t  per  unit  length 
flo'DC  pathlength 

tcmpet eturc  ^iid  wavelength  -  dependent  real  part 
of  the  index  ot  refracilon 

temperature  and  wavelength  -  dependent  taiaginary 
part  ot  the  Index  of  refrectlon 

number  o t  particles  per  unit  volume  In  a • re  range 
between  r  and  r  +  dr 

total  numbei  of  particles  ot  all  alzea  per  unit. 

\'o  1  ume 

chamber  prebsure  (measured) 

apectral  efliclency  tactor  for  extiulctlon  (ratio 
04  extinction  cross  section  to  geometrical  cross 
section^  for  particles  In  size  raniie  between  t  and 
r  4  dr 

spectral  average  extiiictlon  etflclency  factor 
radius  oy  spherical  particle 
spectra!  steradlan^y 
temperatute 

weight  of  particles  In  size  range  between 
r  and  r  dr 

average  weight  for  particle  sire  distribution 
(see  Equation  it)) 

geomet  riisl  dimension 

•nsss  flow  rate  of  propellant 

volurn  concentration  of  pa  tides  fn  aqueous 
suspension 


GREEK 

>  Spectral  extinction  cuefflclent 

t  (r)  spectrsl  ealtt^nce  of  alumina 
parllclas  in  ilse  range  br* 
tween  r  and  r  dr 

T  avarag*  •pectral  ^..rtttanca  for 

pat'tlcla  alza  distribution 
(see  Equation  (b)} 

^  wavelength 

f-  gas  density 

T  sp^fCtral  transmittanca  of 

particle  cloud 

0  al'jmlna-to-gas  mass  ratio  In 

j.aa-par t  ic  1«  system 

d  (r)  particle  size  distribution 
function  (ace  Equation  (2)) 

V  aluoina  mass  fraction  In 

slurry 

SUBSCRIPTS 

\  wavelength 

o  refers  to  coinparlson  source 

radiancy 

g  refers  to  gas  value 

p  refers  to  part  id*,  value 

a  absorption  (l.e.7pe  rafars  to 

vo4.uiretrlc  absorption  coafftclent 
for  particles) 

t  total 

9  scattering 

c  refwra  to  carbon  value 


SUPERSCRIPTS 


indicates  blackbody  value 


SECTION  S 


ALUMINA  CLOUD  EMISSIVITY  -  APPLICATION  TO  DESIGN 


5.1  RADIANT  HEAT  FLUX  FROM  A  TWO  PHASE  MIXTURE 

Calculation  of  the  radiant-heat  intercharige  between  an  absorbing  gas 
with  suspended  particles  and  its  bounding  surface  is  very  difficult.  This  is 
mainly  due  to  the  complicated  nature  of  the  total  emissivity  of  the  mixtuie. 

The  spectrum  of  an  absorbing  is  composed  of  discrete  absorption  bands  with 

the  total  emissivity  (integrated  over  all  wave  lengths)  being  a  function  of 

(1)  static  pressure,  (2)  mole  fraction,  (3)  temperature,  and  (A>  path  length. 
Whereas,  for  a  particle  cloud,  the  spectrum  is  similar  to  that  of  a  solid  sur¬ 
face  and  gas,  with  the  total  emissivity  a  function  of  (i)  static  pressure, 

(2)  temperature,  (3)  mole  fraction  of  gas,  (4)  weight  ratio  ot  particles  to  gas, 
(5)  average  weight:  of  particles,  (6)  distribution  of  particle  size,  and  (7)  path 
length.  Therefore,  due  to  the  difference  in  the  IR  spectrum  of  particles  and 
gas,  direct  addition  of  total  emlssivities  to  obtain  the  total  emissivity  of  the 
mixture  is  not  Justltied.  That  Is,  the  spectral  emissivity  obtained  by  direct 
addition  of  the  absorption  coefficients  for  the  gas  and  particle  cloud  is  valid 
only  if  the  gas  and  particles  do  not  influence  each  other  in  their  emission 
(Reference  5.1). 

Ir.  a  2-phase  mixture,  the  energy  emitted  by  the  gas  will,  in  part,  be 
scattered  by  the  particles  Therefore,  the  spectral  energy  emitted  from  an 
Incremental  volume  of  gas  plus  particles,  per  unit  cioss  sectional  area,  and 
normal  to  a  point  (0)  on  the  surface  of  the  media’s  bounding  surface,  may  be 
written  as 


Please  turn  to  page  5-21  for  foldout  symbol  pages. 
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That  Is,  the  total  radiation  from  a  volume  element  of  the  cloud  is  the  sum  of 
(1)  the  individual  radiating  species, and  (2)  the  energy  reflected  or  scattered 
by  the  particles  which  originate  at  locations  adjacent  to  the  volume  element. 
Since  the  incremental  gas  volume  is  optically  thin,  the  emisslvity  of  the  gas 
and  particles  may  be  written  as 


Xg 

3  . 

»  7  d  X 

'gx 

(5.2) 

xp 

a 

»  -y  cl  X 

'px 

(5.3) 

Multiplication  of  the  above  emissivities  by  the  appropriate  black  body  energies 
and  Including  the  scattered  energy  by  a  defined  coefficient  (7p^)  times  the 
appropriate  black  body  energy  yields  Equation  (5  .1). 

The  amount  of  energy  emitted  from  the  volume  element  that  reaches  a 
point  (0)  is  the  amount  transmitted  throvigh  the  cloud  from  x  to  0. 


where 


r 


Xm 


X 


'tX 


Multiplication  of  Equations  (5.1)  and  (5.4),  and  integrating  between  0  and  £ 
results  in  the  directional  energy  incident  on  a  point  (0)  from  a  nonisothermal , 
hv2terogeneou8  cloud: 


+ 


(5.5) 


For  the  particular  case  of  an  isothermal-homogeneous  cloud,  Equation  (5.5)  reduces 
to 


gX 


7  + 

PX 


7 


si 

PX 


.,0 

Ik 

X 


-\x" 


(5.6a) 


'X 


-^tx^ 


(5 .6b) 
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end  the  emisslvity  of  the  cloud  .Is 


acsuntlng 


'?oi 


p1 

'PX 


I 

[■ 


-  e 


(5.7) 


82 


This  Ir  'he  caje  when  7 
effect  is  rmall  co.apare 
and  Equation  (5  .'>)  becomes 


PX  ^  ^rX  order  of  75^.  When  the  scattering 

a  to  «.mssiou,  then  7pj^  is  neglected  in  Equation  (5.1) 


'Xtn 


V®  + 

&?L .  \£h 


'\x 


1  -  e 


■\x^ 


(5.8) 


Note  that  Equation  (5.8)  Is  used  only  whfrn  multiple  scattering  Is  neglected. 

Also  the  term  is  included  in  the  tran8mi.jslon  term  due  tc  the  inlluence  of 

patii  length.  1?  'V??«  the/i 

PX  ^PX  I  i 

^  ^  1  ii 


'PX 


X  m 


^8X/ 


(5. Pa) 


As  noted  before,  the  above  equations  assume  that  the  emission  of  the  species  in 
the  cloud  do  not  influence  one  aiiorher. 


The  spectral  emissiviry  of  en  opMcally  thin  particle  cloud  (neglecting 
gaseous  emission)  was  discussed  in  detail  in  Section  A.  By  applying  the  same 
physl'tal  definitions  and  empirical  relations  as  was  done  la  Equation  4.7  of 
Section  4,  an  express?  m  for  an  optically  thin  particle  c?oud  may  be  written: 


£ 


Xp 


V 


TTr*"  d  t 


(5.9) 


where  the  particle  density  (number  of  particles  per  unit  volumti)  is  given  as 
N  »  0pg/W;  fg  ^  PM/^T;  M  is  the  gas  molecular  weight  and  0  ■  Pp/Pg 

*  Up .  (Assumed  no  thermal  or  velocity  lag  exists  between  gas  and  particles.) 

In  t!ie  experiments  discussed  in  Section  4,  the  emissivity,  ,  which 
represents  the  errisslvity  of  an  average  size  particle,  was  found  by  use  of 
Equation  (5.9).  Th^  Integral 


TTr 


d  r 


0 


for  the  experiment  was  found  to  be  6.5  x  10 
i  *  2 . 0  inches . 


cm 


with  W  **  5.6  X 


10  gm£ 


and 
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The  exper iment s  were  1 Iml ted  to  atmospheric  pressure  and  only  the 
emisslvity  of  the  particles  at  wavelengths  of  l.A,  1.78  and  2.33 /a  was  reported, 
Tliese  three  wavelengths  of  AI2O3  emission  correspond  to  regions  where  H2O  has  no 
absorption^  thus  justifying  the  A I2O3  spectral  emissi vi cy . 

Froni  the  data  and  experimental  accuracy  given  in  Section  4  it  was  found 
that  the  Sf  lection  of  a  constant  between  I  and  3  /i  is  feasible .  Since  the 
experiments  are  only  applicable  between  1  and  3  ju  and  data  are  nonexistent  at 
other  wavelengths,  one  must  assume  between  0  and  1  jj.  and  between  3 and  »  In 
order  to  determine  the  total  emissjvity.  To  be  conservative,  ex.  was  chosen  as 
1,  between  3  /x  and  <w  and  equal  to  ^1-3  between  0  and  1  /Li.  With  this.  Table  3.1 
was  constructed  and  shows  *^1.3  nnd  ^total*  Comparing  the  two  emissivities , 
one  can  see  the  importance  of  the  unity  emissivlty  assumption  at  low  H ,  ev^en 
chough  the  amount  of  black  body  radiation  between  0  ana  1  /^and  3/i  and  00  is 
small  (see  7r*BB  In  Table). 


TABLE  5.1 


SAMPLE  ALUMINA  PARTICLE  CLOUD  TOTAL  EMISSIVITIES 


FOR 

alumina  in  the 

LIQUID 

PHASE 

Temperature 

1°R) _ 

n 

C 

7.BB* 

(0-1) 

7oBB 

(3-00) 

'■1-3 

^tocal 

4230 

0.438 

X 

lo'O 

13 

19 

0.0007 

0.017 

4500 

0.438 

IC^O 

16 

17 

0.0010 

0.015 

4800 

0.438 

X 

lolo 

20 

14 

0.0023 

0.014 

5025 

0.438 

X 

loio 

23 

13 

0.0041 

0.015 

5190 

0.438 

X 

10^^ 

25 

1 

0 

0.0050 

0.015 

4230 

4.38 

X 

10'" 

,0.025 

0.20 

4500 

4.38 

X 

10^^ 

0.035 

0.19 

4800 

4.38 

X 

lOlO 

■  0.076 

0.20 

5015 

4.38 

X 

10^0 

0.13 

0.24 

5190 

4.38 

X 

10^0 

0.16 

0.25 

4230 

43.8 

X 

■1  r\  1  n 

1  ^ 

0.22 

0 . 37 

4500 

43.8 

X 

’OlO 

0.30 

0.42 

4800 

43.8 

X 

lOlO 

0.55 

0.61 

5025 

5190 

43.8 

43.8 

X 

X 

lolo 

lOlO 

1 

r 

0.7S 

0.82 

0.79 

0.84 

4230 

120 

X 

lo'O 

0.50 

0,59 

4500 

120 

X 

lOlO 

0.63 

0.69 

4800 

120 

X 

lOlO 

0.89 

0.90 

5025 

12c 

X 

lOlO 

0.98 

0.99 

5190 

120 

X 

lOlO 

1 

0.99 

0,99 

*  %BB  (O-l):  Percent  black  body  radiation  in  the  wavelength  range  01  0  to  1 ^  . 
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with  Che  variation  of  Z\  with  temperature  and  wavelength  known,  Equa¬ 
tion  (5.8a)  was  employed  to  determine  Che  spectral  emlssivlty  of  the  cloud. 

(Note  that  7^^  Is  zero  at  the  three  wavelengths  chosen  for  analysis.)  Calcula¬ 
tion  of  ^totai  obtained  from  the  basic  definition  of  emisslvity;  that  is 


e 

T 


0  4 

dX  /  a  T 


dX  /  OT^ 


+ 


d  X  /  O 


d  X  /  0 


where  monochromatic  emiasivity  may  be  taken  out  of  the  Integral  when  constant 
In  the  specified  wave  band. 

For  temperatures  below  the  melting  point  of  AI2O3,  the  data  of  Burch5-2 
was  used  and  the  variation  of  with  X  assumed  similar  to  that  of  K\  (K^  is  the 
absorption  coefficient  obtained  from  experiment  with  data  between  0  and  6  pi  }  .  At 
wavelengths  larger  than  6  can  readily  be  assumed  equal  to  1  based  on  the 

Burch  and  Gryvnak  correlations. 

Fig  are  5-1  is  a  plot  of  the  total  alumina  cloud  emlssivlty.  The  total 
emissivity  is  plotted  versus  temperature  for  constant  Sf  Where 


Ni 


7.57  X  10 


11  0P  M  i 


(Particles/f t  ) 


and 


P  psia 
i  «  ft 


In  order  to  compute  the  total  AI2O2  particle  cloud  emissivity  for  a 
given  combustion  chamber,  one  must  know  the  particles- to-gas  weight  ratio  (0), 
pressure  (P),  temperature  (T) ,  and  gas  molecular  weight  (M) .  Then  assuming  W 


and 


0  (r)  TT  r^  dr 


equal  to  Reference  5.2,  Figure  5-1  may  be  used  to  find 


^rotal * 
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FIGURE  5-1.  TOTAL  Al^O^  PARTICLE  CLOUD  EMISSIVITY 
(TRAKSPARENT  GAS)  VERSUS  TEMPERATURE 
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5.2  CALCULATION  OF  THE  RADIANT  HEAT  FLUX  IN  A  NOZZLE 


The  radiative  heat  fluxes  are  to  be  calculated  for  various  parts  of  the 
100  FW-2,  large  gimbaled,  and  simulator  nozzles  with  emphasis  on  the  following 
cases  . 


(1)  Inlet  radiant  heat  flux  to  the  100  FW-2  nozzle. 

(2)  Radiant  heat  flux  throughout  the  gimbaled  nozzle . 

( 3)  Radiant  heat  flux  throughout  the  t ransient  simulator 
nozzle . 

Therefore,  in  this  subsection,  the  purpose  is  to  devise  an  approximate  method  to 
Compute  the  radiant  heat  flux  for  the  above  cases  using  the  results  of  subsec¬ 
tion  5.1. 


An  analytical  method  for  calculating  the  radiation  heat  flux  to  surfaces 
adjacent  to  the  nozzle  and  produced  by  rocket  exhausts  containing  alumina  particles 
has  been  developed  by  Morizumi  and  Carpenter  (Reference  5.3).  The  exhaust  plume 
was  essentially  Isobaric  and  at  a  temperature  below  the  melting  point  of  alumina. 

By  dividing  the  plume  into  segments  and  devising  an  expression  for  the  effective 
emissivity  and  temperature,  the  approximate  radiant  heat  flux  from  the  plume  cloud 
may  be  found  by  a  numerical  solution.  Calculation  of  the  radiant  heat  exchange 
in  a  rocket  nozzle  is  more  romplicated  due  to  the  nonlsobaric  flow,  wall  irradia¬ 
tion,  and  complicated  geometry.  Therefore,  a  more  simplified  approach  to  the 
problem  was  selected. 


The  radiant  heat  exchange  inside  an  enclosure  with  gray  walls  and  filled 
with  an  absorbing  medium  is  most  easily  found  by  using  the  network  Mthod  (refer¬ 
ence  5.4  and  3.5).  That  is,  the  radiant  heat  flux  to  an  Incremental  area  can 
be  written  as 


where 


4  1  -  e. 
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emissivity  of  absorbing  media 
transmissivity  of  absorbing  media 


(5.10) 


(5.11) 
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When  either  the  al sorbing  medium  has  nonunlfoim  emission  and  transmission 
or  the  enclosure's  shape  is  unsymnetr ical  with  respect  to  Equation  (5  11)  Is 
almost  Impossible  to  integrate.  Also  the  coefficient >  for  an  absorbing  gas 

with  suspended  particles  may  be  a  very  complicated  function  of  /3 ,  ^  ,  and  0  .  For 
examole,  in  a  rocket  nozzle  the  combustion  products  are  undergoing  acceleration  and 
thus  the  pressure  and  temperature  (depending  of  course  on  the  specific  heat  ratio) 
are  a  function  ot  axi*!!  position.  Since  ^  P  **  ^tid  “  Y  ^  "  for  the  absorping  media 
are  strong  functions  temperature  and  pies‘iure,  integration  i.*»  made  impossible  by 
conventional  means.  However,  for  an  enclosure  with  a  sitnpie  shape  and  absorbing 
media  with  uniform  emission  to  a  local  area  A  Equation  (5,11)  may  be  simplified 
considerably.  By  inspection  of  Equation  (5. it)  one  may  recognize  the  equality. 


cos  0  0 

_ I _ 2 

7TS  2 


d  A  -  A-  F, 
n  1  1  -n 


(5.12) 


Also,  let  represent  the  equivalent  emissivity  of  the  absorbing  media;  i  e., 

—  sin  cos  jS  d  )3  d  ^  (5.13) 
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€  m 
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If  the  absorbing  media  consists  only  of  an  absorbing  gas,  Equation  (5,13) 
may  be  recognized  as  the  gas  emissivity  that  is  plotted  in  standard  text  books 
(reference  5.6  and  5.7),  That  is 


a 

=  1  -  e  » 

Defining 

Equation 

-7  L® 

7«  3  ,  equivalenc  cranstnissivicy , 

(5.10)  and  (5.11)  may  be  rewritten  as 

and  assuming  isothermal  conditions. 

e,  r 

*'  -1  IS  -  ^ij 

(5.14) 

k 

F.  T.  J.  -  J 

i-n  1 -n  in  T 

:j,  ■  E  "  (5  15) 

1  m. 

n»l 

, 

With  the  absorbing  media  consisting  of  an  absorbing  gas  with  suspended 
particles,  the  specification  of  the  combined  (particles  and  gas)  emissivities  and 
transmissivities  can  only  be  based  on  experimental  knowledge  (previously  discussed). 
The  presence  of  particles  also  increases  the  number  of  variables  which  effect  € 
and  Therefore,  additional  assumptions  are  required  before  Equations  (5.l4) 
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and  (3.15)  can  be  applied  to  the  three  cases  discussed  previously.  These  as* 
sumptions  are  presented  below  In  detail. 

(1)  When  particles  are  suspended  in  an  absorbing  gas^  the  mono* 

chrnmar Ic  emiaston  of  the  resulting  mixture  may  be  expressed 
according  to  Equation  (5.7),  The  transmissivity  has  been 

Investigated  for  an  AI2O3  cloud>and  it  has  been  found  that 

7  *  The  product  (J^  -  J^)  Is  small  in  a  rocket  nozzle 

which  has  essentially  Isothermal  walls.  That  Is,  when  the  walls 
of  a  nozzle  are  approximately  Isothermal,  and  radlosltlea  (J)  of 
neighboring  Incremental  areas  are  nearly  equivalent  for  large 
values  of  the  view  factor,  It  may  be  assumed  that  the  first  term 
in  Equation  (5.15)  Is  small  and  may  be  neglected. 

(2)  The  emissivity  of  a  particle  cloud  is  expressed  In  Equation  (5.7), 

alteration  of  which  will  approximately  give  the  total  cloud 
emissivity.  Assuming  ^  \  “  .  the  total  emissivity  of  the 

particles  may  be  found  under  a  variety  of  nozzle  conditions  by 
use  of  Figure  5-1.  Using  the  curves  of  McAdams,  the  total 
emissivity  of  various  absorbing  gases  present  in  the  particle 
cloud  may  be  found.  Combining  the  particle  and  gas  emissivity 

by 

e  -  1  -  ;  e  -  1  .  (5.14)  (5.15) 

H  o 

and  “  1  -  e'  (5.16) 


the  mixture  emissivity  is  found.  As  was  disec  >ed  In  Section  5.1, 
this  assumes  (1)  the  particles  and  gas  do  not  influence  each  other 
in  emission,  and  (2)  it  is  possible  to  add  the  total  absorption 
coefficients,  just  as  it  is  possible  to  add  the  monochromatic  co¬ 
efficients  for  noninteracting  substances.  The  combined  emissivity 
as  given  In  Equation  (5 . 16) ,  neglects  the  scattering  coefficient 
included  in  Equation  (5.7).  In  the  experiments  performed  during 
this  contract  period,  emphasis  was  maintained  on  the  emission 
properties  of  clouds.  Time  was  not  available  for  an  intensive 
investigation  of  the  scattering  coefficient  7  ^  and  the  dependence 
of  path  length  on  the  validity  of  Equations  (5.1A)  and  (5.15). 
Therefore,  it  is  proposed,  at  a  later  time,  to  concentrate  interest 
on  the  phenomena  of  light  scattering  by  a  two  phase  mixture 
approaching  optically  thick  conditions. 


(3)  To  eliminate  the  complication  of  a  nonisothermal  absorbing  laedluni 
in  Equation  (5.11),  it  is  assumed  that  the  static  temperature  of 
the  media  at  the  axial  position  of  A^  represents  the  average  temper¬ 
ature  at  which  the  integrated  emission  of  the  nonlsothera&al  gas  to 
the  point  in  question  is  evaluated.  This  assumption  is  reasonable 
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it  the  view  factor  from  A,  to  the  gas  volume  at  is  large, 

add  the  emission  {<--)  from  this  gas  volume  is  noL  sufficiently 
different  from  nelgnboring  volumes. 

(4)  The  pressure  for  evaluating  the  emlssivity  of  both  the  particles 
and  gas  is  taken  at  the  static  pressure  of  A  .  The  justification 
of  this  choice  is  similar  to  that  of  the  static  temperature. 


V  / 


TK.* 

^  t 

is 


mean  beam  length  at  which  both  the 
evaluated  will  depend  on  th^-  axial 


_ _ _  -  I 

L  11..  Ltr 

location 


Slid  ^as  e-iii  i  s  i  V  1 1 y 
in  the  nozz I e . 


(a)  Convergent  Section:  Her-  the  length  is  assumed  equivalent  to 
that  of  a  semi  “Inf  Ini  te  t:yclinder  with  Irradiation  to  the 
center  of  the  base.  This  choice  is  based  on  the  fact  that  as 
the  inlet  angle  of  the  nozzle  approaches  90  degrees  the  con¬ 
vergent  sec  t ion  will  become  a  semi-infinite  cylinder.  For 
Isrge  inlet  angles,  this  assumption  seems  reasonable . 

(b)  Throat  :  A  path  length  is  chosen  which  is  1.4  times  the 
diameter  of  the  throat.  This  is  based  upon  the  irradiation 
to  a  point  on  the  curved  surface  of  an  infinite  cylinder,  the 
volume  of  which  represents  the  approximate  absorbing  media 
radiating  to  the  throat . 

(c)  Divergent  Sec t ion :  In  the  divergent  section,  the  mean  beam 
length  is  assumed  equivalent  to  the  geometry  described  in 
(b),  but  the  diameter  is  taken  as  the  nozzle  diameter  at  the 
point  in  qu e s r i on . 

(6)  In  using  Equations  (3,14)  and  (5.15)  to  compute  the  heat  flux  from 
a  particle  cloud,  it  must  be  assumed  that  the  particles  are  uni¬ 
formly  distributed  and  the  distribution  of  large  and  small  particles 
is  constant  for  a  given  nozzle  cross  section.  The  particle  density 
found  from  Assumptions  (3),  (4),  and  (5)  have  already  Included  this 
phenomena . 


(7)  Assuming  chat  the  particles  in  a  nozzle  do  not  cross  the  gas  stream 
lines  is  in  essence  equivalent  to  Assumption  (6).  Howev'er ,  the 
effect  of  this  assumption  upon  the  accuracy  of  the  calculation  Is 
largely  dependent  upon  nozzle  geometry  and  conditions.  In  very 
large  nozzles,  particle  slippage  is  greatly  reduced. 

(8)  Along  with  Assumptions  (6)  and  (7),  it  must  be  assumed  chat  no 
thermal  or  velocity  lag  exists  between  the  particles  and  gas.  This 
assumption  is  feasible  in  the  convergent  section,  but  In  high 

acce lerat ion  regions ,  one  may  question  its  applicability. 
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Using  the  above  assumptions,  Equations  (5.14)  and  (5.15)  become, 


'll 


01 


(5.17) 


Wlien  calculating  the  radiant  heet  ilux  to  a  black  wall,  one  may  eliminate  the 
assumption  of  a  gray  absorbing  media  (momochromatlc  emlsaivity  is  independent  of 
wave  length)  by  ubing 

q  -  E  -  a  E,  (5.18) 

^1  T  m  T  1 

where  is  evaluated  at  the  media  temperature  end  at  the  surface  temperature 

The  radiative  heat  transfer  coefficient  is  defined  as: 


h 

r 


ad 


where 


1 


7, 


1  -  2 
—  m 


1  +  1 


5.3  RESULTS 


(5.19) 


(5.20) 


5.3.1  100  FW-2  AND  GIMBALED  NOZZLE 

Since  the  100  FW-2  and  gimbaled  nozzle  walls  are  ablative,  the 
emissivity  may  be  taken  as  unity.  With  the  propellant  conditions  calculated  at 
various  stations  in  the  nozzle.  Equations  (5  18  through  5.20)  were  used  to  fine 
the  radiative  heat  flux. 

The  radiant  heat  fluxes  in  the  chamber  and  throat  regions  of  both  tha 
100  FW-2  and  gimbaled  nozzles  represent  the  heat  flux  from  a  black  particle  cloud 
Therefore,  all  the  assumptions  of  the  preceeding  ^^ection  are  eliminated  except 
Number  (3).  That  is,  the  radiation  from  an  opaque  particle  cloud  Is  independent 
of  geometry,  and  one  must  only  make  a  choice  of  cloud  black  body  emission  tempera¬ 
ture  in  the  calculation  of  the  radiant  heat  exchange. 

Table  3.2  presents  the  radiant  heat  Interchange  at  the  chamber  and 

throat  location  of  the  100  F^'-2  nozzle.  The  chamber  conditions  are:  T  «6350^, 

P  *610  psia,  and  0“O.45. 
c 
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TABLE  '■>.2 


RADIATION  HEAT  FLUXES  IN  THE  100  FW-2  NOZZIL 


LOCATION 

WALL  TEMP  (^R) 

h  ( 3tu/ in^  see  ) 

r 

(Btu/in  sec) 

Chamber 

1  500 

0.0010 

5.4 

Chamber 

JOUO 

0.0015 

3 . 1 

Throat 

1500 

r  00089 

4.3 

Throat 

3000 

o.oou 

4.0 

Figure  5-2  presents  the  radiation  heat  rlux  that  exists  at  various 
positions  in  the  gimbaled  nozzle.  Calculation  of  these  heat  fluxes  are  made  by 
utiliz;.ng  the  procedure  of  subsection  5.2.  The  figure  sht'ws  the  effect  of  wall 
temperpture  on  radiation,  which  is  expected  sinc  e  the  d  if  f erenct^  in  T*  is  being 
compareo.  The  dotted  lines  are  regions  of  uncei  udinty ;  that  is,  regions  that  laay 
not  be  strictly  classified  as  either  the  chamber,  throat,  or  divergent  section  of 
the  nozzle  (see  As sump t ton  (5)) . 

The  i.aipor  tance  of  the  radiative  h^^at  flux  in  the  inlet  section  of  the 
ginba lied  nozzle  is  shown  in  Figure  5-3.  The  parameter  h,./h|-  was  obtained  by 
calculating  h  as  explained  in  subsection  5.2  and  using  the  heat  transfer  coef¬ 
ficient  plotted  in  Ftgure7-d.  The  twe  heat  transfer  coefficients  are  then  added 
to  obtain  (h^  +  h  ).  Note  that  the  sum  (h^  +  h  )  is  not  defined  as  the  total  heat 
flux,  si  nee  the  pr.^ticle  impingeir.t_nt  contribution  is  omit  ted. 

The  Increasing  effect  or  importance  of  radiation  with  burning  time  is 
explained  by  an  investigation  of  h  .  Comparison  of  Figure  5-2,  5-3,  and  5-4 
shows  that  h  decreases  with  burning  time  due  to  the  combined  effect  of  increas¬ 
ing  the  redevelopment  length  of  the  boundary  layer  (pressure  gradient  and  re¬ 
development)  CO  the  throat  and  low  wail  temperature 

Figure  5-4  prerents  the  variation  of  the  radiant  to  total  heat  flux  with 
axial  location  in  the  exit  portion  of  the  gimballed  nozzle.  Wall  temperatures  of 
4000  and  lOCO^^K  are  presented.  The  radiative  and  conv^:uive  heat  transfer  coef¬ 
ficients  were  calculated  using  the  same  procedure  as  Figure  5-4. 

5.3.2  SIMULATOR  NOZZLE 


The  simulator  walls  were  made  of  molybdenum.  With  the  emisslvity  of 
oxidized  molybdenum  found  from  Pvcference  5  8,  and  the  thermochemi cal  properties  of 
the  plus  ^^0^3  slurry  in  nozzle  obtained  from  a  thermodynamic  coini  utor 

program, ^Che  radiant  heat  transfer  was  completed.  The  method  of  calculation  was 
outlined  in  subsection  5.2. 
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FIGURE  S-2  PP£DJ.CTED  RADIATION  HEAT  FLUX  VERSUS  NOZ7-LE  POSITION  FOP  GIMBALLED  NOZZLE 


U.  3 


The  importance  of  particle  cloud  radiation  In  tne  simulator  nozzle  Is 
shown  In  Figure  5*!;*  In  which  it  can  be  seen  that  the  amount  oi  radiation  heat  flux 
is  roughly  doubled  with  the  addition  of  Al^O  particles.  In  the  calculation  of  t\u 
ratlo(h  /h  h  ),  the  effect  of  Che  particles  on  the  convective  transfer  is 
essentliily  omlEted.  That  is,  the  convective  heat  transfer  coefficient  of  Section  ? 
was  Simply  added  to  the  calculated  radiacioii  component. 

Figure  1-6  is  sltniia.*  to  5-5  except  the  total  heat  flux  (h^)  is  found 
from  experiment  and  the  raJistion  Is  based  upon  the  experimental  flame  temperature. 
Comparison  of  Figures  5-5  and  5-6  results  in  a  large  discrepancy  for  the  inlet 
section  of  ^  *  0.  This  deviation  is  attributed  to  mainly  the  difference  in  the 
calculated  and  experimental  convective  heat  flux. 

5.3.3  CONDITIONS  WHEN  UNITY  CLOUD  EMlbSlVlTY  APPLIES 

Black  body  conditions  of  a  particle  cloud  In  a  rocket  nozzle  are  shown  in 
Figure  5-7;  that  is,  determination  of  nozzle  conditions  required  to  produce  a  cloud 
total  emlsslvlty  between  0.9^^  and  1.0  are  plotted  in  Figure  5-7,  and  the  use  of  the 
figure  requires  the  same  assumptions  as  were  made  in  the  construction  of  Figure  5-1. 

In  using  Figure  5-7,  the  path  length  in  an  actual  nozzle  is  obtained  by 
use  of  Assumption  (5)  and  rewritten  here  aa : 

Path  length  *  local  diameter  in  inlet  and  outlet  sections  of  nozzle 
*  1.4  times  the  throat  diameter  in  the  throat  region. 

The  parameter  is  the  product  of  the  ratio  of  particle  to  gas  weight 
and  gas  molecular  weigh*.  Typical  gas  molecular  weights  are  18  for  H-  -  C  and 
21  for  100  FW-2  propellcint.  ^ 

5.4  CONCLUSIONS 

The  total  alumina  cloud  emissivity  is  presented  for  general  cloud  chermc- 
dynaroic  properties . 

Predication  of  the  emisslvit ies  are  approximate  and  should  be  u.sed  onlv 
for  an  oider  of  magnitude  estimate.  However,  tor  a  cloud  with  neglible  multiple 
scattering  as  compared  to  self  absorpt ion^ Figure  5-1  represents  a  good  conservative 
prediction  of  cloud  emission. 

The  radiation  in  actual  rocket  nozzles  has  been  found  empirically.  The 
method  used  in  the  calculation  was  approximate  although  the  importance  of  alumina 
cloud  radiation  as  a  mode  of  heat  flux  was  shown;  and  found  to  be  as  large  as 
50  percent  of  the  total  heat  flux. 

Nozzle  conditions  to  produce  blackbody  effects  In  the  particle  cloud 
have  been  developed  by  applying  existing  experimental  data.  It  was  found,  for  example, 
that  with  thermodynamic  properties  of  T  5200^R,  P  *  350  psia,  alumina  loading 
(  ;)  »  0.5,  and  gas  molecular  weight  (M)  =  20.0;  the  cloud  emlssivity  may  be  taken 
at  unity  for  path  lengths  greater  than  14  inches. 


5-16 


URE  5-5.  PREDICTED  RATIO  OF  RADIATION  TO  RADLMION  PLUS  CONVECTION 


<N 

O 
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STATIC  PRESSURE  (PSIA) 


R07252 


rlGURE  5-7  PATH  LENGTH  REQUIRED  FOR  GIVE!:  PRESSURE  AND 

OF  AN  PARTICLE  CLOUD  TO  PRODUCE  BLACK  BODY  CONDITIONS 


5-19 


C-4 


5.5  REFERENCES 


5.1  Eckert,  E.  R.  G.,  Drake,  R.  M.  Heat  and  Mags  Transfer.  Second  Edition 
McGraw-Hill  (1959) 

r> _ .1- _  rt  v  Ai _ AL _ «.j _ _ 

u,  cu.  ai  ni.umi.ua  rvuaui  aiiu  i-.mi.  c  a  v.».  a.%. 

Publication  No.  U-2627,  Coucract  NOnr  3606(00),  May  1964. 

3.3  Morlzuini,  S.  J,  and  Carpenter,  H.  J.,  ‘‘Thermal  Radiation  from  the 
Exhaust  Plume  of  an  Aluminized  Composite  ^ropellant  Rocket.'*  AIAA 
Aerospace  Sciences  Meeting,  New  York,  January  20,  ^964. 

5.4  J.  T.  Bevans ,  R.  V.  Dunkle,  "Radiant  Interchange  Wi  .i  an  Enclosure" 
Trans  ASME ,  February  1960 . 

5.3  A.  K.  Oppenhelm  "Radiation  Analysis  by  Network  Method",  Trans  ASME, 
Vol  78,  1956. 

5.6  McAdams,  W.  H.  Heat  Transmission.  3rd  Edition,  McGraw-Hill,  (1954). 

5.7  Kreith,  F.,  Principles  of  Heat  Transfer,  International  Textbook  Co., 
(1961). 

5.8  Goldsmith,  A;  Hirschhorn,  H.  J.:  and  Waterman,  T  E.  Thermophysical 
Prop«:irtlea  of  Solid  Materials,  Volume  III,  WADC  Tech.  Report  58-476, 
November  1960 


NOMENCLATIVE 


A  Area  .  f i " 

4  2 

E  Enimisslve  Power  (  7  T  )  ,  Btu/br  ft 

F  Form  Fa^rtor 

h  Radiant  heat  transfer  coefficient,  Btu/tn^  sec^F 

r 

J  Radiosity,  Btu/hr  ft^ 

K  Absorpt i on  coefficient,  1 / f t 

/  Path  length  In  absorbing  medium,  ft 

M  Mach  Number 

M  Gas  molecular  weight 

N  Particle  density,  lb  part/ft^ 

'V 

N  Cloud  absorption  coefficient,  1/ft 

P  Pressure,  psia 

2 

q  Radiant  heat  flux,  Btu/ln  sec 

r 

r  Particle  radius,  cm 

r  Rec  ove  ry  f 4c  t  or 

Intensity  from  absorbing  media 

R  ElacKbody  Intensity  trom  absorbing  media 

Universal  gas  constant,  ft  Lb/lb  -  mole-^R 
°  m 

s  Path  length  (different  from  ft 

T  Temperature, 

U  Velocity,  ft /sec 

W  Average  weight  of  single  particle,  grams 


a  Abs  irptivl ty 

Emlsslvltv 

C  Emlsslvlty  of  single  particle 

representative  an  av«.rage  olre 

Weight  ratio  of  particle  to  gas 

C  Dens Ity,  Ib/ft^ 

^  (r)  Distribution  function  of  particle 
size  ,  ( 1 /cm) 

/3  Polar  angle,  radians 

y  Azimuth  angle,  radians 

r  Transmissivity 

a 

%  Absorption  coefficient,  l/fc 

s 

y  Scattering  coefficient,  l/ft 

y  Specific  heat  ratio 

Subsc  r Ip t 

X  Wave  length  (monochromatic) 

g  Gas 

p  Particle 

t  Total 

1  Wall  condition 

r  Radiative 

m  Mixture 

ad  Free  stream  adabaClc 

s  Stagnation 
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SECTION  6 


SAMPLING  OF  MICRON -SIZED  PARTICLES 


Particles  in  the  size  range  of  0,5  to  6  microns  In  diameter  were  used 
in  the  modeling  and  in  the  hot  firings.  In  both  casen,  there  was  a  need  to 
capture  and  observe  the  particles  of  the  gas  streams. 

Picking  up  the  polystyrene  particles  by  implngeiuent  on  sticky  tape  in 
the  impacting  sampler  allowed  a  check  of  the  sticky  tape  as  an  impingement  de¬ 
tector  on  the  nozzle  model.  The  particle  samples  were  then  used  in  the  develop¬ 
ment  of  a  microscope  technique  for  counting  the  impinged  particles.  It  also 
indicated  whether  or  not  the  particles  were  dispersed  individually  by  the 
atomizing  system. 

Sampling  of  the  rocket  particles  allowed  a  check  on  whether  or  not 
the  particles  were  melted  and  or  how  well  they  were  dispersed  or  agglomerated 
in  the  rocket . 

The  Mav-We  I  chman  partirle  analy?eri  which  was  used,  con.*;  is  ted  of  four 
brass  sampling  nozzles  mounted  in  a  square  ring  at  right  angles  to  each  other 
in  cascade.  Opposite  each  nozzle  was  an  impingement  slide  where  the  particles 
were  collected  from  the  turning  gas  stream.  The  largest  nozzle,  first  stage, 
received  the  rocket  exhaust  and  collected  relatively  large  particles,  and 
passed  on  the  part  not  collected  to  the  next  nozzle.  The  fourth  nozzle  haU 
the  smallest  cross  section  and  collected  the  finest  particles.  Its  collection 
chamber  was  connected  to  a  gas  aspirator.  Only  the  fourth  stage  wes  normally 
used  in  order  to  provide  a  complete  sample  of  all  particle  slzeo  on  one  slide. 
Nozzle  throat  dimensions  for  the  fourth  stage  were  14fnn  x  0.6mcL  in  the  shape 
of  a  rectangular  slot. 
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The  sampling  apparatus  was  developed  especially  for  the  rocket  plume 
using  the  May-We Lehman  apparatus.  It  was  (nounled  in  a  device  (Figs  6-l&6*3)  having 
a  neoprene  sealed  shutter  to  (1)  maintain  steady  vacuum  conditions  in  the  noezle 
system  before  sampling,  (2)  pick  up  a  small  enough  aaciiple  to  show  discrete  parti* 
cles,  and  (3)  omit  sampling  extraneous  dust  before  the  firing.  The  shutter  was 
actuated  by  a  solenoid.  ?>ampllng  time  was  a  function  of  the  duration  of  the 
current  pulse  to  the  solenoid  (approximately  I/IO  second)  and  the  spring  tension 
on  the  shutter  mechanism.  The  sampling  time  was  held  constant  at  approximately 
1/5  second  for  all  samples.  The  apparatus  was  mounted  In  the  center  jf  the 
rocket  plume  approximately  30  feet  from  the  nozzle  exit  . 


The  exhaust  particles  were  collected  on  clear  glass  microscope  slides. 
Different  surface  treatment  of  slides  included  double-sided  Scotch  tape, 
translucent  Scotcli  tape  No.  810,  standard  transparent  Scotch  tape,  and  the 
cleaned  surface  of  tne  glass  slide  itself.  The  clean  glass  slide  produced  the 
most  readable  samples. 


From  an  examination  of  the  particle  samples  several  conclusions  could 

be  drawn: 


(1)  Not  all  the  particles  issuing  fr  m  the  throatleas 
motor  were  molten  when  conditions  were  set  for  a 
chamber  temperature  of  5100®F. 

(2)  There  was  very  little  agglomeration  of  particles 
in  the  throatless  engine. 


(3)  The  engine  with  the  molybdenum  nozzle  for  particle 
impingement  studies  produced  a  wide  variety  of 
particles:  spherical,  crystalline,  agglomerated. 

This  is  expected  since  it  is  known  that  a  con¬ 
siderable  amount  of  alumina  ran  down  the  nozzle 
entrance  and  was  theti  broken  up  at  the  throat  . 

Figure  6-2  shows  some  of  the  particle  displays  produced  by  photo¬ 
micrographs  of  the  sampler  impact  area  from  various  firings. 


6-2 


MOLY  NOZZLE  ^  *  0.3  ALjC^  lOOOX 


MOLY  NOZZLE  0  =  0,3  AL2O3  lOOOX 


MOLV  NOZZLE  0  =  0.2  AL2O3  lOOOX 


THROATLESS  ENGINE  0  =  0.2  AL2O3  lO.OOOX 

FIGURE  6-2.  SAMPLED  ALUMINA  PARTICLES 
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SECTION  7 


CONVECTIVE  HEAT  TRANSFER  ANALYSES 


7 . 1  INTRODUCTION 

Although  the  emphasis  in  this  program  has  been  placed  on  the  enigAa 
of  the  condensed  phase  and  radiation  energy  transport  cnechanisms,  there  is  still 
much  to  be  learned  to  define  accurately  the  convective  process.  In  most  casei,, 
the  condensed  phase  and  radiation  contributions  to  total  heat  transfer  in  «olid 
propellant  nozzles  are  of  the  same  order  as  the  errors  commonly  associated  with 
the  convective  predictions . 

References  7.1  to  7.4  contain  rather  general  treatments  of  the 
important  aspects  of  the  nozzle  heat  transfer  problem.  Some  o2  the  more 
important  conclusions  which  can  be  drawn  from  these  and  other  studies  of  nozzle 
heat  transfer  include: 

(1)  For  the  nozzles  of  interest  in  this  program,  the  convective 
heat  transfer  is  a  consequence  of  turbulent  boundary  layer 
development . 

(2)  Chemical  effects  uue  to  (1)  dissociation  and  recombination, 

(2)  nozzle  wall  ablation  or  pyrolysis,  and  (3)  corrosion 
are  not  negligible. 

(3)  Physical  modification  of  the  boundary  layer  due  to  the 
surface  blowing  is  important. 

(4)  Surface  roughness,  resulting  frorr.  the  condensed  phase 
interaction  with  the  nozzle  wall  or  physical  erosion, 
may  significantly  influence  boundary  layer  development. 
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(5)  It  is  a  good  approximation  to  assume  that  radiative  and 
convective  effects  are  additive  when  the  radiative  flux  is 
significantly  less  than  the  convective  flux. 

(6)  The  direct  addition  of  condensed  phase  impingement  and 
convective  effects  is  regarded  as  a  first  approximation 
which.  In  some  cases,  may  be  very  poor. 

(7)  Non*ideal  combustion  and  non*equilibrlun>  nozzle  expansion 
considerations  can  reasonably  be  avoided  for  the  present 
purposes . 

The  interaction  effects  among  the  various  modes  of  energy  transfer  are 
omitted  from  consideration  in  this  work.  The  consequences  of  this  can  be  illus¬ 
trated  by  the  following  example.  The  convective  heat  transfer  to  an  ablating  or 
pyrolyzing  nozzle  surface  will  lead  to  the  injection  of  gaseous  ablation  products 
into  the  boundary  layer.  The  injection  tends  to  reduce  the  convective  transport 
by  altering  the  physical  boundary  layer  parameters,  surface  shear,  and  boundary 
layer  chemical  composition.  Chemical  reactions,  usually  endothermic,  between  the 
combustion  and  ablation  products  can  further  reduce  the  convective  heat  flux  to 
the  wall.  Since  these  effects  are  proportional  to  the  mass  addition  rate,  which 
depends  on  the  convective  heat  transfer,  equilibrium  will  be  attained,  and  the 
net  convective  transport  will  be  less  than  for  the  case  of  no  injection.  Under 
these  conditions,  if  there  is  also  radiation  heat  transfer,  mass  addition  rates  will 
increase  and  the  convective  heat  transfer  will  be  reduced.  Consequently,  simple 
addition  of  radiation  and  convection  is  conservative »  On  the  other  hand,  the 
Impingement  of  condensed  phases  on  the  nozzle  surface  will  tend  to  increase  the 
boundary  layer  turbulence  level  near  the  wall.  This  will  tend  to  increase  the  con¬ 
vective  heat  transfer  so  that  simple  addition  of  the  effects  will  not  be  conserva¬ 
tive.  Clearly,  when  these  coupling  effects  are  neglected  and  the  simplified 
theoretical  results  are  forced  to  fit  experimental  measurements,  the  effective  heats 
of  ablation  of  the  wall  material  (Q*eff)  will  be  higher  than  expected.  Preliminary 
results  obtained  from  current  work  at  Phllco  suggest  that  film  cooling  and  chemical 
reaction  effects  could  increase  the  apparent  heat  capacity  of  charring  and  non¬ 
charring  ablators  by  factors  of  2  to  5,  or  more,  respe  tively. 

The  primary  objective  in  presenting  these  introductory  remarks  has  been 
to  establish  the  context  within  which  the  convective  heat  transfer  results  apply = 
They  also  serve  to  emphasize  the  need  for  caution  in  the  interpretation  of  the 
results.  Following  a  brief  description  of  the  calculation  method,  the  results  for 
three  specific  nozzle  designs  will  be  presented  and  discussed. 
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7.2  METHOD  OF  CALCULATION 


The  method  of  ':alruiat: Ing  rocket  nozzle  turbulent  boundary  layer 
growth  and  convective  heat  transfer  described  in  Reference  7.4  has  been  used 
exclusively.  The  actual  calculation»  were  accomplished  on  an  IBM  7094  digital 
computer.  Reference  7.4  provides  an  excellent  description  of  both  the  theory 
:r.d  the  computer  program  input  requirements.  Briefly,  the  method  features  the 
simultaneous  solution  of  the  boundary  layer  integral  momentum  and  energy 
equations.  One-seventh  power  profiles  of  velocity  and  temperature  are  used 
to  approximately  determine  the  boundary  layer  shape  parameters.  The  skin- 
friction  coefficient  and  Stanton  nueber  are  evaluated  as  functions  of  boundary 
layer  thickness  from  the  best  available  semlempirlcal  relations.  The  computer 
program  provides  a  number  of  useful  options,  including  choices  of  (I)  one 
dimensional  flow  or  arbitrary  Mach  number  distribution  and  (2)  enthalpy  or 
temperature  driving  potentials.  Although  a  number  of  criticisms  ntay  be  raised 
against  this  method  of  calculation,  It  is  ona  of  the  beat  available,  and  it 
was  found  to  be  both  convenient  and  economical  for  the  present  purposes.  The 
Important  Inputs  to  the  program  and  the  aaaumptlona  made  will  be  diacussed 
briefly  before  presenting  the  results  of  the  conputations . 


All  heat  transfer  calculations  require  the  Initial  daflnltion  of  the 
chamber  conditions,  combustion  product  properties  throughout  the  nozzle,  and 
the  nozzle  geometry.  Subsidiary  thermochemical  computer  programs  are  used  to 
calculate  chamber  conditions,  given  the  original  propclli.nt  composition  and  the 
desired  chamber  pressure.  The  combustion  products  are  then  expanded  Isentrop- 
Ically  for  both  the  frozen  and  equilibrium  flow  asaumptlona.  These  calculations 
provide  Che  chemical  composition,  ratio  of  specific  heats,  and  specific  heat 
data  required  for  Che  convection  analysis.  A  secondary  hand  calculation  la  uaed 
to  relate  pressure  and  mach  number  when  wall  static  pressure  distributions  are 
available  from  experiment.  Normally  the  equilibrium  expansion  assumption  is 
preferred,  in  which  case  the  iser.tropic  exponent. 


evaluated  in  thr?  throat  regioc,  is  taken  as  an  average  value  to  represent  the 
flow.  The  equilibrium  specific  heat, 

C 

P.eq 


,  O  H  t 

■  r^ip 


v:qlculated  throughout  the  expansion  is  not  the  appropriaco  valu*^  to  use  in 
boundary  layer  calculations.  However,  such  values  ca>  be  used  to  determine 
an  approximate  average  va?ue  when  the  nozzle  surface  temperatures  are  known 
')r  assumed.  When  an  average  specific  heat  is  so  detertp'fned,  the  temperature 
driving  potential  is  used  to  c/zlculate  heat  transfer. 


The  specification  of  the  nozzle  wall  temperatures  is  normally  an 
Important  and  complex  problem.  Most  of  the  calculations  were  performed  under 
the  assumption  of  constant  wall  temperature,  primarily  for  convenience.  Actual 
wall  temperatures  vary  with  time,  axial  position,  and  wall  material.  Additionally, 
gas  properties  and  friction  coefficients  were  evaluated  at  the  film  reference 
temperature . 

The  Prandtl  number,  stagnation  viscosity,  and  the  temperature  dependence 
of  the  viscosities  were  determined  from  the  data  given  in  Reference  7.6.  The  free 
stream  stagnation  temperatures  were  taken  equal  to  the  ideal  values  for  the 
measured  or  time  averaged  chamber  pressures.  The  heat  transfer  program  assumes  a 
recovery  factor  of  0.89  in  calculating  adiabatic  wall  temperatures. 

In  all  the  calculations,  the  boundary  layer  is  assumed  to  grow  from  an 
initial  momentum  thickness  slightly  greater  than  zero.  Similarly,  the  initial 
ratio  of  thermal  to  velocity  boundary  layer  thicknesses  is  taken  to  be  slightly 
greater  than  unity.  The  effects  of  these  relatively  arbitrary  specifications  die 
out  rapidly  and  do  not  seriously  compromise  the  results.  It  should  also  be  noted 
that  during  the  initial  development  of  the  boundary  layer,  laminar  flow  is  normally 
expected  whereas  the  calculations  regard  the  entire  boundary  layer  development  as 
turbulent.  Consequently,  the  first  heat  transfer  coefficient  calculated  is  not 
accurate.  However,  the  first  point  is  expected  to  be  directionally  correct,  as  it 
should  be  higher  than  the  succeeding  point  or  points  due  to  the  importance  of  the 
boundary  layer  development  length  effect. 

Consistent  with  the  recommendations  of  Reference  7.4,  a  boundary  layer 
interaction  exponent  or  0.1  was  chosen  for  all  convection  calculations.  This 
quantity  provides  for  the  effective  increase  Iri  the  heat  transfer  coefficient  v;hen 
the  boundary  layer  energy  thickness  exceeds  the  momentum  thickness. 

A  few  additional  assumptions  concerning  the  nozzle  geometry  and  boundary 
layer  starting  points  are  required.  These  are  briefly  discussed  along  with  the 
results  of  the  calculations  for  eacu  particular  nozzle  examined. 
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Calculations  were  performed  for  three  different  nozzles.  The  results  are 
presented  separately.  Note  that,  except  '.’here  noted  otherwise,  the  symbols  which 
appear  In  the  figures  which  follow  represent  calculated  points,  not  experimental 
measureaiencs . 


7.3.1  MOLYBDENUM  CALORIMETRIC  NOZZLE 

The  molybdenum  nozsle  geomet.'y  and  thermocouple  locations  are  shown  in 
Figure  3-3  and  in  Table  7.1.  Tlie  convective  heat  transfer  for  Run  lA,  frozen  and 
equilibrium  expansion,  and  Run  45A  are  shown  in  Figure  7-1  for  comparison.  The 
boundary  layer  was  assuaued  to  start  on  the  conical  (55  degrees)  inlet  section  si  a 
radius  equal  to  the  chamber  (grain  port)  radius.  It  was  necessary  to  fit  a  smooth 
curve  to  the  available  thermocouple  data  to  provide  an  adequate  wall  temperature 
distribution.  Heat  transfer  coefficients  were  calculated  at  a  single  tine  near  the 


7-4 


TABLE  7.1 


MOLYBDENUM  SIMULATOR  NOZZLE  GEOMETRY 


Axial  Diatance 

Z,  Inchaa 

Radius.  Inch 

Area  Ratio 

Thertaocouple 

Nuadtex 

0 

3.0 

23.0 

- 

0.418 

2.405 

14.8 

k 

0.618 

2.119 

11.5 

2 

0.808 

1.847 

8.7 

3 

0.998 

1.575 

6.35 

4 

1.198 

1.289 

4.25 

5 

1.290 

1.159 

3.44 

- 

1.430 

0.991 

2.4 

6 

1.689 

0.793 

1.61 

7  . 

1.990 

0.667 

1.14 

8 

2.314 

0.625 

1.0 

9 

2.638 

0.667 

1.14 

10 

2.680 

0.680 

1.18 

- 

2.998 

0.773 

1.53 

11 

3.318 

0.867 

1.92 

12 

3.658 

0.966 

2.39 

13 

3.98b 

1.063 

2.9 

14 

4.328 

1.162 

3.47 

15 
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FIGURE  7-1.  PREDICTED  AND  EXPERIMENTAL  Cor  ECTIVE  HEAT  TRANSFER  COEFFICIENT 
ViaiSUS  NOZZLE  AXIAL  POSITION  F'lR  MOLYBDENIM  CALORIMETRIC  NOZZLE 
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end  of  the  run  when  chamber  and  wall  temperature  conditions  were  reasonably  stable 
U.3  seconds).  It  should  be  noted  that  these  results  must  be  Interoreted  as 
corresponding  to  the  theoretical  flame  temperature  with  a  recovery  factor  of  0.89, 
For  the  purpose  of  calculating  viscosity,  the  composition  of  the  combustion  gases 

»2.  Also  shown  In  Figure  7-1  are  con- 

theo^*Mr-l^‘‘"if?'^f  calculated  directly  from  test  data  of  Run  lA.  As  shown,  the 
hi  h  ^  coefficients  for  froaen  expansion  In  the  throat  region  ^re  50  percent 

experimentally.  In  the  entrance  and  high  area  ratio  exit 
locations,  theoretical  and  experimental  coefficients  are  more  nearly  equal. 


7.3.2  lOOFW-2  NOZZLE  INLET 

Convection  coefficients  calculated  for  the  soft  ablator  portion  of  the 
lOOFW-2  nozzle  Inlet  ate  presented  in  Figures  7-2  through  7-6.  Each  of  these 
figures  represents  a  particular  time  during  the  firing,  and,  consequently.  It  was 
necessary  to  determine  the  appropriate  ablator  geometry  at  each  time.  The  Initial 
and  final  contoura  were  known  from  the  Inlet  design  and  post  test  measurements 
(Reference  7.5).  At  5  seconds.  It  was  assumed  that  the  exposed  portion  of  the 
original  contour  was  unchanged.  At  the  end  of  firing,  90  seconds,  the  poet  test 
measured  contour  was  assumed.  The  19.5,  50,  and  75  second  contours  were  con¬ 
structed  by  assuming  that  the  Insulation  surface  regression  was  a  linear  function 
of  the  time  It  was  exposed.  In  all  cases,  the  heat  transfer  calculations  were 
stopped  at  the  point  where  the  soft  Insulation  met  the  graphite  substrate,  and 
they  were  started  at  the  point  where  the  propellant  grain  met  the  soft  Insulation. 
The  assumed  geometries  are  shown  In  Table  7.2.  The  two  curves  shown  In  each 
of  the  figures  are  a  consequence  of  assuming  constant  wall  temperatures  of 
1500  and  3000°R.  Considerably  less  difference  would  be  observed  If  the  coro- 
buslon  gas  properties  and  friction  coefficients  had  been  evaluated  at  free 
stream  conditions.  The  combustion  gas  viscosity  was  determined  by  assuming 
proportional  contributions  by  the  major  gas  components  using  the  data  of 
Reference  7.6.  The  values  of  the  convection  coefficients  near  Z  «  0  can  certainly 
not  be  regarded  as  accurate  due  tc  the  uncertainty  of  the  grain  hnrning  and  flow 
field  details  in  this  region.  In  general, Figures  7-2  through  7-6  show  the  ex¬ 
pected  trend  of  decreasing  heat  transfer  coefficient  with  Increasing  distance 
from  the  origin  of  the  boundary  layer.  This  trend  would,  of  course.  Indicate 
that  the  linear  Insulation  recession  rate  assumption  was  not  accurate,  unless 
one  of  the  other  heat  transfer  modes  has  a  reverse  tendency. 


7.3.3  GIMBALL  NOZZLE 


A  full  scale  glmballed  nozzle  was  scaled  down  (5.75  to  1),  and  the  wall 
pressure  distributions  were  measured  in  both  the  straight  and  canted  (5  degree) 
positions.  The  convective  heat  transter  has  been  calculated,  using  the  modeling 
pressure  data,  tor  the  full  scale  nozzle.  The  lOOFW-2  propellant  and  chamber 
pressure  were  assumed  for  these  calculations.  Because  of  the  severe  Interruption 
In  the  nozzle  contour  upstream  of  the  throat,  it  was  necessary  to  calculate  the 
heat  trap  'r  In  the  inlet  and  throat-exit  sections  separately. 


COHVECTIVE  HE^iT  IDANSFER  COEFFICIENT  ( 


FIGURE  7-2.  PREDICTED  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT’  VERSUS  AXIAL 
DISTANCE  F(»  lOOFW-2  AFT  CLOSURE  AT  5  SECONDS 
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TIVE  HRM  TRANSFER  COEFFICIEHT  (Btu/IN  SEC  F  X  10  ) 


AXIAL  DISTANCE  (INCHES) 


R07487 


ncuu  7-3.  PtKDiCTn)  COHVECnvi  HKAT  lUAMSFtt  COSmCinTT  VERSUS  AXIAL 
DISTAMCB  FOR  lOOFW- 2  AFT  CLOSURE  AT  19.5  SECONDS 
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FIGURE  7-5.  PREDICTED  CONVECTIVE  HEAT  TRANSFER  COEFFICIEWT  VERSUS  AXIAL 
DISTANCE  FOR  lOOFW-2  AFT  CLOSURE  AT  75  SECONDS 


CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  (Btu,IN^  SEC°F  X  lo'^) 


AXIAL  DISTANCE  (INCHES) 


RO7479 


FIGURE  7-6.  PREDICTED  CONVECTIVE  HEAT  TP.ANSFER  COEFFICIENT  VERSUS  AXIAL 
DISTANCE  FOR  lOOFW-2  AFT  CLOSURE  AT  90  SECONDS 
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TABLE  7.2 


lOOFV-2  NOZZLE  INLET  GEOMETRY: 


5  Seconda  After  Ignition 


.  Position 

Radius 

Axial  Position 

T  nrhi»a 

Inches 

Z.  Inches 

0.0 

20.87 

4.0 

l.O 

19.84 

5.0 

z .  u 

1ft  AO 

X  7  .  V 

11.5 

3.0 

18.35 

19.5  Seconds 

After 

Ij^nlclon 

0.0 

26.77 

4.5 

0.5 

25.92 

5.5 

1.5 

24.62 

6.0 

2.5 

23.35 

15.0 

3.5 

22.22 

30  Seconds 

After 

Ijfjni  tlon 

0.0 

36.22 

7.0 

1.0 

34.97 

8.0 

2.0 

33.67 

9.0 

3.0 

32.02 

IC  0 

4,0 

30.37 

10.7 

5.0 

28.37 

20.5 

6.0 

27.37 

75  Seconds 

After 

Ignition 

0.0 

43.7? 

11.0 

1.0 

42.17 

12.0 

2.0 

41.17 

13.0 

3.0 

41.15 

14.0 

4.0 

40.87 

15.0 

C  /\ 

^  *  u 

40.37 

16, C 

6.0 

39.77 

17.0 

7.0 

39,07 

18.0 

8.0 

38.17 

19.0 

9.0 

37.47 

20.0 

10.0 

36.47 

30.5 

90  Seconds 

After 

Ignition 

0 

50.0 

15.0 

1.0 

49.1 

16.0 

2.0 

48.2 

17.0 

3.0 

47.3 

18.0 

4.0 

46.4 

19.0 

5.0 

45.5 

20.0 

6.0 

44.6 

21.0 

7.0 

43.8 

22.0 

8.0 

42,3 

23.0 

9.0 

41,3 

24.0 

10.0 

41.2 

25.0 

11.0 

41.0 

26.0 

12.0 

40.47 

27.0 

13.0 

39.87 

37.5 

14.0 

39.20 

-  throat  " 


-  throat  - 


-  threat  - 


~  throat 


Radius 

Inches 

17.68 

17.02 

14.20 


21.22 

20.52 

20.25 

14.27 


26.17 

24.87 

23.67 

22.67 

22.17 
14.41 


35.37 

34.17 

32.57 
30. Y7 
2«>,&7 
28.  r, 
2>.  7 

26.17 

24.9/ 

23.8-’ 

14.53 


38.47 
37.60 
36 . 68 
35.59 
34.39 

32.87 
31.49 
30.09 

28.87 
27.80 

25.64 

25.57 

24.57 

14.65 


throat 


The  results  for  the  inlet  section  are  shown  in  Figures  7-7  through  7-9 
for  grain  burn  times  ot  0,  23,  and  58  percent  respectively.  The  geometry  is  given 
In  Table  7-3.  In  Figure  7-7,  the  wall  temperature  (graphite  phenolic)  was  taken  to 
be  constant  at  lOOO^R  to  approximately  repre.«‘ent  the  starting  conditions.  The  two 
curves  in  Figure  7-8  represent  assumed  wall  temperatures  ol  4000  and  SOOO^R. 
Similarly  the  curves  in  Figure  7-9  represent  5000  and  6000^  surface  temperatures. 
Although  these  temperatures  are  arbitrary,  they  are  thought  to  cover  the  range  of 
wall  conditions  which  might  realistically  be  expected.  The  boundary  layer  calcu¬ 
lations  were  started  at  the  propellant  grain  -  wall  coriCour  junction.  The  calcu¬ 
lations  were  all  stopped  at  a  radius  of  7.60  inches  which  corresponds  to  the  first 
throat  rndius  in  the  modeling  nozzle.  The  pressure  m"a8urement8  were  not  used  in 
these  calculations.  Figure  7-10  compares  the  Mach  number  distribution  correspond¬ 
ing  to  the  one  dimensional  flow  assumption  with  the  Mach  numbers  determined  from 
the  model  inlet  pressure  taps.  The  differences  are  not  regarded  as  especially 
significant  compared  to  the  errors  which  could  arise  from  fairing  a  curve  through 
the  three  data  points. 


The  calculated  heat  transfer  coefficients  for  the  throat-exit  section  of 
the  full  scale,  gimballed  nozzle  are  shown  in  Figures  7-11  through  7-13.  The 
geometry  and  Mach  .;umber  distributions  are  given  in  Table  7.4.  Figure  7-11  compares 
the  heat  transfer  for  the  one  dimensional  flow  assumption  with  that  calculated  using 
the  modeling  wall  pressure  distribution  for  a  5000^R  isothermal  wall.  Also  shown 
are  the  results  for  a  non-lsothermal  wall.  Here  the  wal^  temperature  was  taken  at 
6000^R  in  the  throat  section,  and  then  varied  uniformly  to  the  local  free  stream 
static  temperature  at  the  lart  v;all  point  In  the  exit  cone.  This  temperature  distri¬ 
bution  is  quite  arbitrary,  but  it  averages  about  5000^R  while  decreasing  from  OOOO^R 
to  about  4000^.  Figure  7-12  compares  the  heat  transfer  coefficients  determined  for 
isothermal  walls  at  ICOO  and  4000^R  w:  th  the  nozzle  In  the  straight  position  and 
using  the  modeling  pressure  data  •'o  determine  Mach  niimber .  The  5000^  wall  temper¬ 
ature  values  shown  in  Figure  7-11  were  not  ceplo^ted  because  they  are  very  close  to 
the  4000^R  results.  The  one-dlmensiona I  calculations  are  replotted  for  comparison. 
Figure  7-13  presents  the  data  for  the  nozzle  in  the  canted  position.  Modeling 
pressure  data  was  used  to  determine  the  Harh  number  distribution  for  each  of  the 
three  curves.  The  results  are  also  based  on  the  assumption  of  a  500C’^R  isothermal 
Will  condition.  The  canted  nezzie  wall  coordinates  are  referred  to  the  axial 
centerline  of  the  exit  cone.  The  boundary  layer  v?as  started  5.40  inches  upstream 
of  the  15.0  inch  diameter  throat.  This  is  apparently  a  reasonable  assumption, 
beciiuae  it  Is  extremely  doubeful  chat  the  boundary  layer  formed  in  t  ie  nozzle 
entrance  section  could  persist  across  the  major  discontinuity  in  the  contour 
uetween  the  two  throats.  The  initial  Mach  aujber  of  0.80  was  determined  from 
the  model  pre.*9sure  measur ament  ai  sta:iori  4,  and  wac  the  same  for  both  the  canted 
and  straight  configurations. 


7.4  CONCLUSIONS 

In  compf.ring  the  predicted  heac  transfer  coefficients  for  the  molybde¬ 
num  nozzle  with  the  values  obtained  from  Run  lA,  as  shown  in  Figure  7-1,  large 
differences  ere  observed  in  the  throat  region.  Referring  to  Figure  1  of  refer¬ 
ence  7.7,  it  Is  evident  that  transitional  boundary  layer  flow  could  occur  in  the 
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FIGURE  7-8.  PREDICTED  COW^ECTIVE  HEAT  TRANSFER  COEFFICIENT  VERSUS  AXIAL 
DISTANCE  FOR  GIMBALLED  NOZZLE  INLET  AT  23%  BURN  TIME 


CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  SEf 


AXIAL  DISTANCE  (INCHES)  R07C78 


FIGURE  7-9.  PREDICTED  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  VERSUS  AXIAL 
DISTANCE  FOR  GIMBALLED  NOZZ  j:  INLET  AT  58X  BURN  TIME 
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TABLE  7.3 


GIMBAL  NOZZLE  INLET  GEOMETRY 


07L  Grain 

23%  Grai 

T| 

58%  Crain 

Axial 

Axial 

Axial 

Position  Radius 

Posits  on 

Radius 

Position 

Radius 

Z.  Inches  Inches 

Z,  Inches 

Inches 

Z,  Inches 

Inches 

0.0 

9.78 

0.0 

14.95 

0.0 

22.4 

0.5 

9.05 

0.8 

13.80 

0.95  - (3>— 

21  .05 

1.0 

8.45 

1 . 6  — 

12.65 

2.05 

19.45 

1.5 

8.10 

2.6 

11.25 

3.2 

17.85 

2.0 

7.86 

3.6 

9.78 

4.2 

16  45 

2.5 

7.73 

4.1 

9.05 

5.2 

14.95 

3.0 

7.66 

4.6 

8.45 

6.0 

13.80 

3.5 

~  7.60 

5.1 

8. 1C 

6 . 8  - 0 — 

12  60 

4.6 

-Throat-  7.50 

5.0 

7.86 

7.8 

11  25 

6.1 

7.73 

8.8 

y  .78 

6.6 

7.66 

‘».3 

9  05 

7  1  — (3) - 

7.60 

9.8 

1.45 

3.2  -Throat- 

7.50 

10.3 

8.10 

10  8 

7  86 

11.3 

7.7.1 

11.8 

7.66 

12.3  - 0 - 

7.60 

13.4  -Throat- 

•7.50 

NOTE:  Indicates  locatioi 

1  of  pressure  tnp  number 

auown  in  circle. 
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NOZZU  lAblOS  (raCH) 


2^.0 


CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  (Btu/IN  SEC°F  X  10‘") 


FIGURE  7-11.  PREDICTED  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT 

VERSUS  AXIAI.  DISTANCE  FOR  GIMBALLFJ)  NOZZLE  THROAT 
EXIT  STRAIGHT  CONFIGURATION 
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COKVECTIVE  HEAT  TRANSFER  COEFFICIENT  (Btu/ln'  SFC 


AXUl  DISTANCE  (INCHES)  R07481 

FIGURE  7-13.  PREDICTED  CONVECTIVE  BEAT  TRANSFER  COEFFICIENTS  VERSUS  AXIAL 
DISTANCE  FOR  GIMBALLED  NOZZLE  THROAT  -  EXIT  CANTED  (5°) 
CONFIGURATION 


7-22 


I 


I 

[ 

[ 


j|^|  doeboooo^J^J^M^MMMM*^ 


t 

s 

9 

«4 

44 

i 

1 


t  •!  e«p<^lQ04r>^ASoo^'^01^4cQOcb«r^ 

2  ^  ®  ®  «  ao  r^WV  r-‘  r«:  r-'  r.*  ®  »  SC 

(^  (^  (^  (^(^(^ 

S  1 


II 


o 


®«/s»M9>^'«<*\<«4rt>i9^^0^<£P4>rs,00 


^  M  ^  r<4  rt 


I 


I 

r 

r 

I 

1 

I 

1 


1^ 

U4 

»-] 


O 

H 

M 

X 

UJ 

I 

H 

<• 

O 


z 

O 

H 

CQ 

K-l 

f5 

cn 


Di, 

U 


hJ 

csl 

M 


3: 

o 


•:3 

to 


o 


o 

J 

»« 

w 

5 


a 

3 


1^ 


Ann  ®r^r^r^r«r>.r^r^t^<CQ09'9'^(S|<^«^ 


^  0  0  0  000 
O  I 


•^  •-•  «^  fM  ^ 


o<j»«r>2j^a»,4»M^ao^oor^ 


^  »W  »/'  00  ®  ^ 

O'  Ov  «>J 


o  o  o 


3 

o 

u 

J 

« 

•a 

3 

'i 

a 

3 


0  0  0  000 


rM<'^«n4nt/%(npOOOO<rit/'>^>/'r^<'^r^Q 


t» 

i 

a 

a. 


1 

1 

I 


c 


0 


0xr'0<^<~<04>^^00<c 


OoorMOOiTiO^OOMi-'pfMr'OOrstflBr-- 

0''i4^0®<C>0</'ir^,0^0'OrMOoOf'4^r^® 


<«>flO®cOr‘-r^r~f>--«-»~*^t'.  X>:A9'{T^i-t<M-Ox' 


00)0 


o 


u 

m 

«  ca 


3  5 

iH  • 


I 


r-*O^NO< 


;8S8 


M 


0 


throfit  for  diameter  Rtynolds  nuriibers  ranging;  from  about  340,000  to  740,000.  The 
Reynoldii  throat  number  for  Run  lA  was  approximately  520,000,  and  Is  well  within 
the  transitional  range.  This  suggests  that  the  thi^oat  heat  transfer  coefficient 
could  be  about  0.0013,  compared  to  the  measured  value  of  0.0018  Btu/ in^sec  F. 

This  difference  Is  not  unreasonable,  considering  the  potential  influences  of  heat 
transfer,  combustion  phenomina  and  nozzle  inlet  geometry  on  boundary  layer  transi¬ 
tion.  However  j  i  may  also  be  observed  that  the  best  turbulent  boundary  lay?r 
correlation  shown  in  Figure  1  of  reference  7.7  is  in  very  close  agreement  with  the 
experimental  value  for  Run  lA.  Consequently,  the  possibility  that  the  boundary 
layer  in  the  throat  is  not  fully  turbulent  is  neither  positively  accepted  or  re¬ 
jected.  There  are,  of  course,  other  experimental  variables,  such  as  the  distri¬ 
bution  of  the  Injected  water  which  could  account  for  the  low  heat  fluxes  observed. 

For  all  the  other  nozzles  examined  in  this  study,  the  throat  Reyiiolds 
numbers  are  outside  the  range  for  transitional  flow  suggested  in  reference  7.7. 
Further,  the  comparison  of  theory  with  experiment  for  liquid  engine  firings  shown 
in  reference  7,4  indicates  that  the  Elliott,  Bartz,  Silver  method  may  tend  to 
predict  low  rather  than  high  heat  transfer  coef ficier.t s  for  the  nozzles  considered 
here.  Again,  it  is  emphasized  that  there  is  much  to  be  learned  about  the  funda¬ 
mental  nature  of  the  connective  heat  transfer  processes  in  rocket  nozzles. 

The  predicted  heat  transfer  rates  for  the  100FU'-2  nozzle  :'.nlet  show 
several  important  trends  which,  although  they  are  not  unexpected,  should  be 
emphasized.  First,  the  relative  insensitivity  of  the  heat  t’^ansfer  rate  to  the 
choice  of  surface  temperature  is  evident  at  each  of  the  firing  rimes  (Figures  7-2 
to  7»6).  The  results  also  graphically  display  the  competition  of  the  development 
length  and  acceleration  effects  on  the  heat  transfer.  Tn  the  first  few  inches 
(small  Z),  the  boundary  layer  thickens,  and  the  convective  coefficient  decreases, 
roughly  as  tht  distance  to  the  one-fifth  power,  as  suggested  by  flat  plate  theory. 
In  the  region  where  the  main  scream  flow  is  accelerating,  the  coefficients  rise 
again^  and  would  ultimately  reach  a  maximum  at  the  nozzle  throat.  It  is  evident 
that  the  heat  transfer  coefficients  at  any  point  along  the  inlet  surface  will  be 
dependent  on  the  distance  over  which  the  boundary  layer  has  developed  before  reach¬ 
ing  the  point  in  question.  That  is,  the  local  convective  heat  transfer  coef¬ 
ficient  will  exhibit  a  sigr4lf leant  inverse  dependence  on  the  time  after  ignition. 
The  calculated  heat  transfer  coefficients  were  subsequently  used,  along  with  the 
radiation  and  condensed  phase  heat  flux  contributions,  to  estimate  the  effective 
Q  for  the  aft-closure  liner. 

The  heat  transfer  results  for  the  gimballed  nozzle  inlet  appear  to  be 
quite  insensitive  to  wall  temperature  at  high  surface  temperatures.  Figure  7-7 
indicates  that  the  combination  of  low  wall  temperature  and  short  boundary  layer 
development  length  can  lead  to  relatively  higher  coefficients  during  the  initial 
phase  of  combustion.  Since  this  is  the  time  when  the  ablative  mass  addition  occurs 
at  the  greatest  rate,  tht  net  heat  transfer  rates  could  exhibit  the  opposite 
tendency  In  the  actual  motor.  Because  the  calculated  coefficients  are  rather  high, 
the  wall  erosion  rate  Is  likely  to  be  higher  than  might  otherwise  be  expected. 
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The  results  tor  the  throat-cxlt  region  of  the  glmalled  nozzle  are  rather 
interesting.  Again,  a  major  dependence  ot  the  results  on  the  choice  of  the  surface 
temperature  is  not  observed.  The  results  obtained  under  the  assumption  of  one- 
dimenslonal  flow  reflect  the  effect  of  the  Initial  boundary  layer  development,  which 
would  not  appear  if  the  "pipe  flow"  or  "simplified  Bartz"  type  analyses  had  been 
used  to  calculate  the  convective  coef 1 Iclents .  The  use  of  the  Mach  number  distri¬ 
bution,  as  determined  from  the  modeling  pressure  data,  has  the  interesting  effect 
of  completely  eliminating  the  peak  In  the  convective  coefficient  at  Che  nozzle 
throat.  This  suggests  that  the  nozzle  surface  Just  upstream  of  the  geometric 
throat  will  be  subjected  to  an  unusually  severe  environmental  interaction  leading 
to  high  erosion  rates.  Figure  7-13  suggests  that  the  effect  of  nozzle  canting  Is 
to  introduce  about  a  20  percent  'ariatlon  In  the  convective  coefficient  around  the 
nozzle  circumference.  It  Is  also  noted  that  the  effects  are  confined  to  the 
vicinity  of  the  throat  with  the  greatest  deviations  from  the  uncanted  results 
occurring  downstream  of  the  throat. 
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SECTION  8 


PREDICTION  OK  EROSION  RATES  FOR  1 00  FW2  AFT  CLOSURE 


At  the  completion  of  Contract  AF  0^4  (6 1 1 )  -  7C<47  ,  30  October  1962, 
a  preliminary  attempt  was  u’.ade  to  predict  Insulation  erosion  rates  for  the 
100  FW2  aft  closure  using  calculated  convective  and  radiative  heat  transfer  rates 
and  estimated  rates  of  particle  impingement.  In  view  of  the  results  of  cold 
flow  modeling  studies  completed  under  che  present  contract,  a  refinement  of  this 
prediction  was  investigated,  taking  into  account  revised  estimates  of  convective 
and  radiative  heat  transfer  as  well. 


8.1  SUMMATION  OF  HEAT  FLUXES 

The  total  heat  transfer  to  the  aft  closure  surface  was  considered  to 
consi.st  of  three  components: 

(1)  Convective: 

q  =  h  (T  -  T  ) 
cow 

(?)  Radiative: 

Q  * 

(3)  Impingement: 

9^  * 


where 


A  A 
eo  (T  -  T  ) 

p  w 


aih  C  (T  T  ) 
p  v* 


(8.1) 


(8.2) 


(8.3) 


q  >:  convective  heat  vlux  per  unit  area 
c 


convective  heat  tran.sfer  coefficient 
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impingement  heat 
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m 

particle  energy 

ni 

P 

mass  rate  of  par 

C 

average  particle 

r .  I  tAl I  r 

lux  per  unit  nriM 
vl  ty 
uMi  s  tn  lit 

ture 

I  lux  pel  unit  nren 
a f  c ninnioda  t  ion  t  a c  t o r 
ticle  flux  per  unit  area 
heat  capicity 


All  cfllculatlors  ol  instantaneous  heat  i  luxes  were  eonver^cd,  loi  purpo«;es 
of  predicting  erosion  la^es,  to  time-average  heat  fluxes  h,  the  three  mechanisms 
at  any  given  station  along  the  aft  closure  The  average  tlux  at  a  given  station 
was  defined  as 


(q) 


avg 


/  q  dt 

t 

s 


(H  4) 


where 

t i me-average 

t 

e 

total  burn  t 

t^ 

time  at  whic 

back  of  the 

q 

m 

ins tancaneou 

t 

m 

t  ime 

heat  flux 
Ime 

h  a  g^ven  station  is  exposed  burning 
propellant  grain 

s  heat  fiux  (either  qc  .  qr ,  or  qp) 


The  aver-ige  convective  heat  fluxes  were  derived  by  computer  integ.diii-n 
of  Instantaneous  heat  fluxes,  based  on  profiles  of  convective  nec.t  transfer 
coefficient,  **h^,**  calculated  for  several  different  times  during  the  firing.  The 
average  radiative  heat  flux  was  the  same  at  all  points  and  was  calculated  directly 
from  Equation  (8.2).  Inasmuch  as  integrated  values  of 


s 


were  available  directly  from  the  results  of  the  cold  flow  modeling  studies,  only 
hand  calculations  likewise  were  required  to  determine  time-average  values  of  "q 
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The  1 I  low  i  assunipt  i  ons  w^re  tivui(*  i  n  a  1  1  ea  ’  ru  I  a  I  1  i  tis  : 


(1) 

(2) 


,  o 

T  »  T  -  hj/iO  K 

o  P 

^  *  I  . 


( 3 )  a  -  ]  A) 


Calculations  were  perlc^rnied  loi  twt>  possible  values  ol  wall  I  enipera  t  u  re : 

(1)  T  ■  15()0"r 

w 

(2)  T  =  U)()1)‘’r 

w 

The  time- average  heat  t  luxes  computed  for  the  throe  mechanisms  and  the 
two  wall  temperatures  .nrc  shown  in  Figure  8-1  as  functions  of  the  tangential 
distance  along  the  aft  closure,  where  "zerc»'‘  Is  defined  as  the  point  where  the  aft 
closure  meets  the  100- inch  diameter  case,  and  where  the  throat  then  is  located 
at  distar.cc  53.4  Inches.  It  is  seen  that  convective  and  radiative  heat  fluxes 
are  of  the  same  order  of  magnitude,  whereas  1 npingement  heat  flux  calculated  from 
the  nev  mode  I  tng  d.it.4 ,  genera  1  ly  is  so  low  as  to  be  virtually  insignificant. 


The  time-average  heat  fluxes  shown  in  Figure  8-1  can  be  converted  into 
time-overage  erosion  rat^^s  as  follows: 


1  ■  ^total 

'  ‘c  i  Q* 

where 

A  “  time-aveiage  erosion  rate 

«  density  of  virgin  insulation 

■  density  of  porous  char 

■  **effective  heat  of  ablation**  of  insulation 

and  where 


(8.3) 


^Lotal«q+q+q  (8.6) 

c  r  p 

tTime- average  values  are  assumed) 

This  calculation  presupposes  that  the  insulation  ablates  to  a  weak  char  which  Is 
then  mechanically  removed,  so  that  only  the  density  difference  between  the  virgin 
material  and  the  char  is  countea  as  undergoing  the  ablation  process.  Densities 
typical  of  the  100  FW2  Insulation  were  to  be: 

-  80.5  Ib/it^ 

r  -  24.2  Ib/ft^ 

c 
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FIGURE  8-1 .CALCULATEP  AVERAGE  KEAi  FLUXES  FOR  lOO-FW-2  AFT  CLOSURE 


8.2  ESTABLISHING  FOR  THE  WALL  MATERIAL 

Th<'  primary  assumption  involved  in  tht  conversion  of  hef'jt  ilex  rates 
Into  erosion  rates  Is  the  choice  of  a  value  for  the  ’’effective  heat  of  ablation/* 
Q^-  This  quantity,  which  snould  not  be  expected  to  be  the  same  as  any  '’heat  of 
ablation”  measured  in  the  laboratory  under  static  conditions,  crtjst  include  the 
otherwl se- unaccount ed" for  elfects  of  blowing,  mechanical  shearing,  and  the  like. 

Originally,  attempts  were  made  to  predict  erusion  rates  catching  those 
actually  measured,  with  assigned  a  single,  constant  value  throughout.  These 

attempts  proved  unsuccessful;  hence,  a  re  la  C  i<'>nsh  ip  dependent  on  heat  flux 
was  Inves t iga ted .  Inasmuch  as  both  the  radiativ<^  and  tht  impingement  heat  f  lucres 
were  essentially  constant  thr<>ughout  the  region  of  interest  (true,  time-average 
erosion  rates  actually  are  mijaaurable  only  between  distance  points  14.4  Inches  and 
37.4  Inches  in  Figure  8-1),  Q*  was  assumed  to  be  an  arbitrary  power  function  of 
average  ”qc  convec  t  ivr  he^  t  f  lux  .  The  wall  temperature  was  taken  to  be 

1500OR,  since  this  case  would  maximize  thc^*  relative  importance  of  the  convective 
contribution.  The  empirical  relationship  for  Q*  best  fitting  the  measured  time- 
average  erosion  rates  then  was: 

Q  .  (3.22  X  10®)  (q^)  ^  ^  (8.7) 

where 

★ 

Q  *  effective  heat  of  ablation,  Btu/ i b 

2 

■  convective  heat  flux,  Btu/in.  -  sec 

Total  time-average  erosion  rates  predicted  by  Equations  (8.5),  (8.6), 
and  (8.7)  are  compared  it;  Figure  8-2  to  the  measured  average  rates  calculated  as 


expt i . 

where 

/\ 

expt 1 . 
A 

expti . 


^  1 

expt  1 

t  -  t 
e  s 


(8.8) 


time -average  eros ion  ra te 


erosion  depth  measured  at  a  given  station  along  the  aft 
c losure 
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8.3  INTERPRETATION  OF  THE  ANALYSIS 


The  ‘‘predicted*’  curve  based  upon  the  empirical  Q  relationship  Is  seen 
to  fit  the  experimental  curve  very  well,  except  possibly  in  the  distance  region 
30  Inches  to  36  inches .  Actual  Q  values  required  to  fit  the  measured  data 
ranged  from  56,000  Btu/lb  at  position  14.4  inches  down  to  12,500  Btu; lb  at  ^  ^ 
position  37.4  Inches.  The  latter  value  is  in  the  range  of  previous  estimates'^' ^ 
of  effective  heat  of  ablation  under  rocket  firing  conditions,  while  the,  former 
value  definitely  appears  too  high.  Perhaps  the  measured  erosion  rates  ^^^ptl 
are  subject  to  large  errors  in  the  region  of  the  14  4-inch  point  owirg  to  .small 
values  oi  both  *"'  3  :  > 

Upon  examination  of  the  Individual  contr ibutiens  to  total  predicted 
eroclon  rate ,  (each  is  plotted  against  distance  in  Figure  8'2),  the  contribution 
of  Impingement  heat  flux,  "q  again  is  seer  tc  be  negligible  The  question 
naturally  arises  whether  a  mBre  reasonable  fit  of  measured  erosion  ratt s  might 
be  obtained  by  assigning  a  value  much  greater  than  unity  to  the  “particle  energy 
accommodation  factor/'  a,  cr,  in  other  words,  by  taking  mechanical  erosion  to 
be  an  Important  result  of  particle  impingement.  It  can  be  seen,  however^  that 
the  shape  of  the  experimental  curve  would  not  be  matched  even  by  taking  “u,”  to 
be  as  high  as  10;  for  example,  because  the  peak  in  particle  flux  (Figure  S*]) 
occurs  downstream  of  the  last  point  (37.4  inches)  at  which  a  measured  erosion 
rate  value  could  be  determined,  (Beyond  this  point,  the  insula  ion  had  burn.-d 
back  to  the  graphite  backing  during  the  courr^e  of  the  firing) 

In  Figure  8-3 v  measured  and  predicted  erosion  rates  for  the  100  FW2 
aft  closure  are  plotted  versus  local  Mach  number  for  comparison  with  similar 
data  reported  by  Aerojet-General  Corporation®  2  for  the  100  FW2  aft  closure 
Correlation  versus  Mach  number  in  the  region  of  interest  appears  only  to  change 
the  relative  shapes  of  the  curves,  and  does  not  suggest  any  new  or  better 
empirical  approach  to  the  prediction  of  erosion  rates. 


8,4  FUTURE  WORK 

In  conclusion,  a  quite  sophisticated  eiid  detailed  analytical  procedure 
for  studying  nozzle  environment  effects  in  a  giv<.n  nozzle  has  been  demonstrated 
However^  a  real  test  of  the  validity  and  utility  of  the  method  has  nut  been 
completed.  A  larger  number  of  designs,  varying  in  scale  and  geometry^  should  be 
similarly  analyzed  to  develop  more  correspondence  between  the  real  and  predicted 
environmental  effects.  This  v7ould  lend  itself  to  yl  a  better  evaluation  of  the 
particle  Impingement  ''altdity^  (2)  increatjod  realiability  in  the  predictions 
made  by  such  an  analysts,  and  (3)  increased  ability  to  indicate  a  consistent 
method  of  applying  the  analysis  to  various  designs. 
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FIGURE  8-3.  COMPARISON  OF  EROSION  RATES  CORRELATED  AGAINST  LOCAL  MACH  NUMBER 
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SECTION  9 


SUMMARY  AND  RECOMMENDATIONS 


There  are  certain  limitations  and  uncertvtlnties  esaoclatad  with  aome 
aspects  of  the  work.  These  need  to  be  outlined  clearly.  A  straight for^^ard  evalu¬ 
ation  of  the  work  is  Ui^eded  for  Intelligent  planning  of  future  work.  Therefore »  a 
concise  evaluation  is  presented  in  the  form  of  this  sumDary  with  recoimendatlons . 

A  very  useful  particle  sampler  has  beer,  developed.  Small  enough  saraples 
can  be  taken  to  avoid  agglomeration,  and  to  display  all  the  particles  in  proper 
number  distribution  on  .ne  slide.  The  sample  is  not  subject  to  fractionation  by 
air  currents  or  variable  particle  sedlc^ntar ion  rates  such  as  might  occur  in 
simple  exposed  slide  sampling.  This  tool  can  be  applied  ro  motors  of  various  sizes 
and  propellant  formulations,  and  to  sampling  various  portions  of  the  nozzle  effluent 
cross  section.  This  would  furnish  much  needed  data  fur  nozzle  performance  analysis 
where  particle  size  Is  a  factor. 

Further  tests  with  the  sampler  would  be  worthwhile.  Very  large  particles 
might  tend  to  bounce  from  the  glass  slide  upon  Impact.  While  this  does  not  appear 
to  be  the  case  for  the  tests  conducted,  this  factor  should  not  be  overlooked.  If 
this  is  the  case,  multistage  sampling  could  be  used  to  allow  lower  collection 
ve loc it ies  for  the  1 sr ^e  particles 

The  cold  flow  particle  impingement  modeling  technique  nas  a  highly  so¬ 
phisticated  analytical  basis.  The  results  are  consistent  and  reproducible.  More 
varied  designs  than  those  modeled  in  this  program  shculd  be  studied  to  demonstrate 
the  effectiveness  of  the  technique.  Cases  having  more  severe  particle  impingement 
would  be  interesting.  Since  the  cold  flow  modeling  conditions  typically  do  not 
model  the  boundary  layer  thickness  studies  in  which  the  boundary  layer  is  varied 
should  be  conducted.  This  same  technique  of  analysis  and  facility  design  would 
make  a  very  interesting  study  using  dyed  droplets  impinging  or.  blotter  surfaces. 

The  technique  can  be  applied  to  a  study  of  particle  Impingement  through  shock  waves 
and  to  the  detection  of  flew  separation  regions  In  supersonic  flow.  Correct  particle 
size  distribution  data  for  a  given  motor  is  necessary  for  accuracy  in  cold  flow 
modeling  studies. 
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A  coid  flow  modeling  technique  has  been  developed  and  demonstrated  lor 
determining  nozzle  wall  pressure  contours.  By  using  this  data  In  convective  boat 
transfer  analysis^  the  nozzle  material  heat  load  severity  Is  Indicated.  The  severe 
heat  fluK  problems  of  one  particular  design  ha>^e  been  revcjled. 


3200°F. 


EmissLviCy  values  for  alumina  have  been  estabiisiiea  up  to  a  Lemperacure  ot 
The  data  are  applied  to  the  calculation  of  radiative  heat  flux  value^  to 
the  nozzle  wall  for  various  motor  sizes »  pressures,  and  alumina  loadings.  However, 
the  calculation  Is  only  approximate, since  many  assumptions  were  required  to  handle 
the  complexities  in  geometry,  particle  flow  field,  and  thermal  gradients  through 
the  field  of  view  in  the  nozzle.  The  values  calculated  probably  represent  a  lower 
limit  of  radiation,  since  unburned  carbon  and  oropeliant  additives  are  seen  to  add 
significantly  to  the  radiation. 


WhPe  thi*’  work  wee  conducted  mainly  with  a  hydr ogen-oxygen-alumina  flame, 
similar  studies  are  needed  on  actual  propellant  or  simulated  propell^uv  firings. 

More  work  is  needed  on  chamber  window  viewing  technique.  Studies  must  be  made,  as 
in  this  work,  on  flames  with  known  particle  size  distr.* butions,  If  scaling  laws  are 
to  be  inferred.  Particle  size  distribution  in  the  laige  motors  mutt  be  known  if  the 
scaling  laws  are  to  be  applied  to  them. 


Application  of  the  work  of  this  program  to  the  detailed  analysis  cf  heat 
fluxes  to  two  large  nozzle  designs  has  been  demonstrated.  This  allows  conclusions 
to  be  drawn  on  the  relative  importance  of  the  various  heat  "lux  modes,  and  on  the 
Implications  of  various  geometries  for  the  severity  ot  nozzle  attack.  However,  th^ 
study  of  more  cases,  and  the  input  of  the  nozzle  wa-1  reaction  produces,  is  needed  to 
demonstrate  more  ccvivinc ingly  the  utility  of  the  method  in  evaluating  proposed 
designs . 


A  technique  has  been  demonstrated  for  directly  measuring  pa..liclc  impinge¬ 
ment  heat  fluxes.  The  experiment  involved  a  specialized  case  where  a  relatively 
cool  inert  wall  was  used.  The  next  step  is  to  study  particle  impingement  effects  in 
a  similar  manner,  but  with  a  hotter  nozzle  wall  and  one  which  reacts  cc  tl  3  en¬ 
vironment  by  ablating  or  eroding. 
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APPENDIX  A 


PLUME  USING  A  1-INCH  KOCKE 


«  A  U 

I  r\  - 

CEV  rtDTOR 


^2^3 


A.l  INTRODUCTION 

The  prcLCi-''''  *>  i  »»*ni  num  oxide  paj'ticles  In  the  exhaust  products  of 

rocked  motors  warrants  detailed  study  In  order  (I)  to  permit  a  better  under¬ 
standing  of  thrust  losses  due  to  lack  of  equilibrium  between  particles  and  gas, 
(2)  tc  gain  information  on  effect  of  particle  addition  or.  radiant  heat  transfer 
to  rocket  nozzles,  and  (3)  to  understand  the  radiation  emitted  by  the  plume 
f-ciii  the  standpoint  of  plume  detection. 


The  radiant  properties  of  a  cloud  of  particles  whose  sizes  are  com¬ 
parable  to  the  wsvelength  of  light  can  in  principle  be  determined  by  (1)  deter 
mining  from  Mle  theory  the  particle  emission,  absorption,  and  scattering  cross 
sections  as  a  function  of  particle  size,  wavelength  of  light,  and  particle 
material  electrical  properties, which  presumably  are  temperature  dependent; 

(2)  solving  a  radiative  transfer  problem  using  the  above  particle  crosj  sec¬ 
tions  for  the  particular  geo»netry  involved.  (If  there  exists  interaction  De¬ 
tween  the  radiation  and  the  gas  dynamics  then  the  radiative  tvansfer  problem 
must  be  solved  simultaneously  with  the  gas  dynamics,  but  his  additional  cotnpll* 
cation  is  set  aside. > 


The  calculation  of  the  radiant  properties  of  rocket  plume  heavily 
laden  with  aluminum-oxide  particles  is  hampered  by  t'  ^  lack  of  knowledge  of 


^  V.  ^  ^  1  A  ^  M  V.  J  ^  M  ...  .  *  M  ^  .  m  ..  .J  jM  n  ...  ^  .r.  £  ^  ^ 

wtic  W«»  V  c  it  b.  U  UA  C  VlCpCllUdl^CT  \/ A. 


.  1  ^  4  ..  ..  1  ^  4  ..  .. 

A.  *  WCA  A  WA 


AI2O3.  A  major  aim  of  experimental  resear  h  in  this  area  is  che  determination 
of  particle  emission  (absorpt Ion)  and  scattering  cross  sections.  The  determi¬ 
nation  of  this  information  as  a  function  wavelength,  over  all  wavelengths, 
is  a  large  order,  A  second  aim  of  exper itr*enta  1  work  is  to  determine  the  integrated 
effects  of  the  particle  emission  by  determining  the  total  radiition.  This  is  of 
particular  importance  when  the  aim  of  the  work  Is  the  prediction  of  the  radiant 
heat  transfer  of  a  given  rocket  motor  cioud  Total  radiation  (i.e.,  integrated 
with  respect  to  wavelength)  experiments  generally  do  not  give  enough  iaformation 
to  permit  an  understanding  of  detailed  processes  of  the  radiant  emission. 


The  following  discussion  describes  a  method  of  experimentaliy  deter¬ 
mining  the  particle  emission  and  scattering  cross  sections  anil  describes  the 
results  of  attempting  such  experiments. 
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A. 2  experimental  MEAf.UR KMP:nT  OF  OPTICAL  PROPERTIES  OF  A  CLOUD  OF 
SCATTER  INC -EMITTING  PARTlClJiS 


The  complete  radiative  transfer  problem  for  a  given  cloud  of  particles 
of  arbitrary  refractive  Index  generally  will  be  awkward  to  solve.  Furthermore, 


c  o  i>a5  0  L  i  euuC  Lion  oi  i  i  mcrri  i  <i  i 


p‘  ^  A 


radlativ^’  transfer  probleir.  that  Is,  ttaelf,  unverified  by  independent  experiment. 
For  these  reasons,  It  Is  desirable  to  design  the  experiments  In  such  a  way  that 
the  results  can  be  Interpreted  by  recourse  to  the  solution  of  only  a  trivial 
radiative  transfer  problem.  By  performing  tests  at  small  optical  depths  for 
scattering  it  Is  possible  to  obtain  results  Chat  can  be  directly  Interpreted. 

In  the  equations  given  below,  the  presence  of  scattering  by  particles  (l.e.,  a 
change  of  direction  of  propagation  of  Incideni  light)  will  be  accounted  for, 
insofar  as  it  results  in  an  attenuation  of  the  directly  emitted  light  from  a 
given  point  In  space.  However,  the  additional  component  of  radiation  that  re¬ 
sults  by  scatteiing  at  the  point  of  radiation  originating  at  other  locations  is 
not  included,  and  the  results  can  be  expected  to  hold  only  when  multiple  scatter¬ 
ing  is  not  large.  For  this  reason  it  Is  desirable  to  perform  tests  at  various 
optical  depths  in  order  to  observe  the  dependence  of  the  results  on  particle 
concentration.  If  results  are  found  to  be  concentration  dependent,  then  only 
th*^se  obtained  by  extrapolating  to  zero  concentration  can  be  considered  valid. 


Radiant  Intensity  Due  to  Gas  and  Particles  at  Unequal  Temperatures 


The  radiation  generated  at  a  volume  element  of  unit  cross  section 

and  thickness,  dx ,  is  due  to  a  gas  of  temperature,  Tg ,  having  a  volumetric 

emission  cross  section,  and  a  suspension  of  particles  having  a 

volumetric  absorption  (eniTsslon)  cross  section  of  y_„  (cm*l)  is 

pa 


7„  R  (T  .  X)  + 


y  R 

pa 


(T  ,  X  ) 


dx 


(A-l) 


where 

R  (T,  X)  ■  Spectral  radiancy  at  wavelength  X  due  to  a 
blackbody  at  temperature  T. 

This  radiation  will  suffer  attenuation  due  to  the  self  absorption  by 
the  gas,  absorption  by  the  particles,  and  scattering  by  the  particles.  The 
total  attenuation  In  traversing  the  distance  x  from  the  interior  tc  the  surface 
of  the  medium  at  x  ■  0  is 


«xp  (-  y^x) 


(A.2) 
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wh*re 


>t 


(A-3) 


and 


7 


P® 


•caCtaring  cross  section  of  the  particles  per  unit 
volume  (cm* ^ ) . 


Integrating  the  product  of  expressions  (A-l)  and  (A*2)  over  the  length  zero 

to  L  gives  the  following  results  for  the  staradlancy  of  the  gas-particle  mixture; 


Rp  (  X  > 


(l-e  )  R  (T  .  X)  +  (l-«  "  )  R  (T  X) 

P  #  ^  a 


(A.4) 


Equation  (A-4)  applies  at  a  wavelength  where  the  gas  emits  strongly, as  in  the 
case  of  the  sodium  D  line  emission  when  a  sodium  compound  Is  added,  or  at  other 
wavelengths  where  the  gaseous  species  emit  characteristic  radiation. 


At  a  spectral  region  where  no  gaseous  emission  occurs,  Lquaclon  (A-4) 

reduces  to 


Rp  (X) 


P* 


R 


(Tp.X) 


l-e 


•|^l 


P« 


(A-5) 


The  question  arises  as  to  whether  the  particle  temperature  can  be  determined  by 
observation  of  the  radiation  emitted  by  Che  particle  cloud.  If  a  transmission 
test  is  made  Co  determine 


e 


)  L 


and  the  spectral  steradlancy  Is  measured,  the  product 


Les- 


-  R  (T  ,  X ) 

7+7  P 

'pa  '  ps 

la  deCarminad .  However,  further  infornation  on  the  ratio  of 


Je2. 


7  +7 

'pa  'pa 


la  neceaaary  In  order  to  determine  particle  temperature,  and  thia  information 
la  uaually  not  available. 
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If  carbon  is  present  and  Is  in  thermal  equilibrium  with  the  gas,  as  would 

be  expected  due  to  the  small  sireo  of  carbon  parClcies,  then  Equation  (A-4) 

applies  if  7  is  replaced  by  7^  +  7^  where  7  is  the  absorption  cross  section 
per  unit  volu^  due  to  carbon.  iThe  scattering  by  carbon  particles  is  generally 
small  compared  to  absorption;  see  references  A. I  and  A. 2.)  carbon  Is  present, 

Che  determination  of  the  radiation  solely  due  to  MgC  or  AI2O3  particles  is  not 
possible  without  maklrig  tenuous  assuinpt  Ions  •  If  the  radiation  due  to  particles 
other  chan  carbon  is  of  Interest,  then  one  should  eliminate  sources  of  radiation 
of  comparable  magnitude. 

A  direct  method  of  learning  the  radiation  properties  of  a  particle 

cloud  is  to  perform  a  total  radiation  test  by  means  of  a  spectrally  neutral 

radiation  transducer  such  as  a  thermopile.  In  this  type  of  experiment,  the 
total  radiation  is  measured  and  can  be  directly  related  to  particle  cloud  emls- 
sivlty  if  particle  temperature  is  know >  The  application  of  this  type  of  infor¬ 
mation  to  solve  the  large  scale  rocket  motor  heat  transfer  problem  may  or  may 
not  offer  difficulty. 

A. 3  EXPERIMENTS  ON  SMALL  ROCKET  ENGINE 


A. 3.1  EXPERIMENTAL  PROCEDURE 

Experiments  were  performed  on  a  100  Ib-ft  thrust  engine  using  gaseous 
hydrogen  and  oxygen  as  propellants  and  an  aqueous  slurry  of  aluminum  oxide 
particles.  ^  description  of  the  rocket  motor  which  was  designed  with  a  heat 
sink  copper  combustion  chamber  and  a  nozzle  of  porous  tungsten  impregnated  with 
silver,  la  given  in  Table  A.l.  The  e.igine  was  found  to  operate  at  all  desired 
0/F  ratios  without  tendency  toward  unstable  combustion. 

TABfi:  A.l 


Motor  I.D. 

1  i  ch 

Motor  Length 

10  inches 

Throat  Diameter 

0.300  inch 

Exit  Diameter 

0.7  inch 

Exit  Area  Ratio 

6.02 

Contraction  Ratio 

11.1 

L* 

111  inches 

The  propellant  flow  rate  control  and  flow  rate  measurements  were 
achieved  by  use  of  calibrated  critical  orifices.  Pressures  were  measured  up¬ 
stream  and  downstream  of  the  critical  orifices  by  means  of  ressure  transducers, 
and  gas  cemperaturee  were  measured  by  welding  thermocouples  to  the  outside  of 
Che  propellant  feed  lines  upstream  of  orifices  and  wrapping  lines  with  insulation. 
This  method  of  measuring  gas  temperature  is  not  fully  satisfactory,  since  the 
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tubing  therms  I  response  tlmo  is  somewhat  long,  but  conventional  thermocouples 
protected  in  wells  are  probably  almost  as  slow  and  were  found  to  be  a  source 
of  excessive  leakage 

The  slurry  consisted  of  a  mixture  of  60  percent  by  weight  aluminum 
oxide,  40  percent  by  weighL  of  water  and  was  fed  into  the  combustion  chamber  by 
means  of  n  pneumatic  piston.  The  aluminum  oxide  used  in  the  slurry  was  a 
grinding  compound  (Norton  Co  Type  I-C  metallurgical  alumina)  consisting  of 
99,98  percent  y  -AI2O3  of  nominal  1  micron  size.  Particles  were  rounded  but 
nonspherical  and  were  far  more  uniform  in  size  than  most  commercial  grades  of 
alumina  powder  It  is  believed  that  the  use  of  the  metallurgical  grade  alumina 
provided  a  much  greater  control  over  injected  particle  size  than  would  other* 
wise  have  been  possible  The  particles  were  received  in  five  separate  packages 
which  were  blended  together  by  stirring  overnight  in  a  paddle-type  mixer. 
Particles  as  received,  and  also  after  mixing, were  highly  agglomerated  but  could 
be  de-agglomerated  in  small  quantities  by  pressing  between  a  microscope  slide 
and  cover  glass  if  water  or  oil  was  used  as  an  Immersion  medium.  The  unaggloro* 
erated  particles  display  Brownian  motion  when  suspended  in  water  and  observed 
in  a  high  power  (970X)  optical  microscope  with  an  oil  immersion  objective. 

No  attempt  was  made  initially  to  de-agglomerate  the  particles  prior  to  injec¬ 
tion  into  the  rocket  as  it  was  believed  that  combustion  and  shear  due  to  tur¬ 
bulence  in  the  combustion  chamber  might  accomplish  this  More  detail  on  particle 
processing  eventually  used  is  given  below.  Additional  instrumentation  used  on 
the  tests  were  the  usual  spectral  Intensity  measurements  at  X  ■  0,59/i,  frans- 
mission  tests  in  the  same  wavelength  interval,  and  total  radiation  tests.  Total 
radiation  measurements  were  made  by  means  of  a  commercial  rapid  response  thermo¬ 
couple  with  a  KBr  window  that  passes  all  radiation  in  the  0  4/i  to  25M  interval. 
This  device  is  spectrally  neutral  and  receives  more  than  99,5  percent  of  all 
energy  emitted  by  a  blackbody  at  2500^K,  thus  permitting  a  total  or  integrated 
(with  respect  to  wavelength)  radiation  measurement  to  be  made. 

Transmission-emission  measurements  were  initially  made  at  a  position 
bmin  downstream  of  the  exit  plane  of  the  rocket  nozzle  where  the  propellant 
gases  were  expanded  approximately  to  atmospheric  pressure.  The  first  few 
firings  of  the  engine  revealed  the  usual  minor  defects  which  were  corrected 
without  serious  delay.  Initial  firings  v/ere  made  at  an  0/F  of  about  2.8  in 
order  to  maintain  the  chamber  temperature  low  enough  to  prevent  melting  of  the 
particles, thereby  assuring  good  size  control;  i.e.,  sizes  all  equal  to  original 
sizes  of  injected  particles  Under  these  conditions  (Runs  63-2-1  and  J) ,  no 
particle  radiation  or  gas  radiation  was  detected,  and  low  0/F  firings  were  dis¬ 
continued.  High  0/F  firings  (K,  L,  and  M)  revealed  more  interesting  but  un¬ 
expected  results.  When  particles  were  first  injected,  no  immediate  increase  in 
radiation  was  detected  by  the  total  radiation  detector,  but  as  the  slurry  mix 
was  continued  a  gradually  increasing  signal  was  recorded.  This  signal  corre¬ 
sponded  to  the  increase  in  a  shower  of  large  **sparks"  observed  visually  in  the 
plume.  Examination  of  the  nozzle  revealed  heavy  deposition  on  entrance  and 
exit  regions.  It  was  concluded  that  the  shower  of  sparks  was  due  to  the  shed¬ 
ding  of  large  droplets  by  the  shearing  action  of  propellant  gases.  A  second 
observation  was  that  the  transmission  of  the  plunie  as  measured  by  the  optical 
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system  was  nearly  100  percent  even  with  high  particle  flow  rates.  It  was  con¬ 
cluded  that  the  rocket  was  not  *uncrionxig  as  a  slmulatcr  perhaps  because  the 
particles  were  agglomerated  and  were  being  deposited  the  nozzle  in  excessive 
amounts  due  to  large  velocity. 

The  alumina  slurry  was  then  subjected  to  further  processing  in  crder 
to  achieve  de-agglomerat  ion .  It  was  found  that  the  use  of  1/2  perc^"*nt  by  weight 
of  trlsodlum  phosphate  as  a  wetting  agent  followed  by  overnight  mixing 
on  a  ball  mill  [25  percent  volume  of  alundum  balls,  abcut  I  cm  in  diameter, 
added)  achieved  a  high  degree  c f de-agglomeration  of  the  individual  particles. 

The  slurry  mixture  was  vastly  more  fluid  In  consistency  and  ecu  Id  probably  have 
been  further  enriched  with  AI2O3.  Some  tendency  toward  settling  of  the  particles 
remained,  and  although  this  was  not  considered  serious,  the  par  ides  were 
loaded  into  the  piston  immediately  before  firing  in  order  to  minimize  settling. 

The  use  of  a  noncarbonaceous  gelling  compound  to  prevent  settling  would  be  an 
improvement  for  the  slurry  in  f ut ure  tests . 

The  de-agglomerated  particles  were  then  fired  in  the  rocket  (Runs  N, 

0,  P  and  Q)  ,  and  it  was  observed  that  the  shower  of  sparks  was  reduced  but  still 
present  t  that  transmission  remained  very  high  ('v  96  percent),  and  that  simula¬ 
tion  of  rocket  motor  firing  was  not  achieved.  In  a  discussion  that  ensued  to 
determine  the  course  of  action  to  be  pursued,  a  work  was  cited  in  which  a  small 
solid  propellant  rocket  was  fired,  and  the  investigator  found  that  extensive 
nozzle  deposits  were  formed  with  attendant  shedding  of  large,  droplet-sized, 
particles.  It  was  decided  that  further  firing  of  the  BRACE  engine  per  se  would 
be  of  no  value  because  of  its  very  small  size,  and  that  attempts  should  be  made 
to  fire  the  motor  without  a  nozzle  in  the  hope  of  eliminating  droplet  shedding. 

The  firing  sequence  was  modified  to  provide  a  delay  before  slurry  injection  of 
about  1  to  2  seconds.  Optical  measurements  were  made  on  all  succeeding  runs  at 
a  position  14  mm  downstream  of  the  exit  plane  of  the  throat  less  rocket  motor. 

The  first  firing  (Run  R)  conducted  after  these  modifications  were  made 
showed  that  combustion  could  be  sustained  without  the  nozzle  although  the 
photomultiplier  in  the  monochromator  was  overdriven  and  no  transmission  measure¬ 
ment  was  :^btained.  Runs  S  and  T  were  successful,  and  transmission  was  observed 
to  decrease  with  increasing  particle  weight  fraction--no  sparks  were  observed. 

A  series  of  runs  (Runs  U  to  Z)  were  conducted  very  close  to  stoichiometric  0/F 
in  which  gas  and  gas  plus  particle  radiation  were  measured  as  particle  fraction, 

0  ,  %#a8  varied.  A  second  series  of  runs  (Runs  AA,  BB,  CC)  were  also  conducted 
at  higher  total  mass  flows  with  C/F  as  before  and  0  varied.  In  all  cases  the 
radiation  was  only  perceptibly  increased  (over  the  level  due  to  the  gases  (mainly 
H2O  only)  when  the  slurry  was  Injected  into  the  combustion  chamber. 

Alumina  particles  ejected  by  the  throatless  rocket  were  captured  directly 
on  a  greased  slide,  and  also  on  plain  slides,  by  the  cascade  impactor  for  a  number 
of  runs.  Examination  of  these  slides  showed  chat  all  captured  particles  were 
spherical  with  many  agglomerates  formed  on  the  slides  consisting  of  very  small 
particles  and  also  some  larger  single  spherical  particles  ranging  up  to  about 
10  microns.  The  fact  that  the  particles  recovered  after  firing  were  spherical 


indicates  that  the  alumina  had  been  heated  above  its  melting  point  (2320Pk), 

but  no  detailed  information  on  size  distribution  of  recovered  particles  is  possible, 

becarse  of  the  extreme  difficulty  of  collecting  a  truly  representative  sample. 

(It  was  also  observed  that  the  de -agglomerated  particles,  if  stored  for  four  days 
in  an  aqueous  suspension,  would  spontaneously  re-agglomerate.  It  is  therefore 
conceivable  that  some  of  the  large  particles  recovered  on  the  slides  viere  results 
of  re-agglomeration  in  the  slurry  prior  to  firing,  but  this  is  improbable  since 
the  slurry  was  used  very  soon  (four  hours)  after  being  removed  from  the  ball  mill. 

A  number  of  temperature  measurements  by  observation  of  sodium  D  line 
with  and  without  a  background  source  of  known  temperature  were  made  during  the 
first  portion  of  the  firing  when  only  the  H2-O2  were  flowing.  At  this  portion 
of  the  firing,  only  a  trace  amount  of  sodium  is  present,  and  the  concentration  of 
sodium  is  likely  to  vary  between  the  time  that  the  gas  (Ma  D  Line)  emission  is 
observed  and  the  gas  plus  lamp  emission  is  observed.  The  temperatures  so  measured 
are,  therefore,  not  considered  accurate  and  varied  widely  from  2400^  to  3000^ 
about  a  mean  value  of  2680®K  even  though  0/F  was  held  constant.  These  measure¬ 
ments  are  considered  the  best  available  estimates  of  the  gas  temperature  in  the 
ibsence  of  the  dilution  effect  of  the  slurry.  The  raw  data  has  been  reduced  to 
an  intermediate  form  and  is  presented  in  Figures  A-1  and  A-2  and  Table  A. 2. 

Figure  A-1  is  a  plot  of  the  plume  transmission  versus  particle  weight  fraction. 

The  data  shows  a  vide  scatter  around  a  straight  line  which  corresponds  to  a  mean 
diameter  of  0.8  micron.  Transmission  tests  were  also  conducted  on  the  de* 
agglomerlzed  particles  in  an  aqueous  suspension  with  the  same  equipment  that  is 
used  on  the  rocket  motor.  This  data,  shown  in  Figure  A*2 ,  shows  a  curvature  on 
a  log-linear  plot  which  may  be  due  to  forward  scattering  being  received  by  the 
photocell.  As  concentration  is  increased,  the  forward  scattering  increases  to 
cause  a  departure  from  the  usual  exponented  transmission  Law.  This  type  of 
error  is  discussed  in  Reference  A. 3  and  can  be  corrected  for  by  plotting 
[specific  turbidity  or  7  ■  1/L  ln(l/T)  divided  by  concentration]  versus  concen¬ 
tration,  Figure  A-3.  The  value  of  specific  turbidity  at  zero  concentration  is 
free  from  error  due  to  forward  scattering  and  can  be  used  to  determine  ossan 
size.  The  result  obtained  was  0.5  micron  which  appears  to  be  a  reasonable  value. 


A. 3. 2  DATA  REDUCTION 

a .  Particle  Scattering  of  Gas  Emitted  Radiation 

A  brief  study  was  made  to  determine  if  the  increased  radiation  measured 
while  the  particles  were  flowing  was  merely  due  to  a  change  in  the  angular  dis¬ 
tribution  of  the  radiation  emitted  by  the  gaseous  species.  Since  the  important 
radiating  specie  is  H2O,  the  emission  bands  are  Mvelengths  of  l.I,  1.38, 

1.87,  2.7,  6.3  and  20  microns.  According  to  Penner,^*^  the  last  three  of  these 
bands  account  for  95  percent  of  the  total  intensity  of  the  H2O  emission.  It 
was  estimated,  based  on  the  fact  that  the  transmission  was  42  percent  for  a 
particle  fraction  of  0.365  (Run  Z),  that  under  such  circumstances  the  transmission 
would  Increase  to  83  percent  for  the  2  ,lfi  H2O  emission  band.  Most  runs  were  con¬ 
ducted  at  lower  particle  fractions,  and  transmission  would  be  even  higher.  On  this 
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bASlS  ,  one  cen  9€y  that  the  optical  depth  for  scattering  Is  small  for  all  Impor¬ 
tant  H2O  bands.  Therefore,  the  possibility  that  the  increased  radiatlun  measured 
in  the  normal  direction  is  due  to  scattering  by  the  part ic lea,  causing  a  change  of 
angular  distribution  of  energy  emitted  by  H2O  with  no  net  change  in  total  heat 
transfer,  is  unlikely.  The  increase  in  the  total  radiation  signal  with  the  slurry 
flowing  is  attributed  to  particles  emitting.  This  premise  is  further  confirmed 
by  the  fact  that  at  the  visible  wavelength  of  0.59  micron,  where  no  radiation  was 
obaarved  without  the  slurry,  a  detectable  continuum  was  noted  only  when  the 
particles  are  present. 

b ,  Particle  Cloud  Emissivity 

Since  the  optical  depth  fur  scattering  of  H2O  radiation  by  the  particles 
is  small,  the  particle  cloud  emissivity  can  be  determined  by  subtracting  the 
corrected  staradiancy  ^2^  from  the  measured  steradiancy  due 

to  the  perticlea  plus  gas,  Kpg  (T2).  The  particle  cloud  normal  emissivity  was 
then  found  by  dividing  by  the  st«^radiancy  due  to  the  blackbody  the  tempera¬ 
ture,  T2,  that  was  found  experimentally  with  the  slurry  flowing: 


c 


p 


R 


(T^)  -  (Xp 


(A.S) 


The  exparicnentsl  results  reduced  in  this  manner  are  plotted  in  Figure  A-4. 

The  data  falls  on  two  lines  corresponding  to  the  different  gas  flow 
ratas  of  0.16  to  0.30  Ib/sec  for  runs  AA,  B8,  and  CC .  Each  of  theae  two  series 
of  the  data  ahova  moderate  scatter,  but  the  best  guess  seems  to  be  that  the  data 
falls  on  straight  lines  corresponding  to  the  two  flow  rate^.  Several  unsatis¬ 
factory  explanations  can  be  offered  for  the  variation  of  plume  emissivitles  at 
the  two  different  flow  rates.  One  of  these  is  the  following.  At  the  higher 
flow  rates,  the  flow  is  choked  at  the  exit  plane  of  the  motor  dne  to  heat  addi¬ 
tion,  and  the  exit  plane  pressure  is  higher  than  atmospheric.  The  gases  then 
expand  a  diverging,  supersonic  jet  upon  leaving  the  nozzle  resulting  in  a 
dro  in  the  temperature  (T2)  and  reduced  H2O  emission  measured  in  runs  AA,  BE, 

CC  «ven  chough  0/F  ratio  Is  held  very  nearly  constant  (see  Table  A. 2;  sodium 
lines  %^re  in  reversal  only  for  the  last  three  runs,  lamp  current  held  constant). 
Particle  cloud  emissivity  measured  high  presumably  because  of  the  greater  depth 
and  oisas  flow  particles.  This  picture  would  be  expected  to  cause  a  reduced 
plume  transmission,  but  this  result  was  not  obtained.  It  is  this  latter  point 
Cliat  jeems  to  br  difficult  to  reconcile;  while  the  particle  stetc-diancy  increased 
on  the  last  three  runs  where  lower  temperatures  were  measured,  the  plume  trans¬ 
mission  shows  no  difference  for  equivalent  value  of  0  than  on  the  prevlou^^  six 
runs . 


In  any  event,  the  data  obtained  from  the  earlier  six  runs  are  con¬ 
sidered  more  self-consistent,  and  for  chat  reason,  to  be  the  more  reliable.  The 
difficulty  in  interpreting  the  results  of  the  last  three  data  points  is  attri¬ 
buted  to  an  uncertainty  of  particle  velocity,  cloud  thickness,  and  hence  particle 
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concent  rat  Iona ,  Possible  existence  of  similar  uncertainties  in  the  results  of 
Runs  U  to  >1  cannot  overruled.  Such  uncertainty  can  be  eliminated  In  future 
testa  by  per  forr  Ing  a  measurement  of  static  pressure  Just  upstream  of  the  exit 
plane«and  adjust  ing  the  /low  rate  so  that  atiTK>8pher  Ic  pressure  Is  not  exceeded 
at  this  pc 1 n t 

A. 4  INTLGRATED  RADIATION  KEASUREMENTS 

Measurements  of  the  Integrated  radiation  ( i  .e . ,  the  radiation  emitted 
over  all  wav*^lengths)  emitted  by  the  particles  were  made  by  means  of  a  fast  re¬ 
sponse  t iiermocouple  with  KBr  window  ( Per ki n-E imer  thermocouple  and  standard 
preamplifier).  Radiation  emitted  by  the  rocket  plume  was  collected  by  a  spherical 
cn-axls  mirror  and  brought  to  a  focus  at  the  detector,  which  consisted  of  a  blackened 
gold  atrip  with  a  fine  thermocouple  attached.  The  photosensitive  area  of  the 
transducer  was  0.2  x  2.0  mm, and  the  entire  window  sire  was  about  4«0  mm  in  diameter. 
The  mirror  and  detector  assembly  were  enclosed  in  an  aluminum  box  with  a  single 
opening  to  wnlch  an  iris  diaphragm  was  mounted.  The  box  was  purged  with  dry 
nitroge.i  to  prevent  atmospheric  TOioturc  from  damaging  the  hygroscopic  KBi 
window . 


The  fast  response  thermocouple  uses  a  mechanically  chopped  radiation 
signal  at  13  cps  in  order  to  generate  an  ac  output  which  allows  discrimination 
of  signal  at  a  dc  level  that  does  not  depend  on  ambient  temperature  and  other 
factors.  It  was  found  chat  the  output  of  the  thermocouple  was  relatively  sensi¬ 
tive  to  chopping  frequency  (e.g.,  a  reduction  of  from  13  to  ll  cps  increased 
signal  by  10  percent),  and  for  this  reason  a  strobctach  was  used  to  adjust  chop¬ 
ping  wheel  speed  ininedlatcly  before  cilibretlon  or  daily  prior  to  a  series  of 
tests.  The  thermocouple  signal  was  amplified  in  r he  preamplifier  and  amplifier 
and  was  recorded  on  the  oscilloscope  without  rectli icat ion  as  an  ac  signal. 

The  radiometer  was  calibrated  by  n^ans  of  a  vommercial  blackbody 
cavity  (Infrared  Industries  Model  404)  up  to  a  level  of  steradiancies  of 
4.83  watts  cm"2  stcr“^  corresponding  to  a  temperature  oi  1278^K.  In  calibrating 
the  radiometer,  the  front  face  of  the  blackbody  cavity  and  the  thermocouple  were 
us«d  A9  the  conjugate  focal  points.  in  perfornirig  measurements  on  the  rocket 
plume,  the  optics  were  arranged  so  that  a  point  on  ti.^  axis  of  the  motor  and 
the  thermocouple  were  conjugate  focal  points.  Both  the  image  of  the  blackbody 
cavity  and  the  image  of  the  plume  were  large  compared  to  the  target  area  of 
the  detector.  To  provide  a  check  on  the  stability  of  the  radiometer,  a  sub¬ 
standard  calibration  was  used  on  a  daily  basis  before  each  run.  This  was  accom¬ 
plished  by  deflecting  the  light  from  the  tungsten  ribbon  filament  lamp  into  the 
radiometer  by  means  of  a  front  surface  mlrroi  located  at  the  image  of  the  lamp 
on  the  rocket  axis.  It  was  believed  that  if  a  sharp  image  of  the  lamp  was  pro¬ 
jected  onto  the  transducer  Cher  the  day-to-day  alignment  of  the  mirror  would  be 
too  critical  to  provide  a  meaningful  check  of  the  stability  of  the  i adio-iiecer  , 

An  additional  lens  was  therefore  inserted  in  the  optical  system  in  order  to 
defocus  the  image  of  the  lamp  on  the  thermocouple  and  eliminate  the  need  of  a 
highly  critical  alignment  of  the  mirror.  The  irradiancy  at  the  thermoccaple 
with  the  defocusslng  lens  in  position  was  very  close  to  the  irradiancy  due  to 
the  rocket  plume  itself. 
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Erratic  performance  of  the  radiometer  was  noted  during  the  first  series 
cf  tests,  its  sensitivity  seeming  to  decrease  by  a  factor  of  four.  However, 
during  t’  c  next  three  days  of  testing,  the  sensitivity  remained  acceptably  con¬ 
stant;  thL'  substandard  check  was  0.35,  0.36,  0.39  Inches  of  deflection  for 
i6  amps  lamp  current.  After  this,  the  radiometer  was  recalibrated  with  the  black- 
body,  aid  since  Its  calibration  had  changed,  all  equipment  was  again  moved  back 
to  the  t^st  cell  for  a  check  on  the  substandard.  On  this  test,  a  deflection  of 
0.39  inch  at  16  amperes  lamp  current  was  recorded.  Although  the  cause  of  the 
lack  o:  ^  ability  was  never  ascertained,  because  It  occurred  early  In  the  testing 
and  was  fc llowed  by  a  week  of  stable  performance  and  because  the  useful  data 
were  taken  on  only  the  last  day  after  which  a  primary  standard  calibration  and 
good  substandard  checks  were  performed,  confidence  is  high  that  the  radiometer 
was  performing  well  during  tiie  important  tests 


Tlie  amp  L I f  led  s  Igna I  from  the  radiometer  was  found  r.o  be  1  inear  In 
steradlancy.  The  "herraocoup le  is  designed  for  use  in  a  dual  beam  system  where 
a  constant  check  on  Itt  sensitivity  is  provided  by  use  of  a  reference  signal 
and  ratio  recording.  In  the  present  case,  the  preamplifier  tube  filaments  and 
plates  were  supplied  by  external  power  supplies  which  could  be  a  source  of  If-ck 
of  repeatability.  The  geometry  of  the  optical  bench  did  not  allow  the  use  of 
a  blackbody  calibration  In  situ  as  would  have  been  desirable.  For  these  reasons, 
it  is  estimated  that  the  accuracy  of  the  measurement  was  of  the  order  of 
+5  percent . 


A. 5  COMPARISON  CF  TOTAL  RADIATION  RESULTS  WITH  SPECTRAL  ANALYSIS 
RESULTS 


With  the  knowledge  that  the  particle  emission  cross  section  Ic  inde¬ 
pendent  of  wavelength  and  hence  the  particles  emit  as  grey  bodies,  it  is  possible 
to  reduce  the  data  of  total  radiation  tests , (obtained  by  Dobbins  and  reported  in 
the  third  quarterly  report),  to  spectral  emission  cross  sectlcjn.  The  results 
of  this  reduction  indicate  an  emission  cross  section  fer  the  particles  used  in 
the  latter  tests  (D^ve"^  )  cn  the  order  of  10“^  .  It  is  not  possible  to 
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radiation  test  is  not  sensitive  when  significant  amounts  of  water  vapor  art 
present.  The  main  conclusions  from  the  total  radiation  Is  that  the  particle 
radiation  was  small  compared  to  the  gas  (H2O)  radiation  In  the  simulator  firing. 
The  raa^nitude  of  this  result  is,  however,  in  agreement  with  the  results  cf  the 
spectral  emission  tests  of  Carlson. 


A  number  of  improvements  can  be  made  in  order  to  give  more  accurate 
measurements  of  both  total  and  spectral  emissivity  tests  of  oxide  particles* 
Both  types  of  measurements  would  be  Improved  by  using  a  propellant  formulation 
that  eliminated  the  background  radiation  due  to  H2O >  This  can  be  done  by  se¬ 
lecting  d  propellant -ailuent  combination  that  produces  Jiatomic  species. 
Spectral  measurements  will  be  improved  by  using  a  spectroir.crcer  on  several 
monochrome ters  that  make  observations  at  fixed  wavelengths.  The  sensitivity 
can  then  be  improved.  Transmission  tests  can  be  made  in  oruer  to  deduce  pert- 
Icle  sire  and  scattering  cross  section  Thi3  information  can  then  be  used  Cc 
predict  emission  cross  section  as  function  of  particle  size 
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